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An  engineering  methodology  has  been  developed  for  the  prediction  of  pro- 
pulsive lift  system  Induced  aerodynamic  effects  for  multl-llf t-jet  VTOL  aircraft 
operating  in  the  hover  mode  in  and  out  of  ground  effect.  The  developed  raethod- 
ology  takes  into  account  the  effects  of  aircraft  geometry,  aircraft  orientation 
(pitch,  roll)  and  height  above  the  ground  plane,  cross-flow,  lift-jet  vector 
directions  with  respect  to  the  airframe,  lift-jet  exit  flow  conditions  and  nnn- 
clrcular  nozzle  exit  geometry  for  low  aspect  ratio  nozzles, 
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In  ground  effect,  the  prediction  methodology  proceeds  logically  from  the 
defined  aircraft  lift-jet  exits  through  the  free-jets,  jet  Impingement  points, 
wall-jets,  fountain  bases  (stagnation  lines),  and  fountain  impact  on  the  air- 
frame undersurface,  The  Induced  suckdown  flows  are  computed  from  the  potential 
flowfleld  induced  by  the  turbulent  entrainment  of  both  the  free-jets  and  the 
wall-jets  in  ground  effect  and  from  the  free-jets  alone  out  of  ground  effect. 
The  jet  and  ground  flowfleld  methodology  defines  necessary  fluid  flow  propertle; 
in  the  free-jets  and  wall-jets,  the  jet  impingement  regions  and  fountain  base 
regions,  and  the  fountain  upwash  for  certain  classes  of  multi-jet  fountains. 

The  entire  methodology  Is  dependent  on  .nircrait  and  flowfleld  Kuometry,  and  an 
effort  has  been  made  to  clarify  and  simplify  Hoometric  reciuirements. 

The  methodology  emphasizes  accurate  analytical  and  empirical  modeling  of 
free  and  wall-jet  flows  (especially  turbulent  entrainment)}  ground  wall-jet 
interactions  Including  the  computation  of  multl-jet  stagnation  lines  and  foun- 
tain base  flow  direction}  and  fountain  development  and  Impingement  on  the  air 
frame  under-eurface.  An  analogy  between  certain  classes  of  multl-jet  fountains 
and  eelf-inpinglng  paira  of  circular  jats  has  been  ueed  to  describe  fountain 
geometry,  and  to  define  experimental  test  configurations  which  were  employed  to 
empirically  Investigate  fountain  development  and  impingement  on  the  airframe 
under-surface.  Tachniquea  for  predicting  fountain  induced  forces  on  the  air- 
frame have  bean  developed.  An  investigation  of  the  effects  of  protuberancp-s 
waa  included  in  the  experimental  effort,  and  techniques  for  predicting  the 

Incremental  forces  due  to  the  Interaction  of  fountain  flows  with  protuberances 
were  developed. 

The  Induced  potential  flow  field  due  to  turbulent  Jet  entrainment  is 
computed  through  use  of  the  Douglas  Neumann  potential  flow  program.  The  prograr 
has  been  modified  for  direct  application  to  the  VTOL  flowfleld  computation 
problem,  A User  Manual  for  the  modified  Douglas  Neumann  program  has  been  pre- 
pared. ■ 

The  developed  methodology  represents  a state-of-the-art  approach  for 
the  engineering  prediction  of  ground  flowflclds,  induced  sackdown,  and 
fountain  formation  and  development  for  multl-jet  VTOL  aircraft. 
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An  analytical  program  has  been  conducted  for  the  Naval  Air  Development 
taaiter  (NADC)  by  the  McDonnell  Aircraft  Company  (M(;AU0  to  develop  n 
methodolot’y  for  tlie  prediction  of  propulaion  Induced  aerodynamic,  effects  for 
mu  1 1 1 - 1. 1 1't-Jet  VTOL  aircraft  operating  In  the  hover  mode  in  and  out  of 
ground  effect.  The  methodology  emphasizes  accurate  analytical  and  empirical 
modeling  of  free  and  wall-.let  entrainment,  ground  wall-jet  interaction  and 
subsequent  formation  of  fountains,  and  fountain  Impingement  on  the  undersur- 
face of  the  f uselage/wing.  Induced  flows  about  the  wlng/body  are  computed 
through  the  use  of  the  Douglas  Neumann  potential  flow  computer  program  which 
has  been  modified  for  application  to  the  VTOL  problem. 

These  analytical  techniques  and  empirical  correlations  have  bean  incor- 
porated into  the  resulting  methodology  for  force  and  moment  prediction i 
Tlie  basic  analytical  wo»"k  described  in  this  report  was  funded  under 
contract  //N62269-7fi-C-0086  to  MCAIR.  Moat  of  the  experimental  data  presented 
in  this  report  concerning  the  formation  and  development  of  fountains,  and  the 
data  pertinent  to  the  description  of  the  interaction  of  the  fountain  with  the 
fuselage/wlng  and  the  effect  of  protuberances  were  obtained  as  part  of  a MCAIR 
Independent  Research  and  Development  Program  during  the  calendar  year  1976. 

The  program  manager  of  this  study  was  Dr.  Donald  R.  Kotansky  and  the 
Principal  Investigator  was  Norbert  A.  Durando,  both  of  the  MCAIR  Aerodynamics 
Department.  The  technical  effort  associated  with  the  Douglas  Neumann  potential 
flow  computations  was  performed  by  Deavi  R.  Bristow  and  Philip  W.  Saunders 
of  MCAIR  Aerodynamics.  The  Naval  Air  Development  Center  Project  ’ingineer  was 
John  D.  Cyrus. 

The  authors  wish  to  express  their  apprec  la i; Ion  for  the  guidance  and 
technical  support  provided  by  Chester  W.  Miller,  Chief  Technology  l.nglneer  -• 
Aerodynamics  and  Robert  B.  Jenny,  Branch  Chief  Technology  - Aerodynamics, 

MCAIR. 
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Direction  cosines 

Integer  index  for  angular  measurement 

Momentum 

Jet  momentum 

Jet  momentum  at  nozzle  exit 
Mass  flow 

Total  Jet  mass  flow 

Jet  mass  flow  at  nozzle  exit 

Entrained  jet  mass  flow  ■ fti,.,,  - i^i, 

J i J 0 

Direction  cosines 

Normal  distance  above  ground  plane,  number  of  panels  in 
Douglas  Neumann  potential  flow  program 

Normal  distance  to  stagnation  line  l\i  ground  piano 

Normal  distance  to  U /2 
max 

Nozzle  proHSurc  ratio 
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■J  v/21 


T /T 
o'  « 
Je 

U 


X,  y,  z 


x-|^ I y I 
Ax,  Ay,  Az 


.1 

.] 
1 

Direction  coulnes 
Pre.HHure 

Radius  meaBiire;cl  from  real  or  apparent  jet  impingement 
stagnation  point,  radius  from  fountain  virtual  origin 

Ground  radius  at  edge  of  impingement  region 

Radial  correction  to  apparent  jet  impingement  stagnation 
point 

Jet  eitlt  Reynolds  number 

Radius  measured  along  upp  surface  or  wall 
Jet  radius 

Radius  of  free  jet  from  jet  centerline  to  location  of  1/2 
Jet  centerline  velocity 

Radius  to  1/2  jet  centerline  velocity  at  impingement  point 
Distance  between  Impinged  jets,  nozzle  centerline  spacing 
Nozzle  temperature  ratio 
WaU-jet  velocity 
Freustream  velocity 

Jet  exit  velocity.  Jet  velocity,  fountain  velocity 
Cartesian  coordinates 

Input  data  coordinates  for  Douglas  Neumann  program 
Coordinate  translatlono,  displacements 
Jet  centerline  coordinate 
Conical  fountain  semi-apex  angle 

Jet  impingement  angle  measured  witli  respect  to  ground  surface 

Wing  leading  edge  sweep  angle,  direction  in  ground  plane  of 
projection  of  negative  of  jet  main  velocity  vector  at  impinge- 
ment point 

Momentum  correction 

Flap  deflection  angle 
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C 

0 

0 


K 

i 

P 


y 

Subscripts 
1,  2,  3.  4 
A 
c 
e 
F 
1 
J 

max.  , m 
N , n 
o 
R 
s 

at 


Normalized  coordinate  in  jet  deflection  plane  (Fig.  ll-l) 
liluler  angle 

Stagnation  line  alope,  fountain  or  jet  azimuthal  ang].e 
measured  from  vertical  about  virtual  origin 

Fountain  velocity  or  momentum  directional  angle 

Normalized  coordinate  in  Jet  deflection  plana  (Fig.  Il-i) 

Density 

Jet  deflection  angel  at  nozzle  exit  measured  with  respect 
to  fraestream  direction 

Source  strength  density 

Euler  angle 

Azimuthal  angle  measured  around  periphery  of  jet  impingement 
region.  Ft-  ajj  90°,  (|)  ■ 0 in  direction  of  horizontal 
component  of  jet  mean  flow 

Computational  polar  angle  for  stagnation  line  computation 
between  a pair  of  wall- jets 

Ruler  angle,  radial  momentum  parameter  used  in 
Reference  2 

Associated  with  lift-jet  1,  2,  3 or  4 
Aft,  Area 

Centerline,  crosn-flow 
Jet  exit  station,  equivalent 
Forward,  fountain 
Impingement  region 
Jet 

Maximum  value,  mean  value 

Normal  to  surface,  normal 

Total  (fluid  properties),  conical  geometry 

Reference,  rectangular,  radial 

Suckdown 

Stores 
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1.  INTRODUCTION 


The  doflicn  VTOL  alrcrnfl  equipped  with  powered  llfu-jeC  systems  re- 
t|uli'OH  a knowlu.clne  of  the  L’.ompHeatcui  flowfield  produced  l)y  Huch  devices 
over  a wide  ranp,e  of  aircraft  operating  eoiui  1 1 lone , Of  parilcniar  importance 
are  thi'  flowfield  1 ntoraellons  resulting  from  the  operation  of  these  typea  of 
alrc.ral't  In  the  hover  mode  out  of  j;ruund  effect  and  Ln  proximity  tej  the  ^;round 
where  atrontt  Interact  ions  f ri>.i|uei»tiy  o(!cur  between  the  lft-  |et  .streams,  the 
.ilrl'ranie  surlaces,  and  the  pp'oeiul,  Tliese  interactions  nsiially  result  In  tlie 
introduction  of  forces  on  the  airframe  which  can  be  positive  (providing  addi- 
tional lift)  (jr  negative.  Thus,  a propulsion  .system  designed  without  taking 
Into  account  the  inrheced  forces  may  not  provide  sufflcl.ent  thrust  for  an  nde- 
tpiatuly  controlled  takeoff  without  u reduction  In  desired  payload.  In  addition 
to  these  induced  net  loads,  situations  are  encountered  where  unfavorable  moments 
are  produced  ou  the  airframe  resulting  in  significant  stability  and  control 
p robioms . 

The  induced  forces  (and  moments)  in  and  out  of  ground  affect  usunriy 
result  from  one  nf  two  rensonabiy  well  defined  viscous  flow  phenomena!  Jot 
ciurainmo.nt  and  the  fornwtl.on  of  Jet  flow  fountain.s.  det  entrainment  causes 
othi’rwise  siatLc  air  to  be  set  Into  motion,  resulting  in  locally  reduced 
static  iiresHurcs  on  nearby  airframe  under-surfaces  thus  Inltoducing  negative 
aerodyiiainle.  loads,  Tl>u  Jet  entrainment  effect  occurs  hotli  In  and  out  of  re- 
gions influenced  by  tlic  presence  of  the  ground  but  is  frecpiently  accentuated 
as  tlio  distance  between  tlie  Jet  source  (nozsile  exit)  and  tile  ground  is  reduced, 
largely  because  of  the  proximity  of  the  additional  entrainment  resulting  from 
the  ground  waU-Jcts, 

Although  the  Jet  entrainment  effect  occurs  both  In  and  out  of  ground 
effect,  the  formation  of  Jet-flow  fountains  rciiuires  an  impingement  sur- 
face and,  therefore,  is  peculiar  to  operation  close  to  tile  gnuind,  The  forma- 
1 ion  of  foimtalns  l.s  also  configuration-dependent  in  that  multiple  Jets  are 
requlre.il  and  tlie  Jet  impingement  points,  relative  streiigths  of  the  Jets,  .ind 
let  Impingement  angles  arc  strongly  influential  factors.  Tlie  upward  convec- 
tion of  ]e.t  flow  in  the  fountains  usnaJly  results  in  a po.sltlve  aerodvnamic 
lift,  cauised  hy  the  positive  pressurization  of  airframe  under-surfaces  con- 
taining and  deflecting  tlie  fountain  flow.  liecaese  of  the  upward  convection 
of  till'  llft-Jet  flow  in  tile.  fmintnlnM,  a degradation  nf  propulsion  svsieni  per- 
foniiaiici.‘  f ret|uent  ly  rcHuits  through  exliaiist-gas  Ing.entlon.  In  this  resiiect, 
rount.'i  I IKS  can  he  detrimental  to  VTOL  aircrfiil  perl  oritiance  In  j’ronml  proximity, 
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Generally,  the  aerodynaralcH  of  VTOL  aircraft  ara  Influenced  by  a number  of 
geometric  parameters  and  by  the  viscous  (turbulent)  nature  of  the  local  flow- 
field.  Significant  VTOL  hover  flowfield  interactlona  and  their  resulting 
effects  on  aerodynamic  performance  in  ground  effect  are  summarized  in  Figure 
1,  The  comprehensive  theoretical  prediction  of  these  flowfields  and  resulting 
aerodynamic  forces  for  arbitrary  aircraft  configurations  is  undoubtedly  many 
years  off,  and  It  may  never  be  economically  feasible  to  rely  solely  on 
Ihooretleal  predictions.  However,  it  is  possible  wltii  present  knowlmlge  t;o 
synthesize  useful  solutions  with  empirical  information  for  well-defined 
flowfields  about  airframe  components  of  significant  currant  interest.  It  is 
to  this  objective  that  the  work  described  herein  was  directed:  the  prediction 

of  lift  Jet  (primary)  and  induced  flowfields,  including  induced  forces  and 
moments  about  arbitrary  aerodynamic  shapes  in  and  out  of  ground  proximity  for 
a well  defined  set  of  flow  conditions.  This  study  wau  accomplished  through  the 
synthesis  of  new  and  existing  analytical  techniquea  and  empirical  correlations 
in  use  at  the  McDonnell  Aircraft  Company  (MCAIR)  and  the  McDonnell  Douglas 
Research  Laboratories  (MDRiJ . 
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FIGURE  1 

VTOL  FLOW  FIELD  INTERACTIONS  IN  GROUND  EFFECT 

The  combined  metlioJology  wus  developed  bv  pt’occodlug  logicfii  ly  I rom  the 
lift-, let  exits  through  the  viscous  flowfield  inc.l.udlng  the  tree-lets,  w, ill- 
jets,  wali-jct  stagnation  lines,  fountain  formation,  and  fountain  InipingemeiU; 
on  the  airframe  In  ground  effect.  In  and  out  ot  g, round  ellecL,  let  entrain- 
ment induced  flowfields  f>re  obtained  from  computations  of  induced  potential 
flows  about  the  airframe  •is  entabiisned  by  emiiirlcnl.  free  and  wall-jeL  L'lUraln- 
iiuaU.  eharacterl  St  Lcs  where  appropr  lat.e. . Although  the.  viscous  flows  do  Intlueiu'e 
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the  induced  potential  flow  through  specification  of  let  surface  entrainment 
veiocities,  the  induced  potential  flow  effect  on  the  complete  flowField  is  not 
iterated  or  altered  to  compute  additional  Interactive  flowficld  effects.  The 
basic  separation  of  suckdown  and  fountain  forces  shown  In  Figure  2 is  Inherent 
In  tlie  formulation  of  tlie  resulting  me.tliodoiogy . 


SEPARATION  OF  FORCES  ACTING  ON  VTOL  AIRCRAFT  HOVERING 

IN  GROUND  EFFECT 


- — - 

2.  FREE-JE'l  AND  WALL-JET  FLOWFliiLDS 
AND  J1:T  entrainment  0HARACTER.1RT1CS 

Till'  Mn.i  I y f.  I r;i  1 iiiid  (imp  l.r  li'/i  I nuulullnf’  of  the  v:1sc:oiik  (t  urlviilr.'nt)  I'lnws 
liu'ludinp,  till'  rroi'-Ji'lH,  ]i't:  Imp  1 iip,i.!incnt  ri'pluiiK  ntid  wull-Ji'l-n  Is  tli'scr  Ibcd 
in  till'  In  1 low  Inc,  Hiibni'i't  limn.  TIiIh  iU'mt  1 pi.  ion  of  the  flnlil  I’lec  linn  Ion  oi‘ 
llu'ne  flow  plienoniena  provIdoM  a iiieaiiM  for  till'  pri’ti  k'.  t 1 on  iif  the  propertien 
of  tlu'.so  flown  Im'ludlnc  flow  c.i'oiiiet  ;•  < , veloeiLy  tl  1 s t r i liii  I i on  , i:v  Miuil  L nill  fhi>: 
ilintrlbutlon  and  iiians  averaged  tomparature,  Total  prcinHUfea , if  reqirlrod, 
are  obtained  from  knowledge  of  the  velocity  dlatr ibiit iimsi  in  the  constant 
preHHure  (ambient)  free- let;  and  wall.-;]et  reginnH.  The  resulting  prediction 
methodology  Incorporates  u number  of  theoretical,  empirical  and  somi-emiiir leal 
descrip  t Ions  of  the  Lndividua.l  flow  regions  involvedi  The.  emphasis  technic- 
ally, liowever,  is  on  the  use  of  high  quality  empirical  data  whlidi  have  been 
tni'or porated  Into  the  prediction  methodology  wherever  posalble  in  preference 
to  purely  theoretical  pred fc tions. 

Eree-Je.t  Flows 

It  Is  expected  that  for  VTOL  aircraft  In  hover,  the  llft-Jets  will  retain 
free-Jet  characteristics,  i.e.,  they  will  be  to  first  order  unaffected  by  pres- 
sure gradients,  for  some  fraction  of  the  distance  between  the  nozzle  exit  and 
the  ground  plane.  Experience  based  on  wind  tunnel  tests  .ind  proctfcal.  considera- 
tions have  liullcfited  tlint  tlu;  tllstance.s  between  nozzle  exit  and  ground  plane 

(11/15.  to  be  ennsidered,  lie  in  the.  range 

-I 

O.!  < 11/1),  ll 

,1  e - 

Consi'ciuent  ty , interest  in  the  properties  of  free  tnrbulunl  Jets  Is  focuse.il  on 
the  same  range  of  nozzle  diameters  downstream  of  the  exit. 

Tile  velocity  discunt  InuLty  at  tlie  exit  of  a )et  no'zzJ.e  Is  smoothed  out 
by  vlseouH  effects  witliin  a thfn  shear  layer.  The  sliear  layer  is  highly  un- 
Htahle  and  develops  from  the  nozzle*  lip  as  a free  turbulent  ffow,  I'Tir  some 
dl, stance  downatream  of  the  nozzle  exit  plane,  the  jet  retains  a potential  eore 
flow  i.urroundcd  by  a slieath  of  highly  furl)uli>iU  fluid,  where  mixing  hefween  the 
let  and  .surrounding  fluid  takes  place.  This  core  flow  region  Is  eventually 
"t'oiisumud"  by  the  spreading  turbulent  mixing,  .'tnd  the  Jet  eventually  bucumus 
fully  turbulent.  The  mass  flow  in  the  turhulont  region  Increases  along  tlu! 
let,  and  flow  continuity  requires  that  this  excess  mass  flux  be  supplied  from 
' !u;  surrounding  medium.  The  rate  at  which  tlu*  Jet;  entrains  extern, 'il  ai  r 
Is  found  to  he  different  In  the  core  flow  (or  iiofeiUlal  core)  region  from 
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tluit  in  tl\o  fuD.y  turbuliniL  roRlon.  In  thu  former,  entruinmanL  in  duo  to  u 
Hhoar  layoi'  whicli  HeparatoB  two  rpRlona  of  aiiproximattU y conHlank  vuioolty, 
tn  l:tu'  Inttc'r,  an  the  turlnilont  ri'ulnn  BpriwiilH,  tlu;  Jot  i^oiuorlint'  volurity 
(.k’cayH,  Tho  d J ,h l.anc.t'  from  tho  a.t  whlt'lt  tlio  flow  hocomc‘H  fully  Lur- 

liulunt  fa  tluTufoi-u  ai'  Important  iianiMiot or  In  tlio  ca  1 rn  1 a 1 1 on  of  alrfraiin.’ 
Huukdown  foruoH  tluo  to  fri‘0“  )ot  i nt ra  1 nnion l , 

iiu'  axial  il(‘Vo  1 opiiuMil  of  froo  Itirhiilool  lull'.,  and  uonsoiiuuni  ly  lliult  un- 
tralniiK'iit  dlatrihution,  Kenornlly  dopondn  on  tlio  following  paraiiiotor.s : 
o Jet  oxit  Ruynold.s  mimhur  - Hnj 
o Noxxlo  pruHHuro  ratio  - NPR  “ 
o ■ ■ ■ 

o 
o 
0 


Nozzle  tomporature  ratio  - T ,/T 

0.1 


Nozzle  Hliape 
Kxit  voloelty  prof lie 
Exit  tiirlnil  ence  Intena  lty 

2.1.1  ReynoldH  Numhor  Effeota  - For  the  relevant  Reynolds  numher  range, 
tiirhulont  transport  dominates  the  flow.  The  shear  layer  and  jet  develop  aa 
free  turbulent  flows  unuonstrulned  by  rigid  boundaries,  and  e.onsequent  ly 
Reynold.s  number  effects  are  negl  i.'tlbie  (Reference  1), 

2.1.2  Nozzle  PreSHuru  Ratio  - The  nozzle  pressure  ratios  considered 
range  from  stibsonlc  (Nl’R  " 1,05)  tomildly  nndorexpanded  (NPR  ■ 2,2),  These 
are  representative  of  enrrent  ll.ft-fan  or  lift-jet  propulsion  systems.  At 

the  crli  leal  pressure,  ratio  1,H9,  the  nozzle  chokes,  and  sinc.c  only  convergent 
nozzles  are  considered,  larger  values  of  NPR  will  produce  undurexpaiuled  jets, 
I'or  subsonic  Jets,  the  effects  of  NPR  (and  temperature  ratio)  will  he  reflec- 
ted in  a Jet  density  which  Is  different  from  ambient.  The  effect  of  Jet 
to  nmlilent  density  ratio  upon  entrainment  cbnrncteristlcs  is  tnken  into 
account  by  currently  available  entrainment  data.  For  imderuxpanded  Jets,  the 
exit  pressure  is  relieved  by  a supersonic  expansion,  Overexpansion  results, 
and  the  flow  is  rucampressod  by  a series  of  slioc.k.  waves,  Referencu  (2)  Iden- 
tifies tl\e  following  flow  patterns  and  asuociated  pressure  ratios! 

a,  1,05  ^ NPR  < 1,89  - no  sliocks,  exit  pressure  .approximately  ecpial  to 

amh lent 

b,  1,89  ■■■  Nl’R  s 2,08  - weak  expansion  waves  nnd  normal  shock;4  geni’r.il.ed 

in  Jet 

c,  2.08  < NPR  < .1,79  - "diamond"  shook  pattern 

d,  3,79  ^ NPR  - strong  shocks  develop  iit  centiM'  l ine. 

The  nozzle  pressure  ratios  considered  In  this  study  Includi'  the  first  ihi’ci' 
poHs  ill  1 1 1 L ies.  Cases  (b)  and  (c)  are  eompl  Icated  by  the  fact  tlial  the 
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PIQURE  3 

STRUCTURE  OF  THE  AXISYMMETRIC  JET 


l■ro..-|ot■  dovolopmunt  can  bo  mioquatoly  di-Hr.rlbud  by  Mui  b.-havim-  ,,l 
r!hiractin-Lntlc  Jet  vclncitioH  and  IcnRthN  .•iiul  tlui  hu  1 r-M.tmi  lar  pn.in  rtloH  o 
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I'ully  devtilopocl  ,)otH  uC  NPU  ranRa  (o)  of  Uaforancc  2,  Figure  h preHents  t.oii- 
dlmenHlniml  ve.loc'.:lty  pi'orMuH  for  iixiHymmetr Ir.  free-JelB  as  a function  of  '//Dj 
For  Z/l),  > 7,  the  velocity  prufllea  nre  ho  1. f-Hlm  1 lar , /md  the  Bhiidc.  no  lofigo-*' 

varloH  with  X.  Figure  1 preHi'iilM  the  l)ehavl.or  of  tlu'  Jet  char.-u'.tirl'.Ule 
leiigihii  and  veliU'iLleH  w 1 l.h  axial  distauee  ’I'lie  vi’loelly  |M'nl  lli'U  and 

|ihvBleal  geometry  of  l hi"  free-li'l  eau  hi'  olttalned  from  these  iwo  I Ig.iiree.  up 
to  the  beg,  I Uli  I Ug,  of  the  |et  I mp  I ug.emen  t reg.iou  111  the  I'lowtleld  between  I be 
airframe  under'"Hur ftice  and  the  ground  plane. 


PI0URE4 

NON-DIMENSIONAL  FREE-JET  VELOCITY  PROFILES 
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FIGURE  B 

AXIAL  VELOCITY  DECAY  AND  SPREADING  CHARACTERISTICS  FOR  FREE-JET 
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butloi'iH,  wlicn* 

• • « 

”1.11  “■  1 1 

add  )!  jo 

and  ill  ,,  doiiotoH  tiu!  ontra  Lin'd  iiia.Hd  f low, 
iid  tl 

Moat  aiitlioi'H  c.uii.a  Idar  tlio  data  ut'  KLcmi  and  SpaldliiR  (Curve  C)  to  ropruHent 

flu.'  moat  accurate  iiieaHiireiiieiU  ol'  imiHS  entra  Iniiu'iit  la  the  I'nlly  develuiu'd,  re.Rlmi, 

IIIII.'h  data  (Curve  A)  aprouH  well  wltli  Curve.  C In  tin.'  roRloii  ol  Tally  dovc’lo|H'd 

Tlow,  .'iiid  a I HO  LinMiide.'i  me.iHiiremeiil  a In  the  core  Tin''.'  tei'.  Ion,  The  ehaiiRe  In 

entrainment  rate  between  thii  fore  flow  and  the  Tnliy  developed  flow  rey,lun  In 

ovldeiil  In  the  data.  Carve  A hIiowh  that  the  Initial  entrainment  rate  Ih  not  eon- 

at.'int.  On  the  other  hand,  Klela  and  Fohh'  d.ata  (Carve  II)  ahow  a eunHtant  Initial 

entrainment  rate,  with  a hrnali  In  Hlopc  at  (Z/D,  ) Carve  li  la  In  pood 

O' 

aj'.reemeni  with  the  KLeoa  and  SpaldliiR  d.'ita  In  the  lally  develo|'a'd  flow  replon, 
Cnrvi.'  I'  ,‘ihown  a .Hiiialler  a lope  eh.’inp,i‘  hetwaen  the  eote  and  the  lully  developed 
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FIGURE  6 

ENTRAINED  MASS  FLOW  DISTRIBUTION  FOR  CIRCULAR  JETS 


Plow  fugions.  In  thu  fully  dovylopucl  flow  region,  r.hc  slope  for  Cnrvo  K 1h  mih- 
Htnntliiliy  lower  than  tluic  meaHurnd  l>y  othur  Invioitlga tcn-s . Most  of  the.  eiirvos 
In  Figure  6 e.liange  slope  In  the  range 


(Z/l),. .) 


6. 


In  the  anaiyslH  by  Wygnannki  (Kel'erunce  6),  (Z/l).  ) " b Is  taken  to  he  the 

.1  ^ 

beginning  of  the  fully  developed  flow  region.  The  trannltlon  zcjno  1h  nog  lee  f.ed , 
and  a eorc'  flow  entrainment  rate  In  UHod  up  to  “ h,  with  a fully  dev- 

eloped entrainment  rate  for  (F./D|^_)  b.  WygnanHkl's  entrained  raasH  dlnlrl- 
hutlon  1h  indlcateil  by  t’.urve  I)  in  Figure  b,  111m  Initial  entrainment  rale 
1h  reaHonably  ai’eurate;  the  termination  of  the  Initial  region  at  “ Pi 

|u  In  good  agri'ement  with  Curve  B,  However,  the  entrainment  rate  In  the  fully 
di.'vel  opi.'d  region  In  nxe.eHn  Ive, 

In  the  preHent  Hlndy,  the  data  of  Klein  and  Fohh  (eurve  B)  liave  been  hi'- 
lei'ted  to  cHtabllHh  the  imtrainment  rate  for  H\ibHonle  elreular  Jets,  with  the 
lol  lowing  reM\il  ts: 
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iintl 


dO’Tn;;) 

d 

• • 

d(7/l),^,) 

1 Nozzle  'fi'miie.r 

Mtiiri'  U.illo 
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Fur  0 £ < 


« 0.:»2 
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ratio  are  taken  into  acconnt  in  dutermininn  .]ot  exit  veloe.ity,  V , and  Jet 

J ^ 

exit  dentiity,  p . which  In  turn  affects  the  free-Jet  entrainment  rate  aa 

,1  ^ 

shovai  In  FiKurii  ft.  Buoyancy  uffoctH  In  the  Free-, jets  and  wall-Jeta  are  not 

LonHl.dered  in  the  reHuitinfi  methodolony . 

2,1.4  Noity.le  Shape  - Reference  5 contains  data  on  tlie  character istic-s  of 

subsonic  Jets  exlmustins  from  rectanRuiar  nozzloB.  FiRure  7 defines  the 

no.',;',  le  Rcomotry  and  Rtmer/il  ntructuro  of  the  flow.  Tnree  different  rep, ions 

are  Identified  accordinR  to  the  rate  of  decay  of  the  maxinuim  jet  velocity. 

(See  insert  in  FlRuro  7).  In  the  core  flow  ruRion,  the  maximum  veioc,i,ty  Is 

(•niual  to  tl\e  jet  exit  velocity.  In  the  characteristic  decay  rcRlon,  the 

flow  is  fully  turbulent,  hut  the  decay  rate  of  the  maximum  velocity  with  axial 

distance  depends  an  the  noxzle  us|)Bct  ratio  AR,  In  the  axisymmetr ic,  decay 

.1  ^ 

roRlon,  the  maximum  velocity  decays  as  'v  (i/7,),  as  in  an  axisymmetr  ic.  Jet,  The 
velocity  profile,  Irowuver,  Is  not  axlsymmetric.  Finally,  the,  profile  does 
hccoiiu'  as Isymmetrlc  in  what  has  been  labeled  the  fully  axlsymmetric  flow 
reulnn.  All  memory  of  the  no^r.le  shape  has  then  been  lost,  and  tin*  Jet  senn;- 
ture  is  anaioRouM  to  that  of  a fully  developed  axisymmetrlc  Jet. 

In  lift  Jut  VTOL  aircraft,  .ipproxlmaLely  rentaiiRiiiar  Jets  arc  Rcnc'rated 
prliiKirlly  hv  uronps  of  lift-jet  enRlnes  Ruomctrically  clustered  toRether,  Thus, 
a maxlnmiTi  of  AR|  « 4 is  considered.  For  convenience  In  relatinR  aircraft 
and  let  Rcometrles,  the  noxr.Jo  aspect  r.atlo  fs  defined  to  he,  Ktestor  tlian  one 
wlien  the  lonp,  nozzle  dimension  is  perpendicular  to  the  aircraft  plane  of  sym- 
metry, and  less  than  one,  when  it  Is  parallel  to  the  plane  of  symmetry.  Thus, 
tlu;  minimum  aapue.t  ratio  to  be  eonsldered  Is  .ipproxlmal  el  y AK|  = 0,2). 

Kuferomu  5 defines  the  end  of  the  core  fluw  ruRion  and  the  hcRln- 


nliiR  of  the  axlsymmetric  decay  reRlnn  as  a function  of  no;',;',le  aspi'ct  rni  lo. 
D.itii  lor  the  end  of  the  core  flow  rcRion  have,  lu'en  rcplotted  from  Reler- 
eiice  ) In  FiRure  H,  In  a sliRhtly  different  form,  FiRure.  H also  Ineludes  data 
from  llel'erencL'  7,  which  were  presinited  In  Kelereiici.'  ),  loir  ,i  ni|u,iri'  noz;'!,' 
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!,  iZpc'C) 


(AR  " .1),  t'.lu!  core  Flow  rfl^loa  hna  boen  ahown  to  terminate  at  7.p^/c;  •=  6.2. 
For  a circular  jet,  the  indicated  value  of  6,0  is  aomewhat  higher  than  the 
value  of  6 reported  by  other  autliora.  Although  there  is  considerable  data 
.scatter,  an  lnc.rea.se  in  e.ort:  Flow  lengtli  with  no?.Kle  aspect  ratio  l.s  in- 
dicated, peaking  at  ARj^^  = l.A  and  HiibHe.c|uent  ly  ilecre.as  Ing  towards  6,  the 
vahu!  used  In  ReFerence  6 for  two-d  InuMi.slonal  Jets,  IMgure  9 repeats 
the  average  empirical  curve  For  (k|,|,/c.)  shown  in  Figure  K,  and  Indicates  tiui 
depcndiMice  oF  the  onset  oF  ax  I syniniet  r I c dec.ay  on  noxxle  aspect  ratio. 
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FIGURES 

END  OF  INVISCID  CORE  AND  START  OF  AXISYMMETRIC  DECAY  REGION 

Thu  boundaries  oF  the  various  Jet  regions  depicted  above  tire  used  to 
estahilsh  regions  For  Llie  Free-Jet  entrainment  r it-’S.  Ueferetu'e  .6  Indieates 
that  tlu'  entrainment  rate  lor  an  - 10  Jett  initia.lly  is  higher  than  that 

For  a circular  Jet,  but  ;it  the  onset  oF  nxisymmatric  decay,  the  entralr,  lent 
rates  are  the  same. 

For  rectangular  exit  area  Jets,  two  entrainment  rates  are.  usea.  One  is 
va  1 Ul  In  tiu'  core  Mow  ritgion,  .and  one  Is  valid  in  ilie  lully  developed 

turbulent  region.  In  tlie  core  Flow  region,  it  is  assumed  th.it 

. Ji-  'JL.  (4: 

R 2 II  r , 


wherp  - entrainment  rate  for  roctnnguUr  noazle 
Pj^  - perimeter  of  reetangular  nonKlu 
E - entrainment  rate  for  circular  Jpt 

For  the  fully  developed  turbulent  zone,  tlie  results  of  ReFurence  i will  lie 
uHcd,  Reference  5 shows  the  mass  flow  tllstrllnition  for  a circular  ,)et,  and 
fur  an  aspect  ratio  10  rectaiiKular  jet.  These  data  are  plotted  In 
Reference  'i  In  inches  ilownstre.am  of  the  nozzle  e.xit.  Figure  10  shows  the 
same  data  plotted  in  terms  of  axial  distance  normalized  by  a diameter 


d 


C 


(5) 


Satisfactory  correlation  is  obtained,  indicatlriR  that  for  the  fully  developed 
turbulent  reRion,  an  entrainment  equal  to  that  for  a circular  jet  whose  dia- 
meter la  Riven  by  Equation  (5)  can  be  used. 
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FIGURE  10 

MASS  FLOW  DISTRIBUTION  FOR  CIRCULAR  AND  RECTANGULAR  JETS 


In  order  to  fully  specify  the  entrainment  distribution  for  a rectangular 
jet,  the  li-ngth  of  Lh.-  core  flow  region  must  he  defined.  This  can  be  aivomp  I I shed 
by  the  average  curve  in  Fif’ure  8 which  is  based  on  the  data  of  Reference  li. 

This  ciirvi:  shows  the  length  of  the  core  flow  region  (normalized  by  I he  nozzle 
chord  c)  as  a function  of  nozzle  aspect  ratio,  for  the  range  of  interest 
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2.1.5  Exit  Velocity  Profllu  - Circular  jets  with  non-uniform  exit  con- 
ditions arc  analyzed  by  defining  equlv/aljjnt  uniform  jets.  Tho  analogy  has 
lit’on  described  lit  References  H and  9.  I.et 


/’ 


ni|,|>  = 2ii  I P|V,  r dr 


and 


(6) 


• J 

11.1-2^^  r dr  (7) 

o 

• « 

where  and  Mj  denote  the  actual  jet  mass  and  momentum  fluxes,  respectively. 
For  the  non-uniform  jot  a reference  velocity  is  chosen,  denoted  by  Vj^  and 
shown  in  Figure  11.  This  should  be  the  iwiximum  in  the  profile,  which  in  most 
cases  will  also  be  tlie  centerline  velocity.  An  equivalent  jet  of  radius  r^, 
located  at  n distance  Az^  from  the  origin  on  Pi.gure  11,  Is  defined.  The 
equivalent  radius  and  location  are  chosen  such  that,  for  a uniform  velocity 
dlstrihution  the  equivalent  and  ncfual  jets  will  have  the  same  mass  and 
momentum  fluxes. 


DEFINITION  OF  EQUIVALENT  UNIFORM  JET 

Thus : 

<n 

“ 2"  J r dr 

o * 
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and 


J 


H,  2v  I r dr. 

'l' 

() 


(9) 


Sliu'i'  mL'nu'iituni  In  iniiHi'rviid  In  n l*'t,  l'.(|uati<>n  (9)  may  li<^  wrlitiai  ;ik 
r, 


2 IT 


w 2 , 2 ,,  2 

V r dr  - II  r^  p . 


A momiiticum  distortion  coefficient  la  defined 


M 


n 


1 


2 ,,2  2 
which,  wlien  subatituted  in  the  above  equation,  yielda 


(10) 


jV^ 


r = r 

t! 


A maHH  distortion  coeffieiiint  la  defined  aa: 
m , 


n., 


XL 


■itp,  r 


and  f.qimtlon  (8)  may  then  ho  written  na: 


r ^ V j 
j R J o 


(V)a^  rdr. 

V. 


01) 


(12) 


(I’D 


In  Rcferencea  8 and  9,  the  velocity  profile  in  tlie  core  flow  re^^ion  of  a 
uniform  jet  is  deHcrIhed  by  mi  equation  introduced  by  Wnrrcn  in  Reference  10, 


V 

■—  " exp 


- I'll  2 


2 2 

r - r 

c. 


r - r 
v/2  c 


(14) 


wl'icrc 


r “■  inner  radlua  ol'  the  ahear  layer 
c. 


(for  1-  ' r j V/V|^  ™ 1) 


I 


r , ” radluM  at  which  V/V,,  = 

v/2  K 2 

Usinp,  liquat  loiih  (111  and  (14)  in  (I'l),  the  ec|uivalenl  jet  velocllv  prorile 
may  he  Integrated  to  yield 
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, ^ + V2 

c'.  ^ In  2 


f!  ^ A/, 


Tills  const:  I Lii  ton  nn  ImiilU'.ll.  Ion  for  A?.^  , iirovldcd  tliaL  the  clepiMi- 

ilciu'o  of  r , ,,  and  r on  c.  is  known.  I.n  Ki-liironcu  10,  Warri'n  has  obtained  the; 
re  1 a L Ions  li  1 p Ih'Iwiumi  •‘y/2*  from  a,  soiutlim  of  lIil:  let  nioimnnuiii  i!C|ua- 

tlon  and  an  cmii  1 f I I'a  I valiU'  of  tlio  luldy  vlscamlfy  iiuO' f 1 1' i cn|- . The  solution 
loads  to  a comp  licated,  non-Jinoar  re Lat  lonsli ip  between  *^y2’  ^e’ 

In  llefciunce  8,  Wltze  baa  ropreaented  this  relationship  by  a polynomial,  which 
allows  direct  calculation  of  Ai'yy  Wltze's  crpiatlon  1st 


- Q)  * (!r) 


20.62 


in  summary,  a circular  Jet  with  non-uuiCorm  exit  velocity  has  been  re- 
placed by  an  criulvaient  uniform  velocity  Jet  whose  exit  radlns  Is  related  to 
that  of  the  real  Jet  by  liquation  (11),  and  whose  exit  Is  displaced  from  the 
rc'al  Jet  exit  by  an  amount  given  liy  liquation  (.16).  The  Core  flow  lenuLli 
and  entrainment  distributions  discussed  in  Section  2.1.2  may  then  be  used  for 
this  equivalent  uniform  Jet.  If  the  exit  velocity  d istrilniLion  is  luuawn  in 
dctn.n  , the  distortion  coefficients  and  n,,  may  be  I'.aJ cu.lated  directly.  If 
tlic  oxlt;  velocity  distribution  is  unknown,  Heferunce  8 coiUal.ns  calculations 
of  Mj  and  n,,  fiir  the  thin  boundary  layor  c.use  and  for  fully  dovel(3))cd  pipe 
flow.  These  value.a  may  be  used  to  estimate  the  effects  of  Jjel  i^xit  profile, 
d 1.1  tort  Lon, 

2,1,6  llxlt  Turbulence  Intensity  - The  experiments  ol  Keference  A In- 


2, 1,6  llxH  Turbulence  Intensity  - The  experiments  ol  Keference  A In- 
cluded the  variation  of  turbulence  intensity  at  tlie  exit  of  a sub, sonic,  cir- 
cular no/7,U‘.  The  rms  value  of  the  axial  velocity  fluctuation  at  the  exit 
w.as  Varied  from  1,2  to  1.6?',  of  the  exit  veloi'lty,  witli  no  measurable  effect 
on  let  mass  flow  distribution  or  on  the  termination  of  the  core  flow  rep,l,on. 
The  d,i.ta  does  show  a .small  Influi-nce  of  turbulence  lul.cnsU.y  on  t hi'  rate  of 
axial  velocity  dec.ay, 

■’  • •’  The  .let  I mpiniicment  Keg  ion 

The  aecnrac.y  of  description  and  modeling  of  the;  complex  fluid  flow  inter- 
.ictlons  within  the  Inner  ri;glon  has  been  greatly  enhanced  by  the  exlen.slve 
study  of  free  Jet  impingement  by  Donaldson  and  Snodoker,  Keference  2,  The 
expi.u"  Iment  al  measiircinents  reported  tn  Reference  2 allow  ai'  empirical  dcllnl- 
I l.o.u  of  the  Velocity,  momentum,  and  mass  entrainment  d I st,  r Ihn  1 1 ons  for  the 
imliiilcnt  w.'ill-Jets  Including  the  very  s I p,n  I f h-.aiU  eirec!'  of  oblique  |L't 


Intplnsemenu.  Miuiy  oC  tile  a priori  ii h sump t ion h in  the  )ot  impinKement  rugion 
uaccl  in  Relierenee  .1.1.  mny  he  re.laxetl  by  nrnkinK  direet  uae  of  the  Donaldnnn 
and  Hnedeker  data.  As  will,  be  nbown  in  the  par.-iRraphs  fo.' lowing,  these  dnt.a 
will  he  used  dJre.e.Lly  to  cstablJnh  mnnentnm  illstrLhnt  leas  and  stannatJon  line 
and  fountain  loeatlnn  In  the  grouiul  surfi.e.e  wall-jet  flowf.teld,  Tlie  pictorial 
p,ootiu'try  and  assoc  lati'd  noiuenc  latnre  for  an  uhl  Iqiii'  )el  lmi)ln|',einenl.  are  shown 
in  i'lpnre  II’.  TIsa  view  shown  in  i'lgiiri-  111  Js  In  the  plain'  (i|i  = 0;  of  the  free- 
jet  me;ifi  velocity  and  nuinientam  vectors.  i|'  In  the  a/limithal  ani'.le  about  i la.' 

Jut  Implngomunt  stagnation  point  in  the  grounci  pliine, 


FIGURE  12 

DEFINITION  OF  NOMENCLATURE  FOR  FREE-JET  AND  JET  IMPINGEMENT  REGION 


2.2,1  det  ImpinBement  Angle  - The  Jet  Impingement  nngle,  is  defined 

as  the  angle  between  the  ground  nnrface  nnd  the  Je.t-eenter  1 I ne  extended  dowii- 
wnsh  to  the  ground  plane  snrfae.!,  Kor  normal  (perpend  1 eu  1 a r ) impingement  on 
tlie  ground  plane,  ri||  * 90“,  In  the  geiiernl  hover  ilowileld,  U||  Is  a rmictlon 
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of  ir.any  variables  including  aircraft  roll  and  pitch  angle,  nozale  vector  and 
splay  angle,  and  the  deflection  of  the  lift  jet  centerline  by  cross-winds. 

Jet  path  deflections  due  to  cross-winds  are  described  by  a relationship  pre- 
sented in  Kefercnc.e  12: 

2.6 

^ * + t,  tan  0^^  (.17) 

wliere  L and  C are  distances  (see  Figure  11-1)  which  liave  been  normullKed  by  tlie 
Jet  diameter,  and  0^  is  the  complement  of  the  angle  between  the  Jet  velocity  vec- 
tor and  the  freestreum  velocity  vector.  This  equation  has  been  used  with  good 
results  In  studies  of  Jet-induced  effects  In  transition  flight,  (Kcference  13). 

2.2.2  Ground  Plane  Momentian  Distribution  - The  character  of  the  ground 
wall-JcL  flowfleld  is  governed  largely  by  the  location  of  the  impinged  Jet 
stagnation  iioLnts  and  tl:e  azimuthal  distribution  of  the  Impinged  Jet  radial 
momentum.  For  normal  impingement  of  circular  or  nearly  circular  Jets,  the 
resulting  radial  momentum  distribution  In  the  wnll-Jct  la  uniformly  distri- 
buted about  the  Impingement  stagnation  point.  However,  for  non-vertical  Jet 
impingement,  the  nctual  azimuthal  distribution  of  momentum  is  non-uniform 
and  heavily  biased  toward  the  direction  of  the  horizontal  component  of  mean 
Jet  flow,  Very  detailed  measurements  of  momentum  distributions  about 
Inclined  impinging  circular  Jets  wore  accomplislied  by  IDonaldson  and  Snedeker, 

and  their  empirical  data  are  shown  for  two  values  of  H/l),  and  two  values  of 

Je 

Jet  tmpinge.mont  angle  in  Figure  13.  Sliown  also  In  Figure  13*  are  two  solid 
lines  representing  analytical  curve  fits  of  these  data  which  describe  the 
momentum  distribution  about  the  Jets  ns  a function  of  azimutluil  angle  i|j.  Tlu". 
equations  for  f(J))  arei 


iv,:i 


0”  < <j,  < ^5" 
f(^)  - 1.0  + (0.535 


6<|.  + 0.8fi5) 


(18) 


45"  (f)  < 180" 


f('l)  - l.O  + 


-0.7145 


— 82"^^ 

+0.  3327 


'00"  - n \ 

(10) 


Till!  solid  lines  In  Figure  13  are  shown  fur  h0“  and  75°.  The  fuiii'i  Ion 

r(jj)  Is  selected  to  correspond  to  u value  of  H/1)^^^  of  approximately  7 for 
ilu-  curreiU  study.  In  the  form  above,  f(l>)  contains  « 


jl  as  a paranieUir , and 
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FIGURE  13 
RADIAL  MOMENTUM 

Impinging  Jet 


thuH  bu  UHt’tl  I'nr  nny  .Impltinunumt  anuli'.  In  tlilM  Mtiuly,  was 

limltoil  to  angUiH  greater  than  60°,  The  data  of  Refort.MU’.e  2 hIkjwb  that  at 
Ujj.  ■ 45°,  there  is  cHscmt tally  no  momentum  flow  on  the  "backHlclu"  of  thu 
Impingement  region. 

Tt  Hhould  be  noted  that  the  data  fit  in  l''iBurt!  l i donionacrattiH  a ooh  4i)' 
behavior  near  i|i  ■ 0,  and  that  the  momentum  1 evul  on  the  "frnntalde"  (i(i  ■ ()) 
of  the  implnKemeut  region  1h  roughly  kIx  tlniew  uh  high  an  the  "backHide" 

((}i  ■ For  the  larger  Impingement  anglea,  Thla  dlffera  Hignll:  lean  1 1 y 

From  the  Himple  e.oH  if  behavior  aHHumed  in  Referenee  11. 

KquatloiiH  (I.H)  and  (19)  are  Further  modlFled  to  retiulre  that  the  ar.lmu- 
thally  Integrated  momentum  Flux  be  ecpial  to  the  froe-|ot  Intldmit  momentum 
flux  for  any  Jet  Impingement  angle,  Thla  was  aee.omp  1 I ahed  by  norma  1 lx  l.iu; 
till'  area  under  eaeh  curve  in  Klgure  li  hv  ItH  own  area.  The  modified  g, round 
plane  wall-jet  initial  momentum  d lat r Ibut 1 ouh  are  aa  followH! 


f ' ((f)  «■  Y 1 ('!') 
A 


-Jjoo-rT- 

1 + 0.115 


<tml  f((|i)  la  nivun  by  Equation  (J.8)  or  (1.9).  Too  Liaumatucl  monitinuoii  ('.nrrcu'.- 
tlon  Y ta  Bean  to  be  a funut.lon  of  impingement  onivlo,  only  and  hocomoa 

one  for  Uj|  « 90°. 

I'or  lll't-)eta  of  exit  area  aapeet  ratio  greater  than  2 lr\  eloae  ground 
effect  when  , core  flow  fa  prcaont  at  the  Impingement  nolnt,  (aet;  Elgure 
H)  , a f(i|i)  momou:  hill  flux  tllatr  Ihut  Ion  as  shown  sehemot  I ca  1 1 y In  iMuure 
■ s recommeiuled  as  a fun'.'lli'u  of  ;et  Imp  I ngenu-nt  ang.le.  '' j | • ‘■'quallon 

fur  Af  shown  In  Figure  ui  la  based  on  a simple  eoulrol  volume  momentum 
balance  for  two-dimunsionnl  uniform  impinging  .lets. 
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FIGURE  14 

MOMENTUM  FLUX  DISTRIBUTION  FOR  RECTANGULAR  NOZZLE  JET  (AR  > 2) 
CORE  FLOW  IMPINGEMENT  ON  GROUND  PLANE 


2 2 ^ Stagnation  Point  Dlspl lu'.enient  - K xpu r Imen t.i\  1 u bs (>r va Lions  by 
llona.Ldaou  and  Snudeki'i"  (lUd’erenci'  2)  have  eHtahltshed  the  slill'l  of  the  jut 
stagnation  iiofnt  Prom  the  apparent  Impingement  point  In  the  dlreetlon  » 180" 
In  the  ground  iilani'.  These  results  are  presented  In  Figure  IT  for  elreiilar 
Jets. 

“ • Wall -Jet  Flows 

A desurlptlon  of  the  flow  In  the  wall-jet  region  uu  llu'  Hurfaue  of  the 
ground  plane  has  been  formulated  throutth  the  use  ul  the  I'xtenslve  data  pre- 
senlei.1  In  Referenee  2 aiv.l  other  sourees,  As  would  he  I'xpeeti'd,  I low  proper- 
ties In  the  wall-jets  are  depeudenl  on  the  propei'tles  ol  the  Iree-jet  as  It 


Stagnation  Point  Displacement,  AR;r: 


0 


fN 


0,1’ 


0,4 


-O.Q 


-o,a 


-1.0 


-1.2 


-1,4 


90 


Jet  Intpiiiflomimt  «ji 

FIGURE  1B 

STAGNATION  POINT  DISPLACEMENT 

ImpinQing  Jet 


I'ntiTM  tite  Imp  1 Tlti'refnri',  rn.'tt-jpt  ilc'ViO opmi'nt  I’rom  tlii'  ,)ot. 

I'xU  til  tlu’  pifotiml  pliuu!  imwl  bo  OHtiihl.lHliaU  prior  Lo  ilo  f Livl  t Ion  of  i,lu!  wall- 
jot  riowH,  TItiH  may  bo  nooompl  iHltoil  tliroiip,li  uho  ol  Liu'  roHiiltH  proHotUctl  in 
Soctloit  Jill  Till'  Inltla,!  oomlltloiiH  for  tlio  wa.l  1-Jot  flown  aro  alno  ilopon- 
ik’iit  on  tliu  U'oal  flow  at  l.bi>  o.xlt  of  tlu)  implnKomont  ri'nlon, 

2,1,1  Wa  ’ 1-ilot  Vol no  Ity  and  TIilokncHS  Dlntr  Unit  Ioiih  - Tlio  oxti'nH.l.vo 
I'xpor ItiioiUa  1 iHLiaHnromontH  donor Iboil  In  Uoforonoo  2 li.'ivo  ontabllnhod  tho  onol’nl- 
nonn  of  tlio  non-d  linonn  Iona  I volootly  profllo  nbown  In  I'lpnro  1 1'>  for  ilio  uir- 
bnh'nt  wall-Jotn,  Thin  oharao  t i-r  Int  lo  nliapo  pornlntoil  Ihroonhout  tho  wnll-loi 
ri.’o,lon,  and  wan  fotnul  to  ho  I'nnont.  lal  ly  Inilopondont  of  K aiul  I ovon  for  non- 
viM'tloal  Jol  1 mp  1 npomoiu  , With  roforonoo  to  Klppiro  I'.l,  thin  non-d  Imona  I ona  1 
profllo  In  a,‘tHomod  to  donoribo  tho  Inrhnlont  wall-Jol.  for  U hj.  Study  of 
MDRfi  vlnooun  olllptlo  flowflold  noliitlonn  for  two-rl Initmn  Iona  1 ,)otn  aiul  axlnym- 
motrlo  Jot  Imp tnnomon t data  (Koloronoo  2)  Intn  loil  to  tho  ronult  that  Kj  In 
about  tw.loo  tho  ontor  radliin,  |,of  tho  imp  Iiih  Int;  froo-Jot  . Thoroforo, 

from  I'liptro  4 i 


! 

( 

i 

) 


}■ 


NON-DIMENSIONAL  WALL-JET  VELOCITY  PROFILE 


T.mnirtli  tiui  UHc  of  Rqtmtlon  (22)  mv:l  !>.  utay  bu  ubUuiiiedi  Tlui 

rmUiil  Hlrlft  of  ilio  Implnmid  .)cit  Ntanuation  point,  AUj  may  bn  obtained  from 
Kltpiro  I 'i  aa  a func  tion  of  the  Rround  ImpinRomcnt  aiiRli’,  U jji  Tliia  Ih  a rola- 
tlvoly  amall.  corroctlon, 

TluiH  far,  the  neumetry  of  the  free-Jcil  and  ground  radluH  Kj  at  the  bi'y,ln- 

nlna  ol  tlie  radial  wall-let  have  been  duflnud,  ^ relation  for  N,.  and  II  at 

,'j  max 

H ^ K)  Ih  now  required,  (Sou  KlRUi'e  12,)  Krom  MI1HL  turbulent  flow  field  eal- 
I'ulalloiiH  (Uoferunee  lA)  and  liinlLluR  oaHUH  of  the  data  shown  .In  Ueferenoe 
the  followliiR  relation  was  i-stabllshed  for  vertleal  Ji'l  Implnpement! 


max 


- 0.')'.  V, 


,n 


<nr> 


■1. 


(21) 


I H - l<x 

The  geometry  and  results  of  the  MDKl,  two-dlmouHlona 1 turlnilent  flowflold  eom- 
putatliMi  arc.  shown  In  I^'lnuroH  17  and  IH  (Kefereme  1A),  Kxpe  r Inien  ta  1 ibala 
from  Kefereneo  2 define  the  dependenee  of  f ,|,  „ d Ihe  )et  ImplnRement 
attnle  (h/mblnat Ion  of  Kquatlon  (2'1)  with  a eurvi^  III  of  the  approjirlate 

data  In  Kofereiiec^  2 yields  the  folluwhiR  resiilll 
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max 
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FIGURE  18 

NON-DIMENSIONAL  WALL  JET  VELOCITY  PROFILE  AT  EDGE 


OF  IMPINGEMENT  REGION  (X  - W) 


i 


K(|imtlun  (2A)  folaton  1)  ut  tha  iji  “ 0 lultmitlml  poaltlon  to  tlie  centerline 

nui  X 

vcleclty  of  the  Imp  1 Ilf;  Ing  jet  aiul  the  )et  Imp  liij’cment.  annle,r(j|.  r.xpcr  Imenta  1 
raciiHun'meutH  In  Kererem:e  alw.  cliMMiie  I'lc  In.'liavlor  ,il  IJ  witli  iii  I’or  nniwer- 

IllilX 

Ileal  jel  Imp  i nj-p'iiieii  t : 

II 


'»()■’  - 

■ l.n  - I)./,  ( .--p-  J.'-)  siM 

^ max  ^-0 

Combination  of  Kquntionn  (24)  and  (25)  yielda  aa  a iimctlon  of  ji 


(25) 


and 


JT‘ 


9(1"  - u 90”  - a 

U - V + .170( -,“d“^-)iH  - 0.4( p;^“^^)  ain  ()./2j.  (2(0 


max 


L 


The  beliavloi"  of  U with  K followH  the  ehar/icterlat  le  radial,  deeav! 
max 


U (10-  T.-  II 
max  K max. 


(27) 


Tlie  muaHuremontH  of  Nr  reported  In  Keferenco  2 yield  ri'tpiln.ul  InCormatlon  alioiit 
N,,  , 


^5  / ( N ) 


(2H) 


.'Mul , through  lnl;erpretai  ion  of  nu'aHureim.MU:H  of  Nr^  over  a wide  ranne  i>f  radial 
Htat  Iona i 


(K)  ■ + 0.07  (1<  - Kj,). 


(^9) 


The  edne  of  tin'  wall-jel  alu'ar  layer  may  be  defined  from  Hi|uatlon  (29)  aiul 
I’l  Koro  I (i  wh  I eh  ahowH  that 


N “■  2 Nr 
max,  5 


(10) 


To  I'ompl.eti'  '.he  deaer  Ipt  Ion  of  the  wall-let  veloelly  and  fbleUneMM  dla- 
trlhuilons,  a re  1 a t lonah  ip  for  N,.|  In  recinlred.  Thin  la  obtained  by  equating 
the  toLal  radial  wall-jet  raaaH  flux  at  Kj  to  the  total  iiiaas  In  thi'  fri'o-ji't  at 
the  Imp  I in;eim.!nt  pe'lnt,  l.e.; 


/■"  /: 


r„,ll  «,  .!♦  .IN  - 
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*r4*J.  - s'.tT.  ;• 


»tr«r»tj. 
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Tt\ia  may  bi>  ru-wrlCten 


I'  vH| 
mT  I 


f”  f" 

“5 1 I ",.™,  <TI 

^ O Jo 


d cl,|,  » ;n 

mnx  j1 


l.ot  Ll\o  i nt'.i.')»rii  1 o I'  llu'  ari'a  ol'  I hi*  viiltu-liy  pfol'Ilt.'  shown  tn  IMcnn'n  Ih  lu’ 

i(i‘M  ( oh  :ui  I*  , 

111 


n 

1 


(^)  ■ ) d " 1^. 
iiwix  !) 


TIu'n  HuliHt  tt\it  Inf’  II  from  Kqnatlon  (26),  tlui  following  roHul  t 1h  obta  Inud  1 


max 


■ p Nr*  I*  V,  , X 

III  ml  t 5T  in  .|rl 


I 


2n 


‘K)”  - <»  90*  - (i  . 

.66  + .170  ( -pT"  ^111  - 0.6( Bin  .t./2ld.|.. 


I'or forml np,  tlm  IndtL'.atod  lnto(U‘ntlon  and  Ho.lvlnR  for  N 


5r 


Sr, 


90 “ - a 


. 00 


0 1 • V,  , 1.66  + .;17()(- 

ml  r m ,1  c;  f 


p-*o— 11 2 II  - .H(—  j i-T.- 


I 

Thori; fom,  N,.j,  Ih  rolafud  to  tlio  (iroportloH  of  tho  froo-lot  thronpli  ii'j.j.j  ''''cl 
Vj^,  j and  to  the  Imp InRement  ;mR.U'  <»|j« 

KqnatlonH  (22)  throupli  ('3.1)  doHcrllie  ctu'  noomotry  and  vclooliy  dl,Btrlhu~ 
ulon  In  thn  waU-)ot  ronion.  Althonp.h  the  thlcknoHH  of  tho  wall-, lot  rop.lon 
In  not  a fimotLon  of  if,  tho  voloolty  1b  Inf  I ncmootl  hy  if  (f'.qnatlon  (26))  and 
tliorofo I'o , tho  maHH  flow  rnto  In  llio  wall-, lot  will  vary  with  I for  non-vortloal 
lot  Implnpumont  tiH  did  tho  momontnm  flux.  Tlio  maun  flow  rato  will  alao  In- 
oi'oaso  w.lth  It  aH  tho  wal  l-,)ot  ontralnH  amhlont  fluid.  Thin  will  ho  dlHi  nnnod 
fnrthor  In  Sot  l ion  2. '5,2, 

Vor  1 float:  Ion  of  ground  piano  wari-,)t't  flow  prodlotlon  ronuIlM  havo  boon 
ohtalnod.  Plnuro  19  hIiowh  moanurod  valuoa  of  maximum  wall-lot.  volooltv,  U 

max 

at  tho  ImpInHomcnt  rug  Lon  outer  riulluH,  Rj  , aw  oompuLod  from  laiuatloUB  (/) 
and  (K).  Tho.  agroomonl  Is  hooo  to  ho  very  gooil . Thin  In  Blgu  I I laanl  .alnco 
tho  valuoH  of  II  anil  R,  aro  noodotl  for  oimiputat  I on  of  wall-lot  I'lnralumout 
Vo  I or  It  Ion. 
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FIGURE  19 

VARIATION  OF  MAXIMUM  RADIAL  VELOCITY  WITH  RADIAL  DISTANCE  • 
VERTICAL  JET  IMPINGEMENT 


2,[\,2  Wal.1-.?e.t  Eiitr.'Unmci>L  - Wnll-Jet  eiitralnmont  riitoH  o.nn  bo.  obtrnined 

throusli  uso  of  the  roHul.a  preacntcd  In  Suction  2.!). I for  wn]..1~Jet.  geometry 
and  vul,di',.lty  ui..trlbutlonH.  The  radial  wall-jnt  is  a divo'-j^inj;  area  Ilow, 
and  since  the  cmtirc  flow  oc.cura  in  a constant  otatLc  nrossure  environment, 
Severe  thinning  of  the  jet  would  occur  without  the  thickening  effect  of  masB 
entrainment.  Equation  (29)  exprcBses  the  growth  of  the  wnU-jet  tliicknesB 
with  radius,  R,  with  .initial  value  at  R^ , 

The  total  nuiss  flow  in  the  free-jet  can  be  computed  from  empirical  data 
on  mass  entrainment  including  density  differences  betweiMi  the.  jet  and  the 
ambient  density.  At  the  impingement  point,  a mean  jet  density  can  be  com- 
putied.  Let  tills  mean  density  at  the  Impingement  point,  he  designed  as  • 
This  density  Is  used  to  determine  N,.j.  in  Equation  (11),  I'he  Uottil  mass  flow 
Ln  tlie  radial  wall- jet  in  an  ar.lmuthal  sector  between  tlte  angles  and  is 
the  foi  l owing  1 . 

'I'u 

m f(^„  - R]  = 
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MMI. 


uslnp,  Che  definition  of  I*  onct  Equation  (27)  for  the  variation  of  U with 

til  mfl  A 

R,  thiH  bi  tomus: 

li  ilf|i4  (32) 

ma  K j 


m 


l('l>, 


-l>^)  , 


R1  - N, 


I 


Now  uainn  Iviua'.lona  (2y)  aiid  (20)  l or  N^(K)  and  (i|j  I'esiH'r.tivo  ly 


(33) 


(90' 


1 -•  0.4 


15 


.i„i 


di)i. 


The  untrained  maaa  betw(  un  Huccesaive  racill  R^ 


• 

m 

P . 

• 

m 

P 

(<i,^  - 

< 

-©« 

m 

«—  _ 

m 

may  l)e  expruHaed  aa  follnwa: 


ained 
*’nviih  lent 


(34) 


Through  the  uac  of  Equations  (33)  and  (34),  and  a mict.esalve  aiU  of  computa- 
tions along  the.  rrilus  R,  the  mass  l;.low,  untrained  mass  flow,  and  mean  den- 
sity of  the  wall-jut  may  be  determined. 


The.  last  terra  in  Equation  (34)  represents  the  volume  Clow  rate  oC  en- 
trained flow  per  unit  area  of  the  wa.U-jet  .surfnce.  Taking  this  surface  to 
be  approximately  parallel  to  tlie  ground  plane  (actually  it  is  inclined  at  an 
anglt=!  of  approximately  4°), 


A 


sur  far'!' 


R dR  d(j) 


2 


(<t>B  - 9^) 


(3.5) 


from  Equation  (34) 


'jviJijviitujd 

’’aral'iLuU: 


Vm  A . 

N .surface 


(36) 


Combining  Equations  (33), 
obta.lntid ! 


(35)  and  (36)  tliu  foliiiwing  expression  for  V.,  is 

N 


0.14  R.,  V,  , T- 
L J c:  1 m 

9()'’-u  " 

0.55  + .|7n(-~„-ii-) 

c 

- 

1 - 

.4  (--r-o-LL)  sin  | 

L- 

- 

(37) 
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90”,  Equation  (37) 


For  tlui  Hpeclal 
iHK'.omt's 


caHe  of  vorLloal  Jet  Implngemtint , a 


,11 


o.K.  11,  v,^.,  r „.vi 
" 


Till'  lumu' I' 1 r.i  I valiu'  n I'  l;i  approx  tinal  o ly  Tlu'i'oi'oro  : 


.07  3 V 


U-1  (R,j  + l<^)- 


a.|j,  » 90” 


(38) 


E.J.J  Radla.l  Wi.ill-Jot  Moiiumtnm  - Tlio  azimuthal  distribution  of  rnclln.l 
wall-jet  momentum  in  obtained  at  the  exit  of  the  impinKement  region  from  the 
results  presented  in  Section  2,2..2,  FlRure  20  allows  the  radial  dHpeiulence  of 
the  impinRed  Jet  radial  i„omentum,  nlso  from  Reference  2.  Those  data  cioariy 
show  that  wail  shear  stress  can  be  neRiocted  in  tlie  radial  developmen*.  of  the 
wali-Jets  since  the  inuliol  momentum  is  clearly  eonservtnl  throuRh  a large  range 
of  'fills  greatly  simplifies  the  iiKuleling  of  the  wall.-Jet  development 

region,  since  the  (''■M  momentum  functions  will  establish  ttui  azimuthal  momen- 
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FIGURE  20 

RADIAL  MOMENTUM  R DEPENDANCE 
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turn  dtstflbutlona  for  ull  vulueH  of  radius,  R,  up  to  any  wall-Jet  luturar.tlon 
radius.  It  sliouid  he  noted  that  the  data  of  Figuro  20  Inf.  lude  Impingumont 
angles  aa  low  as  A5“  where  there  Is  ussuntlal  l.y  no  nKimontum  f .1  ax  on  the  "baek- 
aldo"  of  the  imp  ingemt'nt  region  (i|i  “ 180°),  Mathematical  ly  exprc.s.'-ual , the 
data  of  I'lgare  'JO  Indicate  that 


;i'R' 


O')) 


la  c.ontrast,  K«|uatlon  ('J-10)  of  Reference  11  may  be  differentiated  with 
respect  to  R to  yield 


a__ 

3R 


1 Ah  COB  ij) 

h r2 
r 


(^0) 


IndlcatinB  dependence  on  R,  This  is  not  a radial  decay  because  of  ground  wa.1.] 
shear,  but  1m  a result  of  the  geometric  modeling  of  the  jet  Impingement  region 
used  in  Reference  11,  This  is  apparent,  since  in  the  modeling  of  Reference 
11,  wlien  normal  Jet  Impingement  occurs.  Ah  ■ 0,  in  wlilch  case  Equation  (40) 
would  Indlcato  a radial  decay  of  momentum  for  oblique  Impingement  only,  Addi- 
tionally, the  data  of  Reference  2 indicate  that  the  height  of  the  wnl l-Jot 
does  not  vary  significantly  with  azlmutlial  angle  for  oblique  impingement 
(Equation  (28)). 

The  slmpllf Ication  of  Che  wall-Jct  stagnation  line  oomputationH  and  the 
supply  of  momentum  to  the  bases  of  fountains  that  follows  from  the  use  of 
results  biised  on  Equation  (39)  as  opposed  to  Equation  (40)  will  lie  discussed 
in  Section  3.0, 
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3.  VISCOUS  FLOWFlIiLD  WITHIN  THE  INNER  REGION 

An  dlsciinaQd  in  Section  2,  the  aeeuraey  of  deficr.lpHon  (ind  mndelInK 
of  the  vIhcouh  flow  in  the-  toRlun  between  the  niuli'T*Miir faen  of  the  /ilrcn'aft 
and  the  ground  plane  has  been  great  ly  enhaneed  by  the  ext(a'slve  ex))(ir  Imeiita  I 
Htudy  of  )e.t  impingement  by  DunaldHon  and  Siu’didti-r , Ui-ferenee  'J.,  iiietu'  and 
other  data  liave  been  employed  to  define  the  propertlea  of  t.he  lree-]ia  and 
wall-jet  formed  by  the  flow  from  any  Individual  lin--|i‘l  nanxli'.  Analyals 
of  a multi-jet  VTOL  aircraft  lift  system  flowficld  requires  u eomplete  de- 
scription of  the  geometry  of  the  aircraft  including  individual  nor.rJc  vector 
and  splay  angles,  nircreft  height  and  orientation  (pitch,  roll)  v/i th  respect 
to  the  ground  plane,  individual  lift-jet  trajectories,  and  jet  impingement 
points  and  jet  impingement  angles  on  the  ground  plane.  These  geomutrlc.  data 
are  required  for  each  jet  to  enable  the  application  of  tlio  fluid  flow  modeling 
techniques  presented  in  Section  2.0,  After  the  fluid  flow  properties  of  each 
individual  wall-jct  '’•.•j  established,  the  wall-jet  Interaction  stugmitlnn  lines 
arc  determined  from  the  procedures  to  be  rloscribed  la  Sections  3.2  for  /luy 
two-jet  pair  and  in  Section  3,3  Cor  a 3 or  4 jet  VTOL  aircraft  multi-jet  lift 
system.  The  final  result  of  the  stagnation  line  e.ompututd.on  procedure  provides 
infi  rmation  on  stagnation  line  locations  in  the  ground  plane,  initial  fnun- 
taln  flow  directions,  and  momentum  flux  entering  the  bases  of  fountains  lo- 
cated above  the  stagnation  lines. 

Theue  procedures  have  been  developed  for  the  flnwfield  within  the  inner 
region  between  the  aircraft  nnder-sur face  and  the  ground  plane,  fhe  Inner 
region  In  hounded  on  the  sldus  by  geometric,  pianos  normal  to  the  ground  pinni! 
and  pussing  through  liivts  connecting  each  pair  uC  Jet  impingement  points  on 
the  ground.  The  procedures  developed,  hovuver,  ran  be  easily  extended  to 
ground  plane  flow  regions  outside,  of  the  inner  rc'gion  where  fovmtain  Implnr.u- 
ment  on  aircraft  undertiurf nces  c-:in  also  occur  and  produnu  Joads  on  the  all— 
f rumo. 

^ Airframe.  Jet,  and  firound  Plane  Geometry 

The  initial  airframe  goometty,  aircraft  orientation,  and  height  ubovu 
ground  must,  of  course,  be  specified  or  known  by  the  user  to  implement  the 
prcci i'.'. tlon  methodology,  A three  viow  drawing  of  the  aircraft  including  noK- 
xlo  geom.'try  and  nozzle  center  of  rotation  o'  vectoring  lo  ivrc  loi  atlcus  and 
riltr.raft  center  of  gravity  locution  Is  helpful  for  this  purpoS'!,  It  Is  Impoi-t- 
ant  that  a c.onsistent  right  handed  system  of  coordinates  he  used  jr,  the  air- 
frame .such  as  t:lie  conventional  Fusulago.  Station  (x),  Ihitt  hlae  (y),  Water  hliie 
(z)  system. 
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A muthematiciil  procuciure  has  been  developed  to  determine  the  in  torn  uc  Lion 
points  of  the  free-jet  eenterl Ines  wl.tli  the  Rround  piano,  for  an  arbitrarily 
orit'nted  alri'rnft,  InohicllnR  ttio  effeetH  of  Individual,  jut  dolMuutlon.  This 
procedure  lnc..ludeH  the  uffec.tH  of  eroSH-riow  on  the  freu-jet  patli  defiectlon 
and  the  correnpond  Inp  Impingement  point  a.  Tlie  mathemat leal  formulation  of  this 
proce.duri'  la  deHc.rihed  In  Appundlees  I and  I I . A coniputi  r pro^>raiu  has  been 
written  to  cany  out  the  ca  I cn  I .it  Ions , aiul  has  been  clicikcd  out  for  Ncver/il 
cases.  The.  re.qnlred  Keometrle  InpuLs  to  the  prop, ram  ai-ei 

1,  Aircraft  rofareneo  c.g.  I leatlon  - Kuselapa  Station  (K.S.),  Ihitt  Line 

and  Water  L.liv.*  (Wtt..) 

2.  Number  of  jets  (n  maximum  of  four) 

.1.  dnordlnatus  (I’.S.,  11,1.,,  W.L.)  for  the  center  of  thrust  of  each 
no?./ li> 

4,  Vectorinp  auR^u,  splay  unRle,  jet  velocity,  and  jut  diameter  for  nach 
noxzl o 

i.  Aircraft  pitch,  roll,  and  yaw  unKlus;  huipht  of  the  referunce  c.r, 
above  pround;  and  freostream  velocity, 

The  proyrnm  outputs  the  center  of  thrust  coordinates,  thrust  vector  components, 
jut  imp  lUKe.ii'unt  points,  and  jut  ImplnRement  anRio  complements,  u,  rufurencud 
to  a prouiul- fixed  coordinate  system  wliose  orlRin  lien  directly  below  the  alr- 
crall;  c,r,,  and  whose  >axis  is  allpnod  wUli  the  cross-flow  velocity  vector, 

Tlio  jet  imp Inipjnunt  anple  coinpl uments  are  definod  by  tlie  vector  tangent  to 
each  jet  pnth  at  Us  point  of  intersection  with  the  prmiml,  and  tlie  normal 
vector  to  the  ground  plane,  The  absolute  aiiRlu  In  the  );rouiul  plane  (0°  to 
IbO*)  of  the  Intersection  with  the  Remind  plane  of  a plane  normal  to  the 
uround  plane  and  contalniuH  the  t.ingent  to  the  jut  path  at  Its  Intersection 
with  the  Rrmiud  is  ulso  output.  This  anp.le,  c,  defi.nes  the  luiRative  of  the 
direction  In  the  ground  plane  c.f  the  projection  of  the  horizontal  component  of 
jet  velocity  In  the  ground  plane,  Thin  in  turn  defines  the  ■ IHO”  dlrec- 
tton  for  eac.li  jut  implnuemont  rcRion  In  the  Rround  plane  coordinate  system, 

At  the  ImpinRe.nient  points  in  the  ground  plane,  a b,.i.iil.e.  polar  eoordliiate 
system  Is  used  to  delMue  the  properties  of  the  will-jets  ahoct  each  linplnue- 
meiit  point  as  shown  In  KlRure  12.  This  Is  the  local  K,  ji  system  of  coord- 
inates which  apply  to  each  radial  wall-jet, 

J,2  Stauiiul.Lon  l.lne  f’omputatlon  for  a Pair  of  Wall-, Jets 

The  computation  of  tiui  fitannatlon  line  position  lietween  a pair  of  Inter- 
acLliiH  wall-jel.s  Is  accompli  shod  IhrouRh  the  application  of  a coni  ro  1 -vo  1 ume 
momentum  hal.'Uiri'  a I onp,  a differential  senmeut  of  tlu'  st.’ipnat  Ion  line.  Momentum 


32 


t'liixca  cnterint^  ttio  I'lintnil  voUimu  ;iro  rulatucl  to  tlu-.  imp InBOHiont  ruRion  exit 
momentum  diatr.lbutiona  for  onch  Jot  uh  shown  in  figure  13  and  tlio  angular  sec- 
tor emanating  from  each  Jet  Impingement  point  that  Is  mihtondod  by  the  differ- 
ential Htagnal.Um  line  Hegment.  The  result  Is  an  equation  for  the  slope  of 
the  stagnation  line  at  any  point  In  the  ground  wall-jet  nowfli'ldt  i’he  stag,- 
natlmi  line  Is  di‘t  erni  Ined  by  a stepwlsi'  marching  process  using,  these  c.umputed 
slopes,  A similar  development  lor  tlu'  sliqu'  of  llte  stagnation  line  between 
ene  arbitrary  wall-jet  aiul  a unitorm  g, round  ern.ss-- I'l  ew  has  been  ai  comp  I 1 shed , 
3,2,1  Dovolopmunt  of  thu  Stimnution  Line  Hlope  Kguntlons  - Tiie  momentum 
balance  for  the  determination  of  the  local  stagnation  line  slope  Is  applied 
to  the  elemental  control  volume.  In  the.  ground  plane  as  shown  in  Figure  21, 

The  c.ouHorvatlon  of  momentum  normal  to  the  stagfwttion  line  for  the  control 
volume  In  the  ahsunce  of  any  pressure  and  viscous  forces  Is  the  to.i.lowlngi 

IL,  (p  li  * dA)  «•  0 (4.1) 

N £, 


CONTROL  VOLUME  AND  GEOMETRY  FOR  STAGNATION  LINE 
MOMENTUM  BALANCE 


Slde.s  I and  2 of  the  control  volume  are  chosen  normal  to  the  radial 
Velocity  vectors,  tl|^j  and  I com  the  tndlvldual  Jet,  lni|i  l.nytement  points, 

Assuming,  ne  less  In  wall-Je.t  radial  momentum  with  g.eometric  loeal  radius  K 
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for  e^^cll  based  on  the  experimental  results  shovm  in  Figure  20  (or 

■ IB  uxpreBSud  by  fCquatlnn  (39)),  the  momeiit\im  flux  for  eaeh  (’.Ireunifcfuntial 


Huetoc  Is  ! 


ii  u " h (I  N “ r(i|))  M, 


whore  N is  tl>e  cllstaneo  normal  to  the  ground  plane.  Applying  lCquuti<m  (Al) 
to  the  control  volumu  (Figure  21),  and  using  the  fact  that  no  flow  occurs 
through  control  volumu  sides  3 and  4,  tliu  following  rustilt  1')  obtained  ' 


p Rj  sin  (0  - 1..^')  dN 


II  ^ 
P 11,^2 

^ n 


sin  ((|)2'  “ 0)  dN 


where  tli"'  angles  are  further  dullned  in  Figure  .12, 


FIGURE  22 

WALL-JET  INTERACTION  GEOMETRY 


02-0^2 


.-wT  *^'2 


Subatitutlng  Equation  (42)  Into  the  above  result • the  following  la  obtalnedi 


i’i('l'|)  Mji  Htn  (0  - “ *0  . (44) 

Ki'orn  tho  HLninii'try  Hhown  l.n  Klj>urc!  22,  iln.  In  1 low Ui);  is  ubt  n 1 lU'vl : 
lIv.  Hln  (0  - i|i  ') 


tU  aln  ((!>,,'  “ 0) 

I _ -..  ■■■  ■ A. 


SiibHlltutliiR  oquations  (45a  uiul  h)  Into  Equation  (44)  tlioru  roHii.lLHi 


’'aI  "u  ^A2  "U  '*2’  - »> 

-j —a . (Af,) 

1.  2 


Solving  Equation  (46),  the  followinu  ik'uinnl  roHult  in  oblahu'tl  al'tor  Homu 
furthor  Roomotric,  suhatltntionH  are  made  from  tho  Ronmotry  In  Fluuro  Z2i 


, „ dv  « ’^2* 


wh  nro 


S ^ U2/RJ  R ' 

Tlu!  cUh t itu'.t Ion  trncla  In  tlio  abovo  oqnatlonH  botwoon  tin.’  anivloH  i|)  ami  iji'  in 
nnr.t'HHltnted  by  tho  fact  tliat  tlio  i|''  ■ n d iroot Ion  I.h  parallel  to  a .lino 
JoininK  tho  Jot  ImplnRomont  points,  whoroos  i|>  0 U rulatud  to  Lhf  indi- 

vidual Jet  ImplnHoment  HaomctrtoM.  and  i aro  rulaU'd  hy  a almpli.'.  rotation 
about  the  ImpinRomont  poJtit  of  tho  typo  of  data  shown  In  Figure  I'.l, 

Thu  stagnation  line  solution  for  any  Jut  pair  la  obtaiuud  by  a Htupwlau 
mnrt'.hlng  proaess  heginnlng  at  the  point  where  the  Htagnat  Ion  line  rroHHua 
the  line  Joining  the  Jet  impingement  pnlnta.  At  this  eommon  point,  lliu 
radii  Ik  and  as  shown  in  Figure  22  are  re.lati'd  an  iollowat 

1 »i 


IS 


II  In 


TIUm  may  l)t>  oliLaliu'd  1 i‘um  I'ani.U  Inn  (4(>)  liy  Hnlllny,  i|' j ' ■>  0 and 
lliai  i.'i|iial  Inn,  I'nr  normal  Ird  I m|i  1 ni’i'tiuMit , K<inal  Inn  a I mp  1 1 I li'ii  In 


(!,0) 


h (’.nnirintatlona.l  proni‘<>m  Iuih  bot-n  wrUtoii  ti>  dul.urminci  Lbu  lanirdlna tie  of  Uio 
MlaRtiaLlnn  lines  uh  JiftermiiU’cl  Ihi-ounh  tliu  nohitioiiH  of  Kc|untlonH  (A7),  (i!t8) 
and  Thi'  mnaliuniofi  of  tl»i»  proaediiru  and  tlie  Implamontation  ol‘  Llio  pn>- 

cuHH  for  multi-, lot  flowtioIdH  In  dlHcuHHod  furthor  In  Siu’.tlon  3.J. 

3,2,2  Wnll-.lot  1‘air  Staiinaltou  Lino  riol.utloiiH  - Si/i  gnat  ion  1 iiu'  hoIu- 
tlona  obtalnud  tbrounli  HoJ.utl.on  of  liquatlonn  (A7),  (iltB)  and  (^9)  aro  pri'Honiod 
in  KiHiirc'N  23  tlirounh  25,  Fifiuro  23  hJiowh  tho  roHulta  of  a parainotrln  varia- 
tion of  tho  Intoraot loii  of  walJ-jutH  formed  by  two  vortlca.1  JotH  of  varlouH 
moiTiontuni  ration,  Shown  also  in  Fluuro  23  are  Mtagnation  lino  aolutionH  for 
non-vur tioa I but  oquai  Jot  implnp,emnnt.  nnKioH  of  75  and  60°  for  Join  of 
equal  momentum.  Note  that  the  romiJ tN  ahown  in  Kinure  23  aro  Hhown  for  the 
upper  riHlit  hand  quadrant  onlvi  and  that  tho  ordinate  Hcalo  1h  nompnmHod  by 
a factor  of  two. 

Tile  Htamiai  iim  line  patterna  for  multi-let  Intorac.t  Ions  mav  bn  obtahnul 
from  results  sucli  as  Hhown  in  Figure  23  by  appropriate  location,  rotation, 

/ind  application  of  lot’lcaJ  ruitm  of  symmetry  to  individual  .jct  pair  stanna- 
tlon  lino  Ho.lutiouH,  For  exanplu,  n family  of  wtanuatlon  .line  aoJutioiiH  fi^r 
a four  Jet  arrannemeat  Is  shown  in  Flmire  These  wuro  couHtructed  Prom 

curves  hoI  jcLcd  I rom  Flpure  23.  hiluiwlKc,  Fluore  25  prcMtiits  a w.ill-Jct 
inlurai'tlon  staunatioa  1 ine  solul  Ion  for  a tliree  Jet  arranno.mant  with  the 
two  alt  JotH  each  havlna  one  lialf  tho  Jot  momontuui  of  tlie  forward  Jet  (all 
,1c:  ,4  are  Impiiiuinp,  vertically), 

I’iiiure  2fi  shown  a compariHon  of  computed  .and  measured  suip, nation  line 
locations  for  a two-.lot  pair.  The  oxpi'rimontal  dat  i wore  ohlaliunl  as  i.’  i i. 
of  a McDonnell  Aircraft  (lompany  I'roinilHlon  Department  Indi'in'iuMMii  i icn.iri  !i  and 
divelopmeat  prop, ram  (lUAD), 
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FIGURE  23 

COMPUTED  JET  PAIR  STAGNATION  LINES  IN  UPPER  RIGHT  QUADRANT 


FIGURE  24 

COMPUTED  STAGNATION  LINES  FOR  4-JET  ARRANGEMENT 


STAGNATION  LINES  FOR  3-JET  ARRANGEMENT 
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FIGURE  26 

COMPARISON  OF  MEASURED  AND  PREDICTED  WALL-JET  STAGNATION  LINES 
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3 • lit 3 Local  PltHctlon  of  Fountnln  y].ow  Abova  Stagnation  Line  - An 
estimate  of  local  fo '.ntain  flow  direction  directly  above  each  point  on  the 
Htannation  lino  may  be  obtained  tbrounh  further  uaaumptlon  of  conHarvatlon 
nl:  momontiim  wltlimit  viacoua  Ioshuh  in  tlie  turnini;;  region  at  oai:h  point  on  the 
Htagnatinn  line.  ConHider  the  stagnation  line,  geometry  mul  momentum  flux 
deeompoH  It  Ion  shown  In  Flgiiro  27.  The  loenl  fountain  flow  plane  shown  in 
I'lguro  27  Is  assumed  to  bo  i>erpenclleu  lar  to  tin;  ground  plane.  Here  the  moment- 
um flux  Incident  on  <i  stagnation  line  Hogment  from  wail-Jet  number  2, 
in  resolved  into  components  normal  to  and  tangent  to  the  stagnation  line  seg- 
ment. If  it  is  asHumed  that  the  vertical  momentum  flux  In  the  fountain  from 
« « 

wall-jot  2,  Mp^2i  equal  in  magnitude  to  1^'’®  tangential  com- 
ponent of  momentvim  flux  in  the  fountain  from  wall-jet  2,  is  equal  in 

muRnltude  to  lipji  the  angle  K2  is  defined  as  shown  in  Figure  27  anti  may  ho 
found  to  be  geomotric.nlly  equal  to  t|)2'  - 0. 


Plans  Containing 


riwwnn 

RADIAL  WALL-JET  AND  FOUNTAIN  MOMENTUM  RELATIONSHIP 


ami,  1 Iki’wlni'  I'ui'  flow  from  wall-joi  1, 


I’or  llii'  t wn  sliloM  of  till'  foimi  .1  In : 


M * M 4*  M 2 M 

‘^'VNTOT  "fNI  FN2  VN 


KTTO'l'  I'Tl  l‘T2 


Tlif  annlo  ► j,  if*  iloftiunl  for  Liu*  fountain  ..h 
M 

i'NTOT 

tnn  K.,  “ 

N 

I'TI'O'f 

TlirouKli  fuibut  I tot  ion  of  EquutioiiH  (51)  anil  Into  l',f|uaf  ton  (')'>)  t.lncrt* 

rt'HU.l  tHi 


fan  r,, 


M -f  M 


(llvldliu',  throuKh  by 


tan  n',, 


!Wl4,!!o2 

^Vn  ^'*I'N 


now,  uHlnj;  tlio  anj’loH  > . and  tlu*  abovo  booonu'H: 


tan  > 


I’.ot  K,  + cot  c 


IlaHud  on  tiu!  aHHiimptlon  of  n coiiHtunt  convoctivi*  voloclt*/  In  t.lu*  turnini 
procuHH  botwuim  tlu*  wall-jotH  and  each  hUU*  of  tlu*  fountain,  rcnulLlnv,  in  a 
mixed  and  avnranod  convoc t tvi*  velocity  wltbln  tlu*  fountain  lluclf,  tlu*  aiuO  o 
f In  Kt|iiatton  ('i'j)  Ih  uHNumud  to  roprcHcnt  the  fountain  vclni'iiv  compoiu'iu 
ratio  OH  well  ax  the  momontum  ooinpoiicnt  ratio.  Therefore! 


(59) 


f 

t 


tmi  jj- 


Sntot 


FTTOT 


cot  + cot  ^2 


Tl>u  l.ucii.l  Hlnpo  of  Uhi;  Htiifpii.it  ion  .l.inu  In  the  proiind  pl.iinc,  (I,  and 
thf  iinp,lo  Mu'n  dotormlno  tlui  liu-nt  ci  1 rcc  t ton  of  tlio  iipwaicl  f- 1 i.iw  in  tlio 
To  on  I.  I'l  .1  n i'hovo  o;u'li  point,  on  tlio  stiinniition  lino  jih  kIiowr  In  I'lp.iiri'  27, 

i . 2 , A S tagiiiil  1 on  I.  l_iio  S 1 o_p_o H'liL'L’ Jd’*'  S I lu;  1 o W.a  1 ojal 21*1*1' 

CroBB-l'Mow  Intoriictlon  - The  theorotical  clovo.l  opment  prosented  in  llociJon 
3,2.1  hiis  been  modified  for  the  dotern’ination  ol  the  ioe.a.1  aiope;  of  the  stag- 
niition  line  in  the  rt'.j’Lon  of  Intcractiion  of  ;i  .slnf'ie  Wiill-Jet  and  a uniform 
I'.ross-I  low,  A diagram  presenting  the  geometry  of  this  interaction  in  the 
ground  plane  ;md  the  necessary  nomenclature  is  presented  in  Figure  28, 


Flux  per  Unit 
Length  In  Y 
Drection 


aPTi  emu 

FIGURE  28 

INTERACTION  GEOMETRY  • RADIAL  WALL-JET  AND  UNIFORM  CROSS-FLOW 

A control  volume  enclosing  the  differential  line  segment,  dJi,,  is  used 
to  i'pplv  the  laWH  of  conservation  of  momentum  in  a direction  pe.rpendic.ular 
to  tlu'  Hill!  iiegment  in  the  ground  (x,  y)  piime.  As  la  the  case  of  two  radial 
wall-  jetH,  the  resulting  expression  is  an  eciuatlon  for  the  slope  of  the  line 
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Begment,  cU,  In  the  x,  y,  piano.  The  inometitum  illstribution  in  the  radial 
wall-,-)et  i?i  defined  identically  an  it  was  for  the  two  wull-Jet  interaction 
case.  The  new  physical  quantity  required  is  the  momentum  flux  per  unit 
len^;th  (!n  tiie  y 'ilreet  loti)  tif  the  cross-f  low. 

Aloii)>  the  X axis  in  !''ip,ure  UK,  wliore  the  axis  is  e.liosen  parallel  to  the 
ensH-flow  direet'un,  tlie  pitnetrat  Ion  tif  ttie  radial  wall-.)et  into  the  e.rons- 
I'low  is  I’lut'n  hv  the  iiiatpi  I taide  oi  the  wali-le.t  radius  U , 
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'■'hlch  for  a vertically  impinging  )et  becomes 

M., 
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3Y 


(60) 


(61) 


Through  the  use  of  Equation  (60)  or  (61),  the  location  ol  the  stagnation  line 
us  it  crosses  the  "x"  axis  may  be  found,  Tin:  general  equation  for  the  slope 
of  the  stagnation  line  is  the  following: 


tan  0 . 

1 *f  (I)  C.OH  (f> 


(62) 


whe  ri' 


f(<E) 


2.1111 


3ftc 

3Y“ 


(63) 


These  equations  have  the  same  general  form  as  the  corresponding  equations  for 
a pair  of  intersecting  radial  wall-jets.  Tlie  solution  of  the  e(|uatlon  follows 
the  same  procedure  as  for  the  Intersecting  wall-jets  which  has  been  described 
previously . 

Itiu  practical  cliffiiulty  in  the  appiicatlon  of  the  above  results  Lli 
cross-f low/rnd ial  wall-jet  internctions  is  the  determination  of  the  quantity 
3M  /3Y,  The  difflcuitv  lies  in  determination  of  how  much  momentum  in  the 

i’ 

vertical  direction  (normal  to  the.  ground  plane)  should  be  included  in  the 
momentum  derivative,  Mathematically: 


It  la  rucomniBndGd  thrtt  until  pertinent  empirical  ver  1 filiation  data  are  oh- 

taiiiL'd,  N be  taken  to  be  ecnin!  to  the  value  of  N ciomputed  for  the  radial 

' max 

wall-jet  flow  (Kquatlon  (30)).  Because  of  the  uncertainty  Involved  In  defini- 
tion of  tlic  erofis-flow  momentum  maKnltude,  the  ,«taf»nar.lon  line  .solution  for 
the  wa  1 1 - jet/eross-f  low  itiLeraet  Ion  lias  not  been  Included  in  the  ureund  flow- 
field  computer  proKram. 

' • ' (lompulallon  of  .Stagnation  h files  for  Thriu'-o  t;_^J'[ou  r .let  Flowflelds 

Till'  I'omputat  ion  of  wall-jet  sLannallon  lines  In  (lie  p, round  flowfleld 
below  a multi-jet  VTOl,  aircraft  is  accomplished  from  l.he  lesultK  of  compula- 
tions butween  pairs  of  InteruotliiK  wali-Juts.  The  fuimtnlns  formed  above 
the  staKiuitlon  lines  art=!  usHumed  to  convent  the  interacting  flows  upward  and 
away  from  the  ground  surfacu  therefore  preventing  the  mutuai  Intoraction  of 
any  more  than  two  wall-jets  on  any  stagnation  line.  At  points  in  the  flow- 
field  where  two  or  more  stagnation  linos  intersect,  the  fountain  flows  from 
more,  than  one  jet  pair  will  Interact  forming  what  luis  been  called  a central 
fountain,  llowovcr,  this  intoraction  takes  place  above  the  ground  plane  and 
occurs  only  at  a "point"  In  the  ground  plane  surface.  Thu  fountains  formed 
above  the  stagnation  lines  are  three  diiaenaional  unless  the  ground  stagnation 
line  is  a straight  line.  If  the  ground  stagnation  line  is  straight,  the  foun- 
tain velocity  vectors  lie  in  a common  plana,  but  are  diverging  is  nature  forming 
a "fan".  Fountain  flows  will  be  di.jcussed  in  detail  in  Section  5 and  6 of  tlvls  rc(iort. 

3.3.1  (leomctry  fur  Three  or  l''oitr  .let  Sta&nntlon  Line  Computation  - For 
the  three  or  four  jet  case,  the  results  of  Section  .3.2.1  are  applied  to  each 
wall“jet  pair.  The  ground  plane  geometry  required  is  defined  in  Figure  29. 

Hero  the  ground  plane  reference  coordinate  system  is  designated  as  XOP,  YOl’ 
where  the  XOP  axis  Is  aligned  with  the  direction  of  the  f rue-stream  or  croiis- 
wlnd  velocity  direction.  The  location  of  the  implnguinent  points  XOPN,  VOPN 
where  N is  the  jet  number  are  obtained  directly  from  the  results  dc.'^ic.ribed  in 
Section  3.1  and  Appendices  I and  II.  A straight  line  is  then  constructed  be- 
tween each  pair  of;  impingement  points  forming  a triangular  nr  quadrilateral 
ns  shown  In  Figure  29.  The  angles,  of  each  line  joining  the  p.ilrs  of 
Imijingemunt  p.il.nts  are  then  determined.  The  angles,  iic.,  must  be  determLned 
in  the  absolute  sense  from  0 to  3b0°  with  m ■ 0“  taken  conventionally  along 
the  +XOP  axis  and  +ui  measured  counter-clockwise  from  that  reference  position. 

The  sense  of  each  straight  line  joining  the  jet  Impingement  point  pairs  is 
taken  in  the  order  of  the  numbering  of  the  jet  ini|i  Ingeraent  points,  For 
example,  tile,  line  joining  Impingement  points  1 and  2 in  Figure  29  has  assoc- 
iated with  It  an  ii)|  direction  of  nllghtly  more  than  90'’,  and  likewise  the 
value  of  Is  si  igiitly  less  than  360". 
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FIGURE  29 

COMPUTATION  OF  STAGNATION  LINES  arrroin.ii 

The  J.ni.1  Iv Ltlual,  HtnKntitlon  line  HoJ.utliniH  lU'e  then  computed  for  each  wall- 
jot  pair  in  the  donhl.c-polar  coordinate  syatom  Nliown  In  Kluure  22,  Tlie  rc- 
auJ.tH  of  thuHf  computations,  which  determine  tlie  coordiiv-UcH  of  flu*  stagnation 
lines  in  the  local  polar  sysfoitm,  are  then  ttunsformod  to  tlu>  XOP,  Y(H’  ground 
plane  reference  system  throup*  the  use  of  the  angles  uij^}  the  lorm.l  polar 
coordinates  U,  (|)'i  and  the  Jot  impingement  point  coordinates,  XOPN,  YOPN. 

Kacli  stagnation  lino  Js  computed  In  the  positive  sensu  of  0^ ' , thereby  dotar- 
mining  the  portion  of  the  stagniUlon  line  In  the  Inner  region,  Stagnation 
lines  of  "symmetry"  (perpendicular  bisectors  of  the  lino  Joining  the  Jet 
Imp  Ingi'iiicni  point  pairs)  must  be  dutermlni*d  separately  sini'c,  in  this  case, 
tan  n ■*'  and  l>|uaf  Ion  (47)  Is  singular  since  the  denominator  on  the  right 
lianci  side  of  the  equation  becimes  zero, 

l.'i.H  Rotation  of  Hanlc  MomeiU.nni  II I s tr  I hut  Ion , r(i|i)  - In  the  computa- 
tion i>r  the  wall-jet  pair  stagnation  line  solutli'ns,  the  loial  polar  coord- 
inate system  R,  0'  will  not  always  result  In  nn  orientation  sueh  that  ili'  ■ 0 
In  the  tllrectlon  i **’  0 as  defined  for  each  Ji't  Imp  lagemont  point.  In  these 
cases,  the  angular  rotation  of  the  momentum  distribution  f(0)  as  sliown^  in 
I’lg.nt'e  10,  must  ho  uLlllne.d,  In  l'’igure  '10,  the  angle,  *1,  descrlhod  In  Section 
1,1  I'stahllshes  the  projection  of  the  diroctli'n  of  the  negative  of  Jet  mean 
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PIQURE30 

ROTATION  OF  MOMENTUM  DISTRIBUTION 

VL’lucUy  in  tlin  Hround  plane,  This  corresponds  to  the  projection  of  the 
thrust  vector  In  the  uround  plane  If  there  Is  no  cross-flow,  The  negative 
Y direction  then  corresponds  to  the  direction  in  which  cjj  • 0,  The  ij>  used  to 
determine  the  momenttim  flux  distribution  f(i|i)  from  Figure  13  is  then  given  by 
the  I'ollowing  (from  Figure  30)  i 

^ - 1 + 180“  - M 360“  - (.1  - <|.’  1 (65) 

wlten  I ^ 180“,  ■)'  - 360“  - i|i,  and  the  local  momentum  flux  is  obtained  from 

Kquiitlon  (18)  or  (19)  as  f(<ti)  nr  F(«t')  accordingly. 

3,3.3  Determination  of  Wnll-.Tet  Momentum  Flux  Katerina  the  lUise  of  a 
Founta'tn—  A Mtagnatlon  J ine  pattern  For  a typ  lc.nl  Four  jet  rase  is  shown 

Mchomat tcally  in  Figure  31,  The  total  momentum  flux  to  the  base  of  the  foun- 
tain locat'’d  above  a stagnation  tine  segment  such  as  A-Ii  In  Figure  31,  Is 

obtained  by  Integrating  tlie  momentum  flux  in  a sector  subtended  by  the  angles 

and  ijijj’  also  shown  in  Figure  31.  The  geometry  is  defined  by  the  appro- 
priate valueM  of  coordinates  of  the  points  A and  H as  determined  by  the  stag- 

nation line  noluLion  for  wall-j<‘ts  2 and  3. 
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FIGURE  31 

GEOMETRY  FOR  DETERMINATION  OF  MOMENTUM  FLUX  INTO  FOUNTAIN  BASE 

(STAGNATION  LINE)  SEGMENT  A-B 


TIu!  total  I'lux  of  momcrntum  unt«rlan  thu  nrliitrary  fountain  l>aHi*  minmant 
A-U  l«  Rivun  !iy 
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Tliu  oxprusH  1.011!)  cluvu.lopod  In  Rt-'rarmoti  11  for  tlio  mnnu-aitum  l lnx  onti'iMnn 
the  baHLiH  of  tlirea  and  finn*  jot  fountaJna  are.  nwro  tpHtr  let  Ivo  Hineo,  for  tho 
rase  of  non-vei't Ira  1 jot  Implngoment,  the  radial  wall-jet  momentum  flux  1h 
dopondont  on  tho  local  radius  R alnuit  the  Implngemont  point  ■ Tills  Is  Indl- 
ofUed  by  Kciuatlon  (4(1).  Thla  R dopondenre  apparently  iu'cchh Itated  tlie  aMHump 
tlon  of  a "dividing  Mtro.'imUne"  perpondlru far  to  tlu'  lino  r.onnortlnn  the  jot 
Imp  Lngomont  point  pair  for  the  fountain  hano  monii'iitum  flux  dotormln/it  Ion  In 
Rufofearo  II,  Thfnugli  the  use  of  l',(|uatlon  (fib),  tlio  ooord  Inaten  of  e.io.li 
point  In  Lho  Htagnatloii  lino  Hegment  In  tlu'  ground  plane  .■ii'i,'  not  rL't|ulred. 


the  angle  (|i  In  Equation  (66)  is  related  to  the  angle  i|)2’  througli  the  use  of 
Etniatii.m  (63)  applied  at  impingemont  point  ninnbor  three,  l.e. 

1 - Iy.,  + IHO”  - M :i6n°  - i.i.j  - . (67) 

The  fiUMmi  1 a 1 Ion  ol  the  momentum  flux  given  In  Hqu.illon  (66)  refurlrea  value, s 
of  pertliK'iU:  angU'H  In  the  loe,al  flowlMi'ld  only,  and  doea  i\ol  dtqiend  on  l.lie 
value  of  the  local  rndiuB  R about  the  local  implngument  point • This  la  a 
direct  result  of  the  indepondenoo  of  wall-jet  radial  momentum  from  local 
radius  about  the  jet  Impingement  point  (see  Figure  20). 
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4.  CALCULATION  OP  THE  INDUClSI)  PLOW  PTELD 

l'\ir  II  VTOL  r.imf iRuraLiDii  tii  linvur,  tliu  iiiiihh  entniinmont  of  lUo  frue-;) liih 
ami  wall-lorn  ImluL’an  an  uxlortml  flow  I'luUI  Hur rinindiuH  tho  airoraft,  Ho- 
caiiMi'  I 111'  Imliu'iHl  fli'hl  la  ofiHon  r la  1 1 y Ini'iimpri'HH  ILl  t>,  Inviai'l.d,  and  Irrola- 
llonal,  |niU'Uttal  flow  Lht'ory  can  be  employed  l.o  calcn.late  the  reanllanl 
forcen  and  monumls, 

'fhii  iiurtLnenl.  flow  buniularios  are  .1  Iluaf rat ed  In  l■’lv,nl•e  2.  In  accord- 
ance wlr.h  the  mathodoioglcH  of  the  procedlng  SBCtlonB,  the  geometry  and  masa 
entrainment  ratc.a  for  both  free  and  wall-iets  enn  be  enleulatudi  For  nppli- 
c.atlon  to  potential  fow  theory,  the  dlatrihiitlon  of  nmaa  entrainment  rates 
imivldea  presorlhed  normal  velocity  boundary  c.onditlons  on  the  Jot  aurfaecsi 
Slmnltnncona  boundary  conditions  for  the  aircraft  are  obtained  correaponding 
to  flow  tanRcncy  ror(ulremcnta  on  all  solid  aurfaceH.  The  complete  distribu- 
tion of  Noumaim  normal  velocfty  boundary  conditions  la  theoretically  auffl- 
cient  to  determine  the  entire  Indncud  flowfleld,  Incliullnp,  the  aircraft  prea- 
anre  dlstr Ihutioni 

The  DouRlaH  Neumann  Three-Dlmenalonal  Potential  Plow  Compiiter  Program 
(Keference  1ft)  was  solected  aa  the  fundumontul  tool  for  calculating  the  in- 
duced I lelili  Suction  ^itl  briefly  deacrlbes  the  program  ami  explains  the 
motivation  Cor  Its  aelectlon.  Program  modif lent  ions  which  tnilorud  Ulie  Dougina 
Ncunmnn  method  to  VTOL  appltcations  are  dnaerllnicl  In  HiictUm  s,2.  In  addl- 
I l.un  to  allowing  for  arbitrary  Nuniminn  boundary  conditlonH,  theae  modifiea- 
tlona  Himplify  user  Input  rmpilreiiuaits  and  provide  Incrcamid  accuracy  and 
tifflclunc.y  for  configuratlona  wltli  wall-jets.  Section  A.'l  proaunts  the  ro- 
aults  of  numerical  exporlmontatlon  conducted  to  establish  program  user  guidu- 
1 lues  tor  Jet-Induced  flow  problems, 

A (lotnlled  user's  manual  for  the  modified  program  luia  been  written  (.IL'Cer- 
unce  17),  The.  manual  contalna  complete  Input  iuHtructlons  for  llic  program 
))luH  HampJ.e  tnpiit/outpul.  for  a alm|)le  hut  11  luaurntlvc  cane. 

A . 1 The  llouglaa  Neumann  tiothocl 

Thu  Douglaa  Neumann  method  is  hiisod  on  tliu  fac.L  that  any  IncomprosHlbl  o 
potential  flow  velocity  field  can  be  considered  an  being  Inducccl  by  a partl- 
eular  surface  dlatrlbutlon  of  source,  aingularit  fes  on  all  flow  hotindarleH, 

For  example,  In  Jet  interaction  problems  the  aiipropr  l-iUe  bouiularleK  of  tlu' 
Invlucld  flowfleld  are  the  aircraft  surface,  free-Jet  aurface,  and  wall-Jct 
Hurl'ac.e.  The  Douglas  Neumann  aulutLon  philosophy  Involves  Llie  veproHent.at  1 on 
of  a conflnuouH  source  d Istr  llnitluu  over  the  bound. ir  1 es  It,'  a discrete  dlstrl- 
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button  of  (’.nuHtnnt  Btrongth  Hiuirfo  qiuidri  Intoral  ponnlN,  A coiiHUfuioiuui  of 
tiK'  1 Inonr  nnfutri'  of  potonttn!  iMow  Is  tlmf  IrlW!  rnmpipto  vc’loHty  nt  nny 
1'i.old  pol.i\t  Is  simply  tlio  1 tiu'iir  pomhliwii' 'on  ni  I ri'i'  hi  roam  vi'lm  lty  plus  i lu' 
sum  of  tlu’  vi'loi' 1 1 loH  Imlut'C'd  by  isu'h  Indlvldiuil  •loin-ci'  p.ini'l.  Furl  Irthioi-!', 
the  vcliu'U.y  liuluued  bv  a alnglo  (nuiri't'  iiaui'l  la  'lln  ctly  prnport  Iona  1 tu  tlu.' 
paiu'l  Hourei'  nLn.'ngtl»,  l.)i'nMi'(|iu'ii(.  I v , it  Is  pusulble  to  oHlal)!  lull  a uyslcni  ol 
liiioar  alp,ohraU'  oipial  Iona  tolal  itip,  prcai'f  I b"d  normal  voloi  lly  hoimilary  i on- 
ditlonH  at  dlHcrctu  cmitro.l  po.tntH  to  tiio  unknown  panul  Knurc.o  Htrungths, 
Solution  to  the  systom  of  oqiuUlonM  gi'HorutcH  tbu  nomplnte  Hourao  dlHtrlbution 
wliidi  1h  used  to  culoulnto  tlio  oomploto  vul.oulty  voctor  at  any  point  of  intor- 
oat  in  till'  potontini  field.  Uso  of  tlui  Hurnoul.li  equation  convartH  velocity 
magnltudo  to  local  pressure*  coofflcltmt  for  usti  In  force  and  moment  integra- 
tion. 


The  modeling  of  arbitrary  flow  boundarioH  by  a quadrilateral  panel  net- 
work is  illustrated  in  Figure  '12,  where  the  pannul  ed  houiularius  Include  the 
alrc.rafti  free-jots,  and  ground  wall-.jets.  The  panelled  geometry  is  dictated 
by  the  user,  who  must  Input  tlu*  corner  point  coordinates  of  each  panel.  Ks- 
cept  at  boundary  udges,  em'b  corner  point  is  common  l;o  four  adjacent  panels. 
The  I'onr  corner  points  which  define  a panel  noud  not  he  exactly  coplanar,  be- 
cause the  program  automat  teal ly  generates  the  flat  panel  having  the  least 
square  error  at  the  corners. 


FIGURE  32 

DOUGLAS  NEUMANN  PANEL  MODELING 
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Pnnoling  n conCigurutlon  1h  HcintBwhnt  of  mi  art  uml  »‘oqu:lros  knowlodRa  of 
bflalc  protiaduraX  ruXe«»  The  program  allowa  u maximum  of  one  thousand  panals, 
which  will  providi’  adaquatu  tu'curacy  for  typical  configuratlnnH , Panols 
with  largu  taiH')'  ratloH  or  high  aspoot  ration  Mhmild  ho  avuiUod  whore  pus- 
Hllih'.  (5roat('i‘  panoi  doiifiLty  should  In*  appl  Uul  In  n glunH  ol  Htrong  prea- 
Mure  I’radieiU,  typhmlly  wlioro  Hurface  eurvaturo  In  hlgheaL.  Thin  goomeLrlen 
require  larreauiHl  I'aiiel  ileiiHlty,  and  IL  la  important;  that  panel  corner  polntn 
lui  allgnod  on  oppositu  HurfaeoH  o!  a thin  geonie.Lry « 

fit'iectlon  of  the  Douglas  Mau.nmin  program  was  prodicatod  upon  a caruful 
trattooffof  advantagos  and  disadvantages.  A duslrabiu  foaturo  of  tho  program 
In  tlu!  fact  that  no  small  disturbance)  nssumptinns  arc  made,  Thuroforei,  tho 
prediction  acc.urncy  is  limited  only  by  the  number  of  panels;  as  this  niimbar 
increases  without  bounds,  tlie  calculated  numerical  solution  will  approach  the 
exact  analytical  potential  flow  solution.  The  wide  acceptance  and  oxtunslve 
c.lteckotit  by  both  industry  and  government  personnel  provide  a degree  of  con- 
fidence in  the  capabilities  and  coding  accuracy  of  tho  Dougins  Neumann 
method,  Presently,  an  automatic  geometry  paneling  routine  is  being  de- 
veloped by  Dougins  Aircraft  Company  under  contract  to  NASA  Langley  Re- 
search Center,  Upon  completion,  this  routine  can  be  coupled  directly 
with  the  existing  program  to  provide  n simplified  user-or lentod 
package. 

However,  there  ere  two  disadvantages  to  the  Dtniglas  Neumann  method  that 
require  recognition.  First,  tlio  computational  expense  tends  to  Increase 
wltli  the  cube  of  tlio  number  of  panels  N,  while  prediction  accuracy  only  In- 
creases linearly  with  N.  Therefore,  prediction  accuracy  tends  to  vary  as 
computational  I'xpense  to  tho  tlilrd  power.  Second,  numerical  instabllitlos 
can  he  generated  if  the  program  is  applied  to  thin  lifting  geometries  bo- 
cauMe  source  strength  magnitudes  become  excessively  largo.  Fortunately, 
work  is  in  progresH  to  alleviate  these  two  sliortcomings ■ 

'ITie  higlier  order  curves  panel  approacli  of  Hess  (kerc.rence  18)  will  gener- 
ate prediction  accuracies  tliaf  vary  as  computational  expense  to  power  two 
Instead  of  power  three,  and  this  approaeh  Is  presently  being  rormulated  and 
coded  for  the  three-dimensional  program  by  Douglas  Aircraft  Company  under 
cunt r act  to  NADC, 

The  tliLn  surface  anomaly  can  be  eliminated  by  incorporit  Ing  tlie  mild 
Hoiirce-iloiihlct  combination  eor  respond  ing  Lo  tlie  I'laHslcal  Creen'w  liliaitltv 
rel  at  loiisli  Ip  (Keferencu  19).  McDonnell  AlriTiift  Company  has  imullfled  the 
two-d  Imens  limal  Douglas  Neumann  method  by  Implementing  a ri'cnilly  devi'loped 
rirei'n's  hli'iitlty  formulation  (Reference  llO).  Typical  luimi’rical  remillu  ari' 
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praoenced  in  PlRura  33,  wliore  n suporci’illcal  alrfol]  in  incompressible  How 
and  A°  anfilo  of  attack  was  modeled  by  AO  panelsi  A conformnJ.  mapping  solu- 
tion by  the  method  of  Catlierall-SellH  (lleference  Zl)  provides  a virtually 
exact  two“d ImenHlonal  pressure  distribution  llir  purpoaes  of  comparlsout  The 
(lalheral  1 “Sel  1 s value  of  lift  coeff  li' lent  ((1^  » 0,903)  was  prescrlbecl  for 
eai'li  of  three  panel  mi'llunls  to  ))reviMit  the  peaslblllfy  of  arbitrary  Kult.a 
I'oniittlon  selection  c.louitlng  ilu’  comparison.  The  conventional  low  order 
Douglas  Neumann  solut  ion  Is  seen  Im  deviate  s I gn 1 1 leant;  1 y from  the  exact  solu- 
tion near  the  thltt  hlKhly  loaded  trailing  edge,  Uusortlng  to  the  higher 
order  curved  panel  appioaeh  with  the  Douglas  Neumann  program  Is  soon  to  im- 
prove, hut  not  eliminate  the  anomaly.  However,  using  the  Oreon's  identity 
approaeh  with  the  higher  order  curved  panels  virtually  ulimlntites  the  trailing 
edge  prohlem.  Furthermore,  the  tireen’s  Identity  solution  is  at  least  as 
accurate  as  each  of  the  other  two  panel  solutions  ovur  the  entire  remainder 
of  the  airfoil  Hurfnee,  MclX^nnel I Aircraft  Company  is  presently  under  cun- 
trnet  to  NASA  Langley  Research  Center  to  modify  the  three-dimensional  Douglas 
Neumann  pn  gram  to  accommodate  the  (Ireon's  Identity  formulation, 


t)»iu  niKi  •« 


FIGURE  33 

SUPERCRITICAL  AIRFOIL  PRESSURE  DISTRIBUTION 


IL  Is  /ippai'ent  that  the  shorti  omlngH  of  the  Douglas  Niuimnmi  thtei'-d  I meii- 
slunnl  program  should  he  luilanced  agalusi  the  fuel,  that  s 1 c.u  I f leant  Improve- 
ments to  the  program  /ire  currently  under  ilevel  opmcnl  , II  Is  I hui-el'ore  rc/ison- 
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iiblo  to  Holuc't  the  oonvenulorm  1.  DouhIgw  Neumann  program  ns  the  hnslc  poten- 
tial flow  tool  for  this  Htutly,  with  tho  knowlo<lw,n  timt  In  'lui  futura,  the 

Impi'ovoti.tints  being  made!  under  other  independent  studies  can  be  merged  to  form 
ont'  lilglily  flexible,  efficient,  and  nccuratu  luuciiLial  flow  i;umputer  program, 
h , 2 Mod  I f I I'll!  Ions  to  botiglah  Neunumn  l‘riigram 

In  (U'der  to  |■,'U'1 1 Itnte  the  application  of  the  threi—tl  Imens  Jonal  Douglaa 
Neumami  pro  glam  to  typical  v I ri’ raft  / jnt/groimd  inter/ict  Ion  prubie.mH,  the  five 
pl^'!’,l^lm  s'.ioil  I f 1 c,i(  I un<)  to  he  illacuHaed  helow  have  hi'cii  lni)i  I laiien  t ed  and  lully 
chi  ,'l.cil  oul  . TIu'ae  mod  I I 1 f.-|  I I oaa  impoHc  niliiiinal  vniikload  on  the  pni)',  I'aiii 
user  an  I di'  not  impali'  any  existing  program  capab  M I ties. 

Option  to  Allow  Arbitrary  Nenmaini  Houndary  Clniullttons  - In  Llie 
fo ntiul at  Ion  and  coding  for  tlie  uxlatLng  Douglat,  Neumann  program  in  Keference 
111,  f low  tangency  boundary  conditions  urn  assumutl  to  apply  on  ail  boundary 
HurfaceH.  Since  prescviberl  entrainment  normal,  velocities  on  jet  boundar.les 
I'i.olate  tills  rustc  Let' Ion , It  w.;s  necessary  to  modify  tlie  program  to  aecept  an 
arbitrary  pfeacrlhed  vaJuo  of  nornnl  velotlcy  as  tlie  boundary  condition  at 
each  pane,  control,  potiu.  The  r-.'-Jlflcat  Ion  was  ac'coiiipl  Ishid  In  a faslilnn  such 
that  the  icier  neeu  pi.'eHci'Um  valuoH  only  on  genmetrie  HtailloiiH  He.lecteil  by 
him  and  the  program  will  assume  Clow  tungiancy  elsovdiere,  '''or  example,  on  solid 
aircraft  tairfaces  no  normal  voloclty  values  need  he  input  since  the  program 
will  assign  a xcro  valuo  on  those  panels  automatically, 

O|)tl.on  to  Allow  Wake-rreo  tlonf Igiirat i oiii  - Tlai  exlstln};  pragnim 
always  assumes  tho  prusuiico  of  a tnilllng  vortex  v.ake,  character intlc  of  the 
flowfleld  Induced  hy  an  aircraft  lii  forward  fllglit,  For  VTOI,  aircraft  In 
hover,  the  effects  ol  wind  gusts  are  usually  no  mor.'  slgnincant  than  the 
lot- Induced  flowfleld,  thereby  Invalidating  tho  concept  of  a hr  riling  wake, 

The  program  was  therefore  modiflod  to  accept  ellhei'  the  presence  or  abRunce 
of  a trailing  vmtex  sheet  at  the  dlscrctlou  ol  the  user, 

4,2,1  Simplified  (ienmutry  ModlFlcntlon  Option  - The  ex  Is  ling  program  re- 
ciulred  geometric  coordinates  to  lie  Input  at  tlie  raU'  of  two  iiolntn  per  card, 
Therefore,  If  tlie  user  dei'iLrod  to  rutke  a local  Ir.ed  gconietry  modification  he 
was  forced  to  add  nr  delete  an  even  number  of  points  or  else  ffoe  the  pros|'vct 
of  rcpuuchlng  Input  cards  for  every  mihsetiuent  ooliu  (possli>l  ' several  Lluuis- 
aml  values),  Tlie  program  was  thoruroro  modified  to  .lilow  an  .irhltrary  mix- 
ture of  cards  each  having  coordinates  for  elllier  one  or  i.wo  points,  The 
I modified  program  nov  assumes  each  card  to  have  two  |ailnts  unless  a non  lU'ro 

value  Is  placed  In  card  column  70,  at  which  t.ime  tlie  program  automatically 
ai'cepLs  Just  Che  tleslred  single  point  coordinates  on  I he  cai'il. 
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Cciordlnfltt'  TrniiHl'orimiltcin  Out  inn  - Thn  oxiHtlnK  progprtiii  anniiincd 
jLI  geoinctrtc  panel.  coordinateH  t.n  bo  Input  in  a unmmoii  C«rton.lan  coordlnato 
HyMtont.  'niortJrore,  If  cho  itHor  dpulrnti  tn  rhnngo  an  aircraft  attitude  «r 
local  Inn  witli  roapoct  to  a gmuiul  piano,  every  aLrcraft  coi)rdl.nato  had  to  ho 
roconiputial  by  hand  and  ropuvu'lu'd  on  carda.  The  o.oordtnaU’  t ran.arormat Jon 
optlni  I'uablou  Iho  iiHur  to  arbitrarily  move  an  liullvidual  nooniotrv  .aoclloa 
Much  an  a wing,  lunolago,  or  )oi  wliliout:  fllsturhlng  t.ho  rcmalmloi'  of  the 
ni'oim.' I rv . for  oxainpio,  the  roll,  pitch,  and  hclp.hi  of  the  alia  rail  In  FI)', are 
32  could  aaally  be  changed  with  roapoet  to  the  ground  plane. 

The  tvana format  Ion  procodnre  conHlatH  of  rotation  first,  translntJon  sec- 
oiul.  If  tlie  panel  Input  coordinates  of  a .section  are  in  tlie  (x^,  V|,  /!|)  sys- 
tem of  Figure  34,  the  program  rotates  the  coordlnatea  through  Kuler  angles  i', 
tlien  0,  and  finally  li,  whore  these  angles  are  arbitrary  Input  vnlMeH,  Next, 
tiu!  rocatod  coordinates  are  translated  by  arbitrary  Input  values  Ax,  Ay,  and  Aa. 
The  rormulaH  relating  the  Input  coord  Irmtes  (X|,  y^ , to  the  final  trans- 
formed eoordlnateH  (x,  y,  x)  follow. 


FIGURE  34 

COORDINATE  TRANSFORMATION  OPTION;  ROTATION 


X 

» Ax 

+ 

P,,Xi 

+ 

rn^yi 

+ 

"l‘‘l 

(68a) 

y 

■ Ay 

+ 

+ 

m^V] 

+ 

n,,x^ 

(68b) 

" At 

+ 

f.jXi 

+ 

m,,yi 

+ 

n,,n, 

(68c) 
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The  caordlnaie  iranHformutlon  option  is  of  Limited  valno  for  conventional, 
aerodynamic  problwmH,  but  in  eHBentln.1  for  V/STOL  aerodynamicH  l.n  which  roll, 
pitc.li,  and  ground  height  are  commonly  varied  parameters! 

4,2,.')  Modified  Method  of  ImaaeB  Option  - For  tlic  Holiition  to  aircraft- 
let  ground  interaction  problems,  tlie  application  of  a method  of  ImanoB  con- 
cept nenerateH  siRnificanU  improvomentB  in  both  computational  efficiency  and 
predict!  m accuracy,  Thia  concept  has  been  implemented  In  the  liouRl.as  Neumann 
ProRram  an  an  option  dcBlunated  "Modified  Method  of  ImaRCH",  Subject  to  the 
mild  restriction  that  all  wall-jet  pnnela  lie  in  a c.oinr,wn  plane,  the  modified 
method  of  iniaROH  approach  cllminateu  ns  unknownH  the  sonree  strength  de.nHltlew 
on  the  wall  panclM,  HecauBu  the  computlnR  time  required  to  solve  a Hystem  of 
linear  equatJoiiH  1h  proportional  to  the  cube  of  the  number  of  unknown.'^,  even 
a small  reduction  in  unknowns  will  bo  fllRnlf .Leant  tor  Uliu  larRu  numbot  of 
panels  aHBoc.latcd  with  V/8T0V,  aerodynamicH , 

The  basic  Roomotry  panelllnR  for  the  modified  method  of  imaRos  1h  illus- 
trated in  FlRure  Xl,  In  addition  to  the  aircraft,  free-jet,  and  wall-jet 
panelllnR,  the  program  will  uuLomaticnlly  panel  an  imoRG  aircraft  and  im.age 
free-jet  Upon  activation  of  the  modified  method  of  imanes  option.  The  pre- 
Benc.p  of  the  image  panels  vlolateB  no  theoretical  restvlc tions  bucauHe  they 
do  not  1.1.0  within  the.  actual  invlscld  flowfield.  For  ciicli  unknown  soure.e 
density  on  an  aircraft  or  free-jet  panel,  tliere  I.h  an  equal  strength  souree 
density  on  the  correspondinR  imaRc  panel,  liy  symmetry,  the  normal  velocity 
indiu'ed  by  tlie  aircraft  and  free-jet  panels  plus  imaj'.e  panels  at  any  iiolnl 
on  tile  wall-jet  will  be  identically  zero.  Therefore,  the  net  normal  velocity 
on  till'  wall-jet  Is  attributable  only  lo  the  freestream  component  plus  t'le 
component  induced  by  the  wall-jet  source  distributions.  Under  t lie  res tr Let  Ion 
that  the  wall-jet  panels  are  coplanar,  this  latter  component,  can  he  sliowii  to 
be  exactly  twice  the  local  source  density.  .Since  tlie  net  normal  veloclly  on 
the  wall-jet  Is  prescribed  as  a boundary  condition  from  the  known  wall-jct 
ent  ra  fniiuint  velocities,  It  Is  possllile  lo  di'termlne  tlie  source  denalty  din- 


Airplane  Vivj  0 


FIGURE  35 

SAMPLE  PANELLED  GEOMETRY  FOR  THE  MODIFIED 
METHOD  OF  IMAGES  OPTION 

Airplane  Free  Jet  - Wall  Jet 

trLbatlon  on  the  wall-jet  panels  ininitHliafely  wl  Llimit  sol  .’lnr,  I’or  alreral't; 
and  free-Jot  source  densities.  riierefnre,  the  vcliuilty  I li'.KJ  Induced  by 
the  wall-jet  can  be  computed  beforehand.  The  only  reinalnlnK  unknowns  to 
be  solved  are  the  source  densities  on  the  aircraft  panels  and  I'ree-Jet 
panels  where  the  known  entrainment  velocities  •'  ■■■I  of  themselves  deter- 

mine the  source  densities,  The  o'Ctra  expensi.'  recjulri-d  In,  calc',\la  ,ing  the 
image,  panel  .Influence  coefficients  Is  typically  1 nsign It  1 can L because  tl' 
almplu  flowflekl  formulas  of  Hess  (He  ferenee  Ifi)  .are  employed  wheiu‘Ver  tae 
distance  to  ;in  Influenced  point  Is  suftlciently  large  compared  to  the  In- 
fluencing panel  d Imen.'i  1 iins . A seeuiidary  advantage  ol  I he  modified  n'a'thod 
of  (mages  Is  that  the  boundary  conili  lions  will  be  moi'e  nearly  satisfied 
between  wall-)et  jiaiiel  control  point;;  than  In  the  conv'mt  1 ona  1 solution 
approach,  which  Is  a coiiseciuence  of  the  f;:-'.!.  th.it  symmetry  conditions  ap|ilv 
at  all  points  tn  the  w.ill-Jel  Inste.ad  of  control  points  only. 

lU'ianise  ihi'  progr.am  .allows  more  I'n.aii  [UU'  symmetl'v  plane,  Lhe  mod  i I' led 
ineth.od  of  Im.ages  option  doe.s  |■'ot  foree  tin-  lun  r to  paiiol  both  the  right,  anil 
I c f t-hand-H  I de.s  of  .a  symmetrleal  alrerall  above  a \,ial  1-Jet 
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^ ^ Numorlr.il I Exporlmp.nta  with  PouRlas  Neunmnn  ProRrnni 

lU-u'.auHu  pruviouH  (ipp;l i (,'atioii8  of  the  Dmip,inH  Noiimann  prt'grnm  to  Jet  Inter- 
artitm  probli'inH  have  been  limited,  it  was  oonBicleracl  benorioial  to  conduct  a 
Munnillvlty  study  uf  Jot  panel  model  1 ok  el'Fects  on  luilculatod  solutions,  fiiudi 
a stuily  would  re.voiil  liasie  user  panolllnn  rules  iind  Indicate  the  proper 
ilej’rei'  ol’  pre.d  lot  ion  accuriu.y  versus  paiiel  density. 

In  .icoordanei'  with  thi'  philosophy  that  the  slnii'lest  posslhle  case  should 
lie  used  to  Isolate  numerical  t riaids,  a systematic  panel,  stutly  was  I'onducted 
lor  an  ax  Isymmet  r Ic  Jet  exhaustiu)*  I'roiii  an  Infinite  flat  wal  I.  This  simple 
case  is  especially  attractive  because  both  experimental  data  and  an  alternate 
exiu.'t  potential  flow  solution  are  available  lor  the  wa.l]  pressure  distribu- 
tion (Reference  22), 

The  basic  panel  model  is  illuatrtitod  in  Figure  16,  Rntrainment  normal 
velocitiei;  consistent  with  those  of  Reference  22  are  prescribed  on  the  Jet,  and 
flow  tantioncy  Is  prescribed  iilonp,  the  flat  wall,  The  entrainment  normal  velo- 
city distribution,  From  Reference  22,  is  V.,  “ -.012  V for  C)  < 2 < 6 1)  and 

N ,1  e " J 

V.,=  -.114  V,  for  Z > 6 D,.  The  pressure  coefficient  Is  defined  as  the  neRa- 
N 1 e J 

tive  of  the  square,  of  the  rutii)  of  local  spt'cd  to 

Thu  objective  of  the  study  was  to  determine  the  sensitivity  of  the  ciilcu- 
latecl  wal  l pressure  distrllniu  Ion  to  systematic  chaiiRc.'s  In  the  basic  piincl 
nHulcllnu  of  Flpure  16.  la  that  ruRard,  five  se.pataLe.  panel  modclltiRH  were 
established,  The  first  model  liiR  was  Rcnerated  l)V  riimovlnp,  the  Wiill  pam.il  s of 
the  basic  setup  of  Figure  16  and  then  addliiR  an  image  Jet  In  ui.>..ordanc.e  with 
the  metliod  of  ImaRos  (Figure  .ihi).  Illustrated  In  Flgvire  17b  iMrc  the  second, 
Ihlrd,  and  fourth  mndcl.s.  The  second  moilel  Is  simply  the  ha.sic  semp  of 
Figure  16.  To  demonstrate  the  effect  of  reduced  panel  density,  the  third 
nu.'del  was  created  from  tlio  basic  setup  by  a fifty  peicent  reduction  in  panel 
density  both  on  the  Jet  In  the.  axial  direction  and  on  the  wall  in  the  radial 
direction  away  from  the  Jet  axis.  The  fourth  model  was  established  by  extend- 
ing  the.  flat  wall  of  the  basic  setup  across  the  Jet  exit,  tlicrehy  effectively 
c.ipplng  the  Jet.  A total  of  eighty  panels  (twinity  per  quarter  vli'w)  were 
used  to  cap  '.he.  Jet,  and  flow  tangeucy  boundary  conditions  wt'fe  applleil  to 
I'.'ich  such  panel  on  the  side  f.'ic..ing  the  Jet  (x  0 ) , The  ftlfh  model  repre- 
sents a let  emanating  from  a wing  of  small,  hut  finite  thic.kncss  (Figure  'J7c,). 
The  panel  trig  was  established  by  adding,  an  upper  mir'ac.i'  wall  to  tlu'  h.isic 
setup.  The  wall  upper  aurfac.e  panel  i.'orner  polnt.s  project  onti’  the  lower 
Hurl'acr  coiaikO'  points,  in  accordance  with  a ha8.1c  panel.  niU'  of  lle.ss,  I’low 
tangency  was  prescribed  on  both  the  upper  and  lower  .auri'ace  of  tin.'  wall,  la 
order  to  c.lwi  lactcr  l>;c  a wlag-type  configuration, 


(a)  Jat  > Imiga  Jat 


FIGURE  37 

THREE  PANEL  MODELS  OF  AN  AXI8YMMETRIC  JET  EXHAUSTING 
FROM  A FLAT  WALL 


Colculatocl  resaltH  for  the  wall  pressvtre  and  .source  tlensity  diatribu- 
tl.ons  appear  In  FiRures  38  tliraugh  Al,  0,f  the  f,Lve  panel  models,  the  one 
ine.urporat Lng  the  metliod  of  Imagea  1b  of  nocoBalty  the  most  accurate  because 
symmetry  guarantees  tlie  absence  of  leakage  anywhere  on  tlie  wall  plane.  A 
comparison  (Figure  38)  of  the  wall  presauro  diatrlbution  cfilculated  by  tlio 
llouglus  Neumann  program  using  the  Jet  plus  image  Jet  paneling  (Figure  37a) 
demonstrates  excellent  agreement  with  both  experimental  and  analytical  values 
from  Rerercnc.c  22.  The  analytic  curve  was  obtained  by  assuming  the  Jet  to 
be  a concentrated  .source,  filament  and  then  applying  tlie  method  of  images. 
Virtually  all  the  discrepancy  between  the  analytic  curve  and  the  c:alculnted 
DoiigloB  Neumann  values  l.s  attributed  to  the  upproximation.s  associated  with 
a finite  panel  density. 

Krrors  nasocinted  witli  the  introduction  of  the  fiat  wall,  panels  are 
evident  in  I'lgure  39  for  both  tlie  dense  and  spatae  paneling  of  Figure  37h. 
Calculated  pressures  are  quite  good  except  lor  the  five,  or  six  panels  closest 
to  the  Jut.  It  is  noted  that  within  this  region,  the  .sparse  panels  Intro- 
duce. approxinuituly  twice,  the  relative,  error  of  the  densu  panels,  which  is  a 


conscquuiun  of  the  fact  tliat  the  low  order  approach  of  the  existing,  pn'j'.rnm 
gencr;it;cM  prediction  accuracy  tliat  varies  as  the  first,  power  of  panel  ilcn- 
slty.  However,  the  simple  technique  of  capping  the  Jet  eliminatea  most  of 
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Dlitanca  Along  Will  • 2r/Dj, 

FIGURE  38 

PRESSURE  DISTRIBUTION  DUE  TO  AN  AXISYMMETRIC  JET 


Dlitance  Along  Wall  ■ 2r/Djg  gomoim  i 

FIGURE  39 

PRESSURE  DISTRIBUTION  DUE  TO  AN  AXISYMMETRIC  JET  CALCULATED  BY 
THE  DOUGLAS NEUMMANN  PROGRAM 
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the  problem,  but  at  the  expense  of  adding  more  panels.  Corresponding  to  the 
calculated  pressure  distribution  of  Figure  39  are  the  calculated  wall  source 
distributions  of  Figure  40.  The  close  agreement  between  the  source  distri- 
butions is  expected  for  a flat  wall'.  The' disagfeement  in  pressure  distributions 
la  alnipiy  attributable  to  tlio  neglect  of  source  gradients  on  a panel,  wlilch 
la  obvioualy  of  inereaalng  importance  near  the  edge  of  the  Jet. 

I'lgui'e  41  cIcminstrati'M  llial  tlu-  mldltlon  <if  i.lu'  wall  upper  surlare 
thickness  impairs  tlie  calculated  pressure  distribution  within  approximately 
five  panel  widths  of  the  Jet  edge.  Again,  however,  the  results  are  good 
oiitslde  this  range. 

It  is  concluded  that  i 

(1)  The  method  of  images  should  be  used  where  possible, 

(2)  Panel  density  should  Increase  near  the  Juncture  of  a Jet  and  solid 
surface. 

(3)  Prediction  accuracy  Is  proportional  to  panel  density. 

(4)  Wings  with  thickness  require  greater  panel  dena Lty  than  Infinite 
flat  wal Is . 

(5)  Jets  should  be  capped  with  additional,  panels. 

A complete  user's  manual  for  the  modified  Douglnp  Neumann  program  is 
available  (Reference  17),  In  accordance  with  the  iiser's  mmuial  instructions, 
a four  lift-jet  VTOL  attack  aircraft  in  hover  was  panelled  for  the  program. 

The  set-up  incorporated  approximately  one  thousund  panels  to  model  the  left- 
liand-side  of  the  aircraf t/free-Jet/wall-Jet  configuration.  Approximately 
three  man-weeks  (120  hours)  were  required  to  complete  the  panelling.  On 
the  IBM  360,  the  computing  time  for  one  complete  solution  was  fifteen  mimUes, 
Several  values  of  fren-strenm  or  cross-wind  velocity  can  be  computed  for  a 
fixed  geometry  at  no  significant  increase  in  computing  time. 
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5,  FOUNTAIN  ITPWASH  FLOW  CHARACTEIVISTIOS 
Tl>o  study  of  fountnin  flow  chorHcterist los  has  concentVfitud  on  oases 
wl\ero  tl\o  Jets  are  clrrular,  perpendlo.u.lar  to  the  Htound  piano,  and  of  equal 
thrust.  In  addltio\i,  the  pattern  in  whioh  tlie  nor.Klos  are  .irranKod  has  been 
usually  asHiinunl  to  pi>siH>ss  at  least  one  pj.ine  of  symnietry,  Tho  oxlsrcnce  of 
one  or  more,  planes  of  symmetry  slmpllffen  the  I'lowfiekl  poometry  eons  Iderahl  y , 
and  leads  to  a bi'tter  undi'rstand  l.ni>  of  the  fundamen  t.a  1 fountain  upwash  jiroper- 
tiuH.  For  a uonfinurutiun  wLth  two  perpend  U'ular  ,|ets  of  e(|ua  1 thrust,  eerlaln 
asaumpt-tons  about  tht:  flow  tunvinB  proco.ss  result  In  an  analoRy  butween  tho 
fountain  and  a Jet,  Basic  features  of  this  unalDgy  have  been  verified  in  a 
test  proRtam  conducted  in  a parallel  MCAIR  IRAll  effort.  Based  up(5n  the  two- 
Jet  results,  the  annloRy  between  fountains  and  Jets  bus  been  extended  to  con- 
flRurations  with  three  and  four  perpendicular  Jets  of  equal  thrust. 

An  initial,  qualitative  description  of  the  upwash  in  asymmetric.  two-Jet 
fiumtnlns  has  been  developed,  and  a plausible  onuloRy  between  the  fountnin 
and  a conical  flow  has  been  postulated.  A q\iantitative  description  of  the  flow- 
field  is  not  possible  at  present,  however,  due  to  a lack  of  adequate  data  for 
this  type  of  three  dimens ion.al,  fountain. 

5 . 1 Fguiviilence  Between  Symmetric  'rwo-..lct  Fniintains  and  ,Had I 

The  simplest  possible  fountain  flow  is  that  Renerated  by  the  Jet  arranRement 
shown  in  FlRure  A2.  Two  circular  Jets  of  equal  dJamutur  and  thrusl  ImpJnRe  at 
90*  to  tho  ground  plane.  The  Jets  are  separated  by  a distance  R and  located  at  an 
equal  lu'lflht  11  from  the  ground.  Tlio  waU-Jets  fk'w  radially  outward  from  tlie  Impingi- 
ment  point:  hel.ow  each  vertical.  Jet,  and  separate  nloug  the  line  AB  fa  form  che 
fountain.  Because  of  symmetry,  the  separation  line  AB  is  straight  and  is  located 
halfway  between  the  Jet  Impingement  points.  The  entire  flowflold,  IncludlnR  the 
fountain,  is  symmutric  about  plane  u (and  also  about  tho  plane.  containinR  tac 
JoL  centerlines).  Hie  radial  lines  emunutlng  from  the  ImpinRetnent  points 
represent  the  direction  of  flow  In  the  nxlsymmctrlc  waU-Jets.  Oil  flow  studies 
of  Jet  impliiRement  patterns  indicate  that  as  the  separation  line  .\B  is  appro.icbed, 
the  oil  flow  lines  curve  and  become  tangent  to  this  line,  but  the  region  of  curva- 
ture Is  limited  to  the  Immediate  neighborhood  of  the  line.  At  the  line  AB,  the 
w;ill“Jcts  turn  and  leave  the  surface,  forming  tlu'  fountain.  The  turning  proci'ss 

Is  idealised  .as  depicted  in  the  enlargement  cif  the  flow  region  around  point  f 

y 

shown  in  Figure  42.  The  vectors  1)^  and  l).,  are  veloulUes  in  the  wall-jet  profiles 
corresponding  to  Jets  1 and  2,  respectively.  Those  vectors  have  normal  (subscript  n) 
.and  tangential  (suhsi.ript  t)  components  to  the  separation  line  All.  It  Is  assumed 
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FIGURE  42 

IDEALIZED  WALL  JETS  AND  FOUNTAIN  TURNING 
Two  Symmetric  Jets 


that  in  tlu-  turniiiK  process,  the  normal  velocLfy  oomponontn  Ui^  and  U'j,,  are 
rotated  to  bo  pt’i‘pc.i\d Icul.ar  to  the  Rnnind  piano,  and  tho  lannonttal  ooniponontH 
Uit  and  IJ2(.  are  unchanged.  Thus,  the  angle  0 as  defined  In  Figure  A2  is 
preserved  in  the  turnitlR  prooosH.  With  tho  I'Xcoptlon  of  the  Immediate  nelgh- 
borliood  of  tlie  stangatlon  line,  to  a first  approximation  the  statlo  priuisure  In 
tho  fountain  is  taken  to  ho  constant.  ConsequeiiL  1 y , streamllnoH  in  tho  fouiUntn 
are  straight,  and  tho  fountain  velocity  profile  remains  at  an  angle  0 to  thu 
Vertical.  The.  ungLo  ^ is  different  for  different  points  along  All,  and  the  flow  is 
t'l'iuivalont  to  what  would  be  obtained  by  folding  the  ground  piano  along  the  separa- 
tion line  AIJ,  which  causes  the  impingement  points  for  li.’ts  1 and  2 to  rotate  by 
^)0”.  This  procoHH  la  lliustratod  in  Figure  'ti.  The  flow  in  the  fountain  would 
then  bo  analogous  to  a radial  flow  originating  from  a virtual  origin  locatLal  at 
a distance  S/2  below  tbo  ground  piano.  A radial  flow  sneh  as  this  mav  be  generated 
by  two  opposing  axisymmetrlc  Jets  which  Impinge  upon  each  other.  I’lows  of  tlil.a 
type  have  been  extensively  studied  in  Rcferenee  H, 

In  order  to  verify  tin:  iio.stulated  analogy  hetwei!a  a symmeirical  two-|ei 
fountain  and  tiu'  flow  in  a radl.'il  Jet,  a test  was  conducted  as  part  of  a MCAIR 
IRAH  program  (Ri'ference  211),  Properties  of  llie  fouutlaa  flowfleld  were  Investi- 
gated by  means  of  flow  visual  1 xat lou  and  flowfU'ld  surveys  with  a total  pressure 
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FIGURE  43 

EQUIVALENCE  BETWEEN  TWO-JET  FOUNTAIN  AND  RADIAL  JET 

niko.  Tilt'  roHirlts  nry  compared  to  clntn  obtuiiu’d  from  proHHuro  HiirveyN  of  nn 
.'ixlsymmutr I.C  rudlnl  Jut  cretited  by  thu  ImplnKemuiU'  of  two  opposed,  circular  Jots. 
TTui  apparatus  conslNtod  of  a pair  of  1.05  In.  ID  pipos,  siippUeri  witli  hlt>h  prossiire 
air.  For  tliu  fountain  studlua,  the  pipes  wore  arransud  ho  that  tliay  implnsc'd 
at  90“  to  a ni’'>und  board,  as  indicated  in  FlKurc  Ono  value  of  H ^ 2.0d|^^ 

was  used,  and  three  values  of 

S/d^^^  - 3.81,  8.  12 

wore  toHtod.  Tlu'  total  pressure  in  the  pipes  wn«  varied  between 

' ' bop  J'. 

An  oil  film  technique,  woh  used  for  flow  visualization.  A total  pruHSuro  rako, 
located  perpendicularly  to  the  Hymraetry  plane  a In  Fiiture  A2  wan  used  to  Hurvey 
thu  fountain  velocity  distribution.  Tlu*  probuH  were  aliftnud  with  the  local 
fountain  velocity,  UHluf;  Che  flow  visualization  rcHults.  To  ^!eneratu  a radial 
Jot,  thu  pipes  were  aligned  end-to-end,  so  that  the  jetH  Impinged  on  each  other. 
One  v.iluu  of  the  Jut  spaclng-to-d Inmetur  ratio 

Ax/d,  ■ 0.5 
Je 
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Wrtrt  tuHted,  Total  pipo  prunHuroH  of;  '10  ami  40  pal.  worn  Lestod.  Total 
prossurc  aurvoys  wc?ro  rondui'tod  with  thn  rnko  iicrpandloulnr  to  the  plane  of 
tho  Jot,  and  tho  prohoa  were  all.Knud  with  the  radial  voloclty  voolor. 

(Hoo  PlRori'  41.)  Tho  outiuit  of  tho  total  iiroHHiii'i'  rake  was  road  on  water 
iiianomotorn,  and  the  pipe  total  proHHiiroH  wiTi'  roail  on  mercury  manometera . 

I’ll’uro  44  liulloatoa  tho  loat  ar ranj'oraont  lined  for  I'lumtaln  flow  visual  lua- 
tloii,  T1k‘  plane  of  aynimot  ry  hetwei'ii  the  leta  was  replaeed  by  a solid  hoard, 
porpondlcu lar  to  tho  Krmind  board.  A Hheet  ot  paper  waa  attaohed  to  the  Hround 
and  Hynimetry  plane  aurfaecH.  An  oil  film  waa  bruahad  on  the  paper,  and  the 
Jet  waH  then  turned  on.  The  remiltnnt  oil  flow  )uittern  1.h  ahown  on  the  lower 
part  of  Figure  44,  (after  unfoldlnn  the  aluiet  of  paper).  It  1h  evident  that  tho 
oil  flow  lliiBH  on  the  plane  of  nymmetry  are  HtrnlRht  extcnHlons  of  the  oil  flow 
llneH  on  the  Rround  plnne,  confirming  tho  flow  turning  procoBS  deaerlbed  In 
I'lgure  42.  The  oil  flow  does  Indicate  a region  of  aeparatlon  right  nt  the  corner 
between  the  two  plnnea,  but  this  region  1h  very  small.  Figure  44  shows  the 
results  for  S/d|  ■ H.  Figures  45  n and  b indicate  a similar  behavior  for 


S/dj  - 12  and  6, 

respectively.  The  photographs  show  that  oil  flow  lines  on  the  ground  plane  are 
slightly  curved.  This  is  due  to  tlio  fact  that  In  tlie  experimental  setup,  tlie 
ground  pinna  was  located  vertically. 

Tho  oil  flow  data  Ind  Icated  tluit  there  Is  a qualitative  equivalence  l.ctwi>en 
the  flow  in  ii  symmetric , two-Jet  fountain,  and  a radial  flow  I'rlglnatlng  from  a 
virtual  origin  located  at  a distnace  S/2  below  the  ground  plane.  Measured  velocity 
prof  11  OH  will  now  be  compared  to  quantitatively  verify  the  analogy  hetween  tho 
fountain  and  a radial  Jet. 

The  first  question  to  bo  investigated  was  wlioLber  the  fountain  flow  was 
axlsymmetric  about  the  virtual  origin.  For  this  purpose,  fountain  velocity  pro- 
files were  mo, 'inured  nt  constant  values  of  the  radius  K from  the  virtual  origin, 
and  different  values  of  0,  ns  indicated  by  the  sketch  In  l'’lgure  4(i.  The  plots 
in  Figure  46  are  fountain  velocity  profiles,  normalised  by  tlie  Jet  exit  velocity 
Vj^,  .-md  plotted  vs  X (see  Figure  4:!).  Kach  plot  corrcHpomlu  to  a different 
value  of  K,  The  cuvvtis  Indicate  that  fer  a fixed  value  of  K,  profiles  at  0 »«  H“ 
and  20®  are  In  good  agreement.  However,  measurements  at  0 “ 40°  show  euiis I si ent 1 y 
higluT  velocity  profiles.  This  would  Indicate  that  the  foantaln  flow  Is  not 
axl  symmet  r le  about  tlu’  virtual  origin.  The  result  Is  somewh.it  1 neonc  1 ns  f ve , luiwi'ver, 
because  all  flow  patterns  on  the  ground  ho.ard  Imllcale  I hat  the  Jets  |iI',mIui' I nj'  tlie 
rnunt.'ilii  ,'iri'  net  quite  axlsvmmet  rl  r. 
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FIGURE  45 

WALL  JET  AND  FOUNTAIN  FLOW  VISUALIZATION 
AT  Dl FFERENT  VALUES  OF  JET  SPACING 
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n^ure  !\1  tiompa rt'H  thu  fuuntaJu  aiul  racllal  Ici  vnlni'lty  profile;'  aliaiios 
at  si'Vor;il  vaUios  of  Uto  d LMtanoo  K from  llic  fountain  vlrLvia!  nrl>>ln,  i'<r  from 
tin'  radial  axin  of  aymnutry.  Tlu>  vi' loc  1 1 Ins  aro  plottod  In  a l.mi  1 ;■  r I ty 

vari.ibli'M.  Tliat  Is,  at  anv  p.lvoii  valiio  of  K tho  volorltlrs  aro  normal  Izod  bv 
tliL'  maxliiuim  voloolly  at  that  valui'  of  K,  and  the  dlstanoi'  x Is  normal  ixod  bv 
tho  dlstanoi'  ai  whli.'h  tho  voloi'lty  has  dooayod  lo  o'S'  half  of  Its  inaxlmum 

valuo.  The  data  show  that  the  velocity  profllus  aro  approximately  sel f-s tmll ar , 

Htnee  correlation  1h  achieved  for  the  fount,aln  and  Jet  i)roflleH,  for  different 
values  of  U and  values  of  0 up  to  20“,  The  dotted  curves  in  Figure  A7  hnvc  been 

obtained  from  a fuirinB  of  the  correlated  radial  Jot  data  In  Reference  8.  The 

plots  in  KlRure  A7  then  indiiu'ite  that  the  iiiensurod  fountain  velocity  profile 
Is  In  Bo'd  nureement  with  that  for  a radial,  axisymmet rlc  Jet.  Sc!  1 E-slmllarlty 
of  the  foontaln  velocltv  profileu  Implies  that  its  velocity  field  may  be  speci- 
fied by  tic.-  hal  f“vuloolty  width  ''tid  the  local  peak  velocity,  Vm.  In  addition 

to  the  profile  shape. 

figure  A8  shows  the.  dependc'nce  of  J<]y2  R,  for  tho  radial 

jel:  and  for  the  fountain.  The  fountain  d.ata  show  that  linenrly 

with  R.  Independently  of  tlie  Jet  spacing  8/d^^^  and  the  pipe  pressure  Ftu'  the 
radial.  Jet,  Xj^..  also  nrows  lliiiu'irly  and  is  Independent  of  pipe  pressure.  The 
spreading  rates  for  tiie  fountain  .and  radl.'il  Jet  tested  are  quite  different. 
Ho'vi'ver,  it  lias  lieen  shown  In  Reference  8 that  the  spreading  rate  for  a radial 
Jet  depends  on  the  spacing  betvteen  the  radially  opposiid  iioxzle  exit  planes,  Ax^. 
Figure  Ay  sliows  the  spreading  rate  for  i.adlal  .lets,  as  a function  of  (d^yAx). 

The  curve  Is  based  on  the  data  of  Refovence  8 Wlien  Ax  Is  large,  tlu'  opposed 
Jets  develop  as  free- jets  before'  Lliey  Impinge  uptm  eaih  other,  Tlie  spreading 
rate  of  the  resultant  radial  Jet  Is  then  very  high.  As  Ax  decreases  the  spreading 
rate  decrear-'i'S , until  a point  Is  reaehed  at  which  ttu’  space  between  the  let  exit 
planes  tu'gins  to  act  as  an  annul.'ir  no/.xli'.  With  further  decreases  in  Ax,  tlu' 
spi'ead  1 11)’,  rate  he'-.i,ies  Independent  of  Ax.  In  Reference  H,  Wltr.e  refers  to  this 
ty|)e  of  radial  )i't  as  .1  ''.•on;- 1 r;i  1 nod"  radial  |oi.  Tho  indicated  fountain  data 
lie  falrlv  close  to  the  eons  t r.i  liuol  radial  Jel  results  troiii  Ki'ferelice  H. 

Figure  'jO  shows  the  decay  of  tlic  I’oniia  1 1 r.od  pi'ak  velocity  ‘ts  a 

function  of  distance  from  the  fouiilalii  virtual  crlgln,  or  dlslaiic.L'  from  the  radial 
Jet  axis.  Tin.'  data  corf.-late  will  tor  difforoiU  valnos  of  I ho  )ilpi  pri  ro. 

Till'  fountain  ciirvos  mergo  Into  tho  radial  lot  .locav  curvos,  hut  tho  value  of 

di'lu'lld.s  ell  tho  lot  s|i.;c!ii)',  S./d^  . Figure  Mi  sogi’.est  s that  fo  |-  different 

let  upacings  S,  the  fountain  iiiav  have  lo  he  'ilmulatid  hv  conn  1 ra  i iied  radial  lels 
of  dir  f e n n 1 Initial  rad  ins  R | , , 
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FIGURE  48 

RADIAL  JET  AND  FOUNTAIN  SPREADI.JQ  RATES 
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FIGURE  49 

SPREADING  RATE  FOR  FOUNTAIN,  MCAIR  R.\DIAL  JET  DATA,  AND 
RADIAL  JETS  TESTED  BY  WITZE 
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FIGURE  50 

PEAK  VELOCITY  IN  A FOUNTAIN,  AND  A SINGLE  RADIAL  JET 


fjumnuuvi.f! Ing  the  eonclusioriH  reached  from  an  analyHlH  of  the  fountain 
and  radl'ii  Jet  datni 

1.  Tlie  flow  i.n  a symmetric  two-Jet  fn  uln  may  be  rcs'itdcd  iik  a radial, 
flow  from  a virtual  origin  located  below  the  ground  plane, 

2.  At  constant  radius  From  the  virtual  origin,  tlie  fountain  velocity 
profiles  are  ,ipproxlnjtoly  equal  at  angles  of  0“  and  20“  from  the 
vartical,  but  higher  at  A0°  off  the  vertical.  This  resiult,  however, 
in  tnconclualvo  because  of  poBslhie  .-isymmetries  in  the  jot  noazlns. 

3.  Fountain  velocity  distributions  are  sol f-r Imll ar , and  the  shape  is  in 
good  agreement  with  that  for  a radial  jet, 

4.  Fountain  spreading  rate  Is  Independent  of  tin,'  jot  simcing  S,  and  is  in 
good  agreement  with  that  measured  In  Keforem  o H for  ' oouBtra  Inod"  radial 
jets . 

5.  Fountaoa  velocity  decay  l.s  In  reanon.iblo  agroomont  with  that  for  a radial 
jet,  but  the  d'ameter  of  the  equivalent  annular  nozzle  used  to  generate 

a indinl  jet  may  depend  on  the  spacing  F. 

In  spite  of  some  remaining  uncertainties,  tho’  analogy  between  a symmetric 
two-jot  fouHLaln  and  an  axlsymmetrie  constrained  radial  jot  has  been  subs lant  1 ated 
by  the  data,  The  analogy  permits  a definition  of  tho  fountain  upwash  flowflold 
and,  as  deso.rlhod  In  Section  fi,  lemis  to  a molhod  for  prePi 'ling  fountain  Impinge- 
mont  forces  on  the  underside  of  an  airframe.  .Some  ginieral  conclusions  nppardlng 
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t:helMc)Win  thtce  and  Conr-Jc’i  syrarantrlc  fountirlnn  may  also  ha  drawn,  and  thoy  are 
diacusHed  in  the  next  section, 

.'5,2  Extension  of  Two- Jet  I’mmtujn  l^eaujja_  -'.‘.'yL  l''oa‘'"de ^ ..ia>n f 

An  Lilual  1/at  Ion  of  the  flow  t'lrnlni'  process  when  two  erpia  1 axlsymmotrlc 
wall-|L'i;,s  mi'Ct.  to  lorni  a lountaln  has  led  to  an  analony  between  the  tountaln  and 
a radial  flow.  This  etiulva  leiu-e  has  been  substantiated  by  exper  1 mi'O  t . The  same 
Ideal  l/ed  I liiw  Inntlns  proei'ss  will  nowi  be  used  to  deduc.i'  (ountain  con  f I I'jira  1 1 otis 
for  itrranBemeats  ol.  three  and  loar  perpiaull  cal  ar  je' s ol  etpial  thrust..  lieeaiau' 
of  Its  greater  simplicity,  a square  .ar rangemcait  of  four  Jets  will  be  discuasud 
rirat. 

Elfiure  ,5i  La  an  ideailzod  aketch  of  the  (;rmmd  Impingement  process  for 
four  equal  vert  Leal  Jets  at  equal  holRlits  11  from  the  ground  plane.  Alonn  uucb 
of  the  llnea  OA,  OH,  OC,  and  01)  In  Figure  .51,  the  wall-jets  are  turned  such  that 
the  angle  0 defined  In  Figure  hi  la  preaerved.  Therefore,  the  fountain  segments 
generated  above  each  of  these  separation  lines  are  radial  flows.  Unlike  the.  two- 
jet  case,  however,  these  r.idl.il  flows  along  vertical  planes  will  Implnr.e  upon  ench 
other  In  tlie  neighborhood  of  the  vertical  line  througlt  the  point  0.  The  radial 
f.I.ows  are  turned  once  again  In  tills  region,  and  a central  fountain  is  formed 
whose  flow  direction  is  expected  to  be  primarily  vertical,  Examination  of 
Figure  31  shows  tliat  only  n sector  of  the  wall -jots  generated  by  jets  1 and  2 
will  contribute  to  the  fountain  segment  above  the  line  015.  The  rest  will  I'on- 
trlbute  to  the  segments  formed  above  the  lines  OC  and  01),  llie  .analogy  between 
tlie  four-jot  fountain  and  radial  jets  may  then  be  formulated  ns  Indicated  In 
Figure  32.  The  ground  plane  Is  "cut"  .along  the  diagonals  OP  | , OP^,  OP.j,  .and  OP^^, 

and  then  folded  .along  the  l ines  OA,  015,  OC,  and  01),  such  tlial  four  radial  jets  are 

formed  with  origins  at  a distance  (S/2)  below  the  ground  plane.  Note  that  for 
each  radial  jet,  only  the  sector 

a < < 5r 

>7  ' " h" 

Is  relevant.  Indeed,  only  the  sector  (^|-  2 0 " -”)  will  contribute  to  .a  vortical 
fountain  force.  The  four  nidlal  jets  cn.alehco  at  point  0,  forming  the  CLUitral 
fountain.  A flow  pattern  such  •a',  th.at  Indicated  In  Figuri'  32  may  he  geiier.aliMl 
by  four  p.alra  of  opposing  axlsynmoirlc  ji-tn,  jarovUlcd  that  the  radial  jet.s  a.re 
pn.'vented  fiami  mixing  until  tho  point  0 Is  renchod. 

Figatu’  3'5  shows  the  Impingement  points  for  a rectangular  arrangoiient.  of 
four  perpend  1 cu  I ,a  1"  jets  of  equ.al  Ihrust.  'I'he  eijulvalent  yi  r laangenuai  t of  l^arllal 
jets  would  In  tills  case  he  developed  by  oneo  .again  "eat  ting"  the  grouaJ  plane 
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RECTANGULAR  ARRANGEMENT  OF  FOUR  EQUAL,  PERPENDICULAR  JETS 


alonp,  dlagonaiH  0I’|,  0I’2t  01.’.^  ami  OP^t  nnil  then  foliltng  tlu*  plane  along  symmetry 
linos  OA,  OB,  OC,  and  Oil.  The  virtual  origins  for  the  OA  and  OC  fountain 
segments  now  lie  tit  a dlstatieu  (L/2)  below  the  ground  plane,  and  the  origins 
for  fountain  Hogments  OB  and  01)  are  located  at  (S/2)  below  the  ground  plane. 

The  releviint  radial  Jet  sectors  will  now  bes 


Fountain  segments  OA  and  OC:  Oj^  - 0 - ,y 
I'onntiiln  segments  OB  nml  01):  0^^  - o - ~ 
where  0.j  2 "•  tun  ^ (S/L) 

'^lA  “ '2  " ”l2 

Figure  54  depicts  the  Idetilltied  flow  pattern  for  three  etiual  perpendicular 
Jets  in  an  Isosceles  triangle  arrangement.  Throe  fountain  Hegraenls  are  juMienifed 
In  this  eahef  OA,  OB,  and  00.  Because  of  Hvmmetry,  vlie  fountain  MegtiteiUs 
Lie  along  the  perpendicular  bisectors  of  the  triangle  sides.  The  virtual  origin 
for  the  fountain  segment  OB  in  located  at  a distance  (R/2)  below  the  ground  plane, 
and  tin’  virtu.il  origins  for  segments  OA  and  OC  are  located  at  (I./2)  below  tlie 
ground  plane.  Tlie  equlvaUait  radial  Jets  are  now  made  up  by  Lhtee  palim  of 
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ISOSCELES  TRIANGULAR  ARRANGEMENT  UP  THREE  EQUAL, 
PERPENDICULAR  JETS 

QPtt'OSM'U 

opimBod  Jets,  wltli  each  pair  ci'iitered  about  each  of  tiiu  virtual  orisitiH.  Xlu? 
relevant  radial  Jet  sectors  are! 

Fountain  segments  OA  and  OCi  - 0 ■* 

•C  K “Jf 

Fountain  segment  OH:  - G - 

where  from  the  geometry  Jn  Figure  54: 


0^2 


O23  - cos  (2-j|-)  - G32  • 


As  in  the  fi.iur-Jet  arrangements,  a central  fountain  Is  formed  above  point  0, 
where  the  tliroe  fountain  segments  impinge  on  each  other. 

lining  tlie  equivalence  between  fountain  segments  and  radial  Jets,  it  is 
possible  to  compute  the  vertical  momentum  flux  for  various  Jet  arrangemmits . 
This  will  be.  done  In  the  next  Muli.H'n-t  Ion . 

'5 . J Normal  Monienttim  Flux  In  I'nuntaiiis 

Th(‘  fountain  impingement  forie  on  tlie  underside  of  an  airframe  Is  related 


•.  •,'**.  V • vn;  -Pv 
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to  the  total  vertical  momentum  flux  In  the  fcnintaln  flow,  and  calculation  of 
tills  momentum  flux  provides  a measure  <if  the  maximum  fountain  ImplnRcmcnt  force 
nvullahle  for  n particular  .let  .irranRoinont . The  fountain  velocity  Is  not 
vertical,  so  It  Is  Important  to  take  Into  account  Liu'  momentum  flux  components . 

Tlie  analogy  between  fovmtaln  ami  sliigle  radial  JetH,  or  como  I nat  loiia  of  radial 
JetH,  pi'.riiiitK  calculathjn  of  tlie.se  comininenta , 

A.s  dlHCU.M.sed  In  prevloti.a  Hi'ctloiiH,  .symmetr  li’.a  1 Lwo-jiM  roantaln,s  havt 
lu'e  a repl.aced  hv  ef|ii  1 v.a  1 eii  t i.adlal  Ji’l.s.  I''nuntnln.a  for  con  f I p.urat  1 omk  wlLli 
three  and  four  perpendicular  Jets  of  et|ual  thrust  liave  been  divided  into  sogment.s 
and  each  segment  related  to  the  flow  from  an  appropriate  sector  of  a radial  ,)et. 
Some  of  tlie  fount, '.in  segments  contain  portions  of  a central  Fountain,  wli  Ich  has 
been  taken  to  tse  oquivaient  to  the  flow  generated  by  the  Impingement  of  the 
various  radial  jet.  sectors,  fliven  a group  of  two  to  four  vertical  jets,  the 
calculation  of  the  rosultant  fountain  vertical  momentum  flux  requires  the 
following  information! 

1.  A relation  between  the  momentum  flux  of  the  vertical  jets  and  wall” 
jets  wliicli  penerati'  the  fountain,  and  file  radial  Jets  wlilcli  have 
replaced  the  rountnln. 

2.  Calculation  of  the  vertical  momentum  flux  for  tlie  various  radial 

Jet  sectors  which  make  up  the  fountain,  and  for  any  central  fountain 
wlilch  may  have  been  created  by  the  Impingement  of  these  sect-.ors, 
item  [ is  discuHHed  below,  and  Item  2 is  tiu'  I'oriowlng  Hubsecti.on , 

It  has  been  shown  that  when  a circular  jet  impinges  at  9U®  to  a ground 
plane,  the  total  momentum  flux  Is  conserved,  and  at  the  bmiiulary  of  the 
"impingement  zone"  In  Figure  12,  the  total  momentum  flux  per  unit  azihiuthal 
angle  in  the  wall-jet  Is  equal  to  the  total  momentum  flux  per  unit  angle  in 
the.  frec-Jut,  Furthermore,  it  has  been  .shown  In  Section  2i').'3  that  as  tlie  wall- 
jets  develop,  the  total  momentum  flux  per  unit  angle  remains  unchanged.  Indicating 
that  the  vLscou.s  shear  stress  on  the  grouml  plane  docs  not  significantly  affect 
the  wall-jet  development.  It  will  now  he  assumed  that  when  two  ax Isymmet ri c wall- 
jets  Impinge  upon  each  other  and  form  a fountain,  the  tot.il  momentum  flux  per 
unit  angle  Is  also  conserved  In  the  turning  process,  Since  the  fountain  combines 
the  momentum  from  the  wall  jets  on  each  side  of  It,  the  total  momentum  Mux  per 

unit  angle  In  the.  foimtaltl  will  be  equal  to  twice  tlie  momentum  flux  per  unit 

angle.  In  a single  wall  jet.  Since  the  lountaln  and  eciulv.alent  radial  let  have 

the  same  momentum  flux,  It  follows  (see  Figure  A2)  that  It 


ilonote  tlio  niomuiiUum  fluxnH  In  the  vevtleal  JctH,  tlicn 


where  Mj|,  ilennti'H  the  l.otnl  nutiiienHiin  flux  In  the  rmllnl  Ji'l  whieli  ri'pl;iees  the 
I'eun ,1 1 11 . The  flow  In  n rotintiiln  nen.iient  1h  in  the  riiillnl  (ilreetlon  nlioiit  the 
fiej’.meiU  ' s vlnn.il  orlpln,  no  the  a v.i  I I .-ih  I e verileal  iiiiiimai  I.  nm  llnx  Is  nn  I y a 
ftiU'tlon  of  the  totii.l  HiiRinont'H  momentuin  flux.  The  mnunltncle  of  this  fraction 
will,  be  computed  in  the  Following  HubHectlon. 

5.3.1  IJeiil  Noriiml  Momentum  Flux  in  a Pmintnln  Seainent  - In  rignres 
51,  51),  and  54,  titc  waU-;let’H  flowfleld  has  boon  subcltvlded  into  an  inner 
region  .and  an  outer  region.  The  boundary  between  theMe  1h  the  polygon  constructed 
by  eonnee.tlng  the  Jet  Impingement  points.  Since  the  Jets  are  all  perpendicular 
iind  of  equal  thrust,  e.uch  fountain  segment  will  form  along  the  perpendicular 
bisector  of  each  side  of  the  polygon.  Thus,  one  part  c.f  the  fountain  segment  lies  at 
least  partially  In  tlie  inner  region,  /ind  is  denoti'd  the  "Inner  region  sector." 
llie  rest  of  the  segment  lies  in  tlic  outer  region  and  is  correspcaidingly  called  the 
"outer  region  sector."  Thu  liver  region  sectors  from  dific’nnt  Jet  pairs  will 
intersect  to  form  one  or  more  central  fountains. 

i'lie  Jet  ai'i'angenieiits  will  nov./  be  restrli’tod  to  cases  wlicre  all  Inner 
region  sectors  Intersect  wltliln  the  timer  region,  llie  geometries  in  Figures  51, 

5),  and  54  all  satisfy  this  eondltiun.*  Figure  55  shows  a rase  whore  tills  is  nci_t  so. 

At  central  fountains,  the  radial  Jets  simulating  the  varliuis  segments  Impinge 
on  each  other.  I'lie  veitie.al  inotnentum  flux  across  a plane  parallel  to  tliu  ground, 
located  a distance  z above  the  ground,  will  partly  come  from  undisturbed  portions 
of  the  radial  Jets,  and  partly  from  the  central  fountain  or  rountalns.  However, 
for  central  fouiualns  wlifch  lie  within  the  inner  region,  exami  nut  ion  of  Figures 
5 1 and  511  will  reveal  that  there  ,1s  no  vertical  imMiiuiUuin  flux  l-|iroiii’h  flu'  lateral 
area  of  the  Inner  region.  Conset|uenlly , the  total  vortical  momentum  flux  acrosH 
planes  pa. Ml  lei  to  the  ground  plane  Is  liulependent  of  r. , ,ind  can  tlierefoi-L'  In' 
computed  at  x 0.  That  Is,  where  the  grmmd  plane  was  located.  At  the  plane 
X » U,  It  is  as.sumeil  that  the  contra!  fountain  li;is  not  yet  foi-med,  and  conHi'tinont ly 
the  vertical  moiiK'iitum  flux  In  tlu'  region  may  lie  ca  I cii  I aied  from  llie  momentum  dlsirl- 
hut  loll  for  Liu.'  Indl'ildiial  r.ulial  Jet.  si'ctors. 

Tills  Is  not.  a serious  restriction,  sinee  most  praci  leal  designs  would  seek  to  maxi- 
mlxo  the  fountain's  Impingement  forei'  by  local Ing  l bo  ccnlral  ioimlala  within  the 
Inner  region,  when'  It  will  u.sn.-illv  Impact  on  .'i  p, renter  aircraft  planforrii  area. 
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FIQURE65 

ISOSCELES  TRIANGULAR  ARRANGEMENT  OF  THREE  EQUAL  JETS  WITH 
FOUNTAIN  SEGMENTS  INTERSECTING  OUTSIDE  THE  INNER  REGION 

F'Ugaru  56  hHowh  tlu;  wall  Jots  and  equivalunt  rmlial  |ut  sectora  Tor  otu! 
Hemmint  u(!  a multip lo-jot  fountalii4  Tho  fountain  vortical  iTiomontum  In  tln>  Inner 
nasion  Hector  in  c.ilculnted  by  computing  the  vi-rtical  momentum  flux  In  the  radial 
Jot  Huctnr  Hubtundod  by  thu  line  OA  In  Kiguro  56b,  The  fountain  vortical  momontnm 
in  the  outer  region  Is  calculated  by  computing  thu  vertical  momentum  flux  In  the 
radial  jot  Hector  Hubtondod  by  n finite  value  of  y,  and  then  letting  y approach 
Infinity,  Thu  charactet'lHtics  of  constrained  radial  jets  doHcribed  in  Reference  5fi 
will  be  used  to  calculate  the  vertical  momentum  in  the  fountain  sector,  Fur  a 
strip  of  width  dy , the  momentum  flux  In  tlie  x-d 1 recL I on  1h  given  by 
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dx  dy 


(6'1) 


where  w Is  the  radial  jet  velocity  In  the  x-d  1 rei' 1 1 on 


W " Vj^  COH  II 

where  V|^  1h  the  radial  )ct  velocity  In  (lie  R dlrei'tlon,  l'',(|uat  I on  (bU)  niiiv  I hen  he 
wr I tten  aw 
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FIGURE  be 

WALL  JETS  AND  EQUIVALENT  RADIAL  JET  FOR  A FOUNTAIN  SEGMENT 


dM^  ■ coh“  0 I dx  dy. 


From  total,  momoiUum  consorvat ton  In  the  radial  let 


M - 2ffU  I oV,^  dx,  K '>  K, 


where,  ax  before,  Nj|^  in  the  total  momentum  flux  throuBh  the  radial  .iet'a 
annular  nozKle  (of  radiuH  Rj). 

ConHequently , Kquation  (70)  may  be  written  aa! 

• Mil,  2,, 

V.  2 71  R 


M I m;t! 


y “ tan  fi. 


R^j/R  ••  eoH  I), 


till’  total  vertte;il  momentum  flux  becomoH  a Imply 


cnM  0 tlf-> 


. M, 


or 


^’iK 


I sin  (ij.  -■  nln  ii| 


(12) 


wIru'l'  aiul  i)|^'  dunoLi.'  tliu  lnlLl;il.  a\ul  I’liiiil  limits  ol'  I n Lu  grn  L liMi , ri'S|u’cl' Ivi' ly 
Keferrlus  to  Flsuro  56, in  the  limit 
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So,  for  the  totnl  sector,  the  vertical  momentum  flux  la 
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(7.'3) 


wlu'.ro  ,,  Is  positive 


The  simple  result  In  Equation  (73)  may  he  used  to  enl'Milate  the  ideal  vertical 
momentum  flux  for  the  eombln.itlon  o''  radial  Jet  Mectoirs  which  make  up  mu  1 1 Ip  1 e- J et 
fuuu tains. 

I) ,3,2  tde.il  flomeiitum  Kluxes  for  IVo,  Three , ajytl Euuntalns,  and 
Optimum  Ar rantjements  of  Three  and  I'our-.let  Epun tains  - For  a twi>-Jet  fountain,  the 
an«lc  (i|2  In  Figure  56  Is  it/2,  and  Equation  (73),  then  yields 


II 


and  since,  as  sliown  before, 
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wliei'e  Mj  denotes  the  momentum  flux  for  eneh  one  of  tlie  two  Jets,  It  follown  th.il. 
the  Ideal  fountain  momentum  flux  divided  by  the  total  |et  tlirurt  will  hel 
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An  "optimum"  thiri’e  or  four-let  fountain  ananp,eiiieiU  Is  conaideri'd  to  In 
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mu'  Tur  which  tho  total  fountain  vortical,  momentum  fa  a maximum.  The  condltionH 
unclor  which  this  occurs  will  be  derived. 

Considering  first  tho  Ihoscu.Iuh  trlanKulnr  nrranRomnnt  In  riguro  "jA  and 
Hummlnn  the  radial  Jet  sector  contr  Hint  Ions  uslnH  Kt|iiatlnn  (7'l)i 
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now 


and 
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where  M Is  tho  momentum  flux  for  u mIukIo  vortical  Jot;  and  also, 
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Kciuat  Ion  (^.l)  then  hccomea 
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and  dividing,  bv  the  total  vert  leal  )et  im'imniliim! 
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l'’nitti  the  expresalt  1H  for  Ci  ,intl  ‘Vi  Ine  Inch'd  In  Flnufe  J<'*i  thuae  anulea 
depeiul  on  the  ratio  (S/I.).  Ill  florent  lat  Inn  Kt|iiatlon  (7b)  with  reapeet  to  (S/l.) 
.ind  .ael  tinu  tho  rcanl  t ot|ual  to  zero  Jeacla  tot 


I 2 I'OH  1 1 1 ,j  - I’oa  ,, ) 


2'i' 


d 

d'(H7‘l.‘) 


(Bit.-  (-^Y) ' '*2j  d‘(s/i:) 


I 


I). 


Kxpanclinp,  the  derlvatlvea  I e/uln  to  the  i.Hult-  that 


showing  that  tho  liiaxlmum  Ideal  fountain  momentum  flux  la  achieved  when  the  jets 
arc  arranged  in  an  equilateral  pattern. 

Considering  the  rectangulai*  Jet  arrangement  slmwn  In  Figure  5'3  and  summing 
the  radial  Jet  sector  eont rl lint  1 ons  using  Kqufitlon  (71),  the  total  vertical 
momentum  flux  in  the  fountain  becomes 


Mjr 

^ ®12  ^ 

where  geometrical  aymmetry  has  been  taken  Into  account. 

Using  the  relation  between  radial  and  vertical  Jet  momentum  fluxes, 
the  above  may  be  written  as: 

tain  G^2  + 0l4  + 21. 

Dividing  again  by  the  total  momentum  flux  of  the  vortical  Jets 
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By  finding  the  maximum  value  of  this  function,  it  is  possible  to  show 
that  the  maximum  ideal  fountain  vertical  momentum  flux  is  achieved  when  the  Jets 
are  arranged  In  a square  pattern. 

Figure  57  shows  the  ratio  of  Ideal  fountain  vertical  momentum  flux  to  total 
jet  momentum  flux,  plotted  against  S/L,  for  two,  three,  and  four-jet  fountains. 

The  figure  indicafes  that  the  variation  with  (S/L)  is  very  weak,  so  that 
configurations  whicl^  differ  from  the  optimum  equilaf  f^rnl  triangular,  or 
quadrilateral  square  patterns  will  experience  a very  sms’l  reductlo.i  in  fountain 
vertical  momentum  flux. 

5.4  Asymmetric  Twu-Jet  Fountains 

The  dracrlprlon  ( I two-jet  fountain  flowflelds  has  so  far  been  confined  to 
cases  in  which  the  two  jets  that  generate  the  fountains  are  pnrpendicular  and  have 
equal  thrust.  In  Figure  58,  this  simplest  jet  arrangement  is  labeled  "configuration 
A."  A more  complicated  situation  arises  for  configuration  B,  in  which  the  jets 
arc  vectored  as  mirror  images,  wbfle  their  thrusts  are  stil  i kept  equal.  Tlie  flow- 
field  Is  still  symmetric  about  the  plane  3,  and  consequently  the  fountain  remains 
planar  and  the  fountain  velocity  profile  is  symmetric  about  3-  The  discussion  of 
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FIGURE  57 

NORMALIZED  IDEAL  FOUNTAIN  VERTICAL  MOMENTUM  FLUX 
FOR  TWO.  THREE.  AND  FOUR-JET  FOUNTAINS 
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FIGURE  58 

TVO-JET  FOUNTAIN  CONFIGURATION 
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oblique  jet  Impingement  in  Section  2.2,3  has  indicated  that  the  free  jet  impingement 
stagnation  point  on  the  ground  plane  Is  displaced  from  the  projected  jot  asda  by 
a distance  AR^.  Thus,  the  minimum  distance  between  the  jet  Impingement  poiiit  and 
the  symmetry  plane  3 is  given  by 


S'/2  = S/2  - H tan  + AR.j. . 

If  it  is  once  again  assumed  that  In  forming  the  fountain,  the  wall- jets  are  turned 
as  if  they  were  "folded"  along  the  Intersection  between  the  plane  3 and  the  ground 
plane.  It  follows  that  for  configuration  B the  fountain  flow  is  analogous  to  a 
radial  flow  emanating  from  a virtual  origin  located  at  a dlMtnnce  S'/2  below  the 
ground  plane.  Due  to  the  asymmetric  momentum  distribution  in  toe  wall- jet,  however, 
the  azimuthal  momentum  distribution  In  this  fountain  will  differ  from  that  for 
configuration  A.  Assuming  momentum  conservation  in  the  turning  process,  the 
azimuthal  momentum  distribution  In  the  fountain  may  be  calculated  from  the  known 
momentum  distribution  In  the  wall-jets. 

In  Figure  bti  configurations  C and  I)  are  no  longer  symmetric,  so  tlwit  the 
fountains  no  longer  Me  on  a plane.  I'uc  to  n lack  of  data,  it  Is  impossUile 
to  specify  the  fountain  flowfleld  with  any  certainty  In  these  cases,  but  a 
qualitative  picture  of  the  fountain  geometry  has  been  established.  For  configura- 
tion C,  the  "separation  surface"  (analoguua  to  the  plane  3)  is  postulated  tu  be 
conical  with  the  cone  apex  located  below  the  ground  plane.  Fiirurc  59  deplete 
the  conical  separation  surface  originating  from  a hyperbolic  ground  separation 
line.  The  cone  semi-apex  angle  Is  denoted  by  Og,  and  the  cono  apex  lies  on  the 

z axis,  at  a distance  z below  the  origin.  The  fountain  flow  would  then  be  conical, 

o 

with  streamlines  radiating  from  a virtual  origin  at  the  cone  apex.  A conical 
separation  surface  has  several  properties  which  make  it  plausible: 

1.  The  intersection  of  the  cone  and  ground  plane  is  a hyperbola,  and  It 
will  be  shown  that  this  is  compatible  with  the  separation  lines 
calculated  In  Section  3.2,2. 

2.  With  the  fountain  velocity  vectors  located  approximately  along 
generatrices  of  the  cone,  the  fountain  streamlines  are  straight. 

This  la  consistent  with  a constant  pressure.  Free  turbulent  flow. 

3.  As  the  jot  momenta  Mp  and  become  equal,  the  cone  semi-apex  angle  can 

approach  n/2,  and  the  planar  radial  flow  postulated  for  configuration 

A la  recovered. 

Section  3.2,2  has  described  calculations  of  the  wall-jet  separation  lines 
for  a pair  of  vertical  jets  operating  at  different  thrusta.  Results  of  the.se 
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FIGURE  69 

FOUNTAIN  GEOMETRY  FOR  PERPENDICULAR  JETS  OF  UNP.QUAL  THRUST 


cfjlculiitiona  have  been  used  to  obtain  an  expression  for  the  asymptotic  angle  of 
a hyperbolic  fit  to  the  seoaration  lines,  and  thus  to  the  conical  fountain's  semi- 
apex  angle.  The  points  in  Figure  60q  indicate  the  asymptotic  angles  measured  from 

thrj  exact  numerical  results.  The  dotted  line  in  the  same  Figure  shows  that  the 
points  are  closely  approximated  by  the  expression 

. . 1/2 

(1.56) 

tan  n ,~Ai — . (78) 

° I - (M^/Mp) 


The  location  of  the  cone  apex  may  be  obtained  from  tl>e  location  of  tlie 
separation  point  P on  the  x axis  of  Figure  59.  From  the  etagaation  line 
equations  it  may  be  shown  that  for  two  perpendicular  .)etsi 


1 

1 + (M^/MjJ 


(79) 
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FIGURE  eOa 

CONICAL  FOUNTAIN  SEMIAPEX  ANGLES  AS  A FUNCTION  OF 
JET  MOMENTUM  RATIO 


where 


Xi  - 2x^/S  , 


From  the  geometry  of  Figure  59 


z - -X.  ten  a 
0 i o 


where 

(79) 


a bar  again  denotes  normalization  by  (S/2)  , and  using  Equations 


. . 1/2 
(1.56)  (M^/Mp) 

1 + (M^/Mj,) 


(78)  nnd 


(80) 


A pint  of  Equation  (80)  Is  shown  by  the  dotted  line  In  Figure  60b.  Equation 
(80)  docs  not  approach  the  correct  limit 
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Jit  Mominturn  Ratio  , 
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FIGURE  eOb 

CONICAL  FOUNTAIN  APEX  LOCATION  AS  A FUNCTION  OP 
JET  MOMENTUM  RATIO 
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Adjusting  constants  In  Equations  (78)  and  (79)  so  that  this  condition  is 
satisfied  yields  the  relations 


tan  « 


2 (M^/Mp) 


1/2 


o 


1 - (M^/Mp) 


. . 1/2 

2 

1 + (M,/MJ 


Equations  (81)  and  (82)  are  plotted  as  solid  lines  in  Elfiures  60n  and  60b 

« • 

respectively.  For  (M^/M^)  > 0,5,  Equation  (81)  still  provides  a good  approxi- 
mntion  to  the  computer  program  results.  Equations  I'or  the  complete  stagnation 
line  may  now  he  written.  Choosing  the  numarlr.nl  constants  for  the  correct 
limiting  beliavior  as  (M^/M^)  "*■  1,  and  For  a good  fit  to  the  calculated  results, 
the  equations  are: 


. . 1/2 

2 

1 - (M^/Mp) 

' 1 

y =■  - - . . - 

1 - (M^/Mp) 

X 

1 4-  (M^/Mp) 

1 

ai:d 


for  0.5  1 (VMp)  < 1 


_ (1.56)  (M^/Mj.)^^^  1 

[ 

|y2  _ 

'i  - tvs' 

2 

y • * ‘ 

1 - (M^/My)  1 

\ ^ 

1 + (M^/Mp) 

f(J4) 


for 


(lyMj.)  1 0.5. 


These  equations  are  compared  to  the  computer  program  calculatloni  in 
Figure  61. 

The  above  discussion  has  provided  a plausible  qualitative  description  of  the 
fountain  geometry  for  the  case  of  two  Jets  of  different  thrust  impinging  per- 
pendicularly to  a ground  plane.  In  order  to  compute  the  force  generated  by  such 
an  asymmetric  fountain  on  an  airframe,  data  on  fountain  velocity  profiles  and 
spreading  rates  are  required.  This  informat J.on  is  not  currently  available. 

Indeed,  the  whole  description  of  this  asymmetric  fountain  geometry  is  speculative 
at  present  I and  requires  experimental  verification.  When  the  Jete  have  different 
thrust  and  are  vectored  at  different  angles  (configuration  U In  Figure  58)  the 
picture  is  complicated  still  further.  It  may  still  be  possible  to  poatulahe 
a conical  fountain  flow  for  which  conditions  vary  for  each  meridional  angle. 

The  overall  fountain  geomatry  might  still  be  obtained  by  adjusting  cone  parameters 
to  fit  calculated  or  measured  valijes  of  the  fountain  sci  aratlon  Line. 
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C^lcuIttUd  (from  Section  3.3.2} 
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FIGURE  01 

APPROXIMATION  TO  CALCULATED  STAGNATION  LINES  BY  HYPERBOLAS 
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6.  FOUNTAIN  IMPINGEMENT  ON  A FUSELAGE-WING 
A description  of  the  fountain  flowflelds  generated  by  two  to  four  circtilar, 
vertical  Jets  of  equal  thrust  impinging  on  a ground  plane  has  been  given  in 
Section  5,  The  total,  ideal  vertical  momentum  flux  generated  in  such  fountains 
has  also  been  calculated.  This  section  describes  how  this  vertical  momentum 
flux  is  converted  into  a positive  normal  force  when  a fountain  strikes  the  air- 
frame undersurface.  It  has  also  been  found  that  in  addition  to  this  positive 
impingement  force,  the  fountain  induces  a negative  suckdown  force.  This  force  is 
to  be  added  to  the  suckdown  forces  generated  by  the  free-Jeta  and  wall-jets,  and  it 
it  related  to  the  fountain  entrainment.  An  approximate  method  to  account  for  the 
fountain  suckdown  force  has  been  derived,  and  is  also  described  in  this  Section. 

In  a parallel  MCAIR  IRAD  effort,  (Reference  23),  an  apparatus  to  simulate  two 
and  four-jet  fountain  configurations  was  designed  and  fabricated.  The  forces  pro- 
duced on  a representative  aircraft  configuration  by  these  simulated  fountains  vjere 
measured.  The  experimental  forces  have  been  compared  to  the  sum  of  predicted  foun- 
tain impingement  and  suckdown  forces,  and  tlie  results  are  described  in  this  Section. 

6. 1 Two- Jet  Fountain  Impingement  and  Suckdown 

In  Section  5,  the  flow  in  the  fountain  produced  by  two  circular,  vertical  jets 
of  equal  thrust  has  been  shown  to  be  approximately  equivalent  to  the  f!low  in  a con- 
strained radial  jet.  This  analogy  will  be  used  to  calculate  the  pressures  induced 
by  fountain  impingement,  and  the  terms  "two-jet  fountain"  and  "radial  jet"  will  be 
used  interchangeably,  Based  on  available  data  on  the  pressures  generated  by  the 
impingement  of  two  axlsymmetric  wall  Jets,  a simple,  relationship  between  local  jet 
dynamic  pressure  and  the  pressure  generated  on  the  surface  upon  which  the  Jet 
impacts  has  been  obtained.  This  relationship  is  used  to  calculate  the  force  generated 
by  the  radial  jet  impingement  on  an  airframe  undersurface.  ITie  calculation  of  two- 
jet  fountain  suckdown  forces  i,s  also  discussed  in  this  subsection.  The  simpler  case 
of  suckdown  produced  by  a radial  jet  impinging  on  a rectangular  wing  Is  discussed 
first,  and  the  more  complex  case  of  a trapeaoidal  wing  is  discussed  later, 

6.1.1  Fountain  Implngamaut  Force  - The  constrained  radial  jet  used  to  represent 
the  two-jet  fountain  is  an  axlsynmetric,  turbulent  flow  which  spreads  out  from  an 
annular  nozzle  and  whose  width  grows  linearly  with  distance  from  the  nozzle  radius 
(See  Reference  8 and  Figure  48).  Referring  to  Figures  42  and  43,  the  analogy  between 
radial  jet  and  fountain  implies  that  the  radial  jet  development  on  each  side  of  the 
symmetry  plana  is  approximately  equivalent  to  that  of  the  corresponding  wall  Jets 
which  merge  to  form  the  fountain.  It  has  been  postulated  (and  experimentally  confirmed) 
that  when  the  wall-jets  turn  to  form  the  fountain,  the  velocity  components  normal 
to  the  plane  of  symmetry  are  rotated  90®,  and  the  tangential  components  are  unchanged. 
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(See  the  enlarged  sketch  in  Figure  A2) . Il  in  poHtulatod  tliia  the  Hume 
turning  process  occurs  when  the  radial  jet  impinges  on  the  underside  of  an  air- 
craft planform.  Figure  62  depicts  the  impingement  of  a rndiul  Jet  on  u rectangular 
plate,  for  a completely  aymmotricul  situation.  The  sectors  of  the  radial  Jet  which 
do  not  Impinge  on  tlie  plate  ate  assumed  to  remai.n  und  Isturhed , Tlui  deflected  sectnr 
Is  first  split  along  the  line  II,,  and  then  deflo(M:uil  such  that  the  local  any,Je  0 
as  defined  in  I'l.gure.  62  is  preserved  In  tlu'  turning  process.  The  direction  of  the  I 1 ow 
along  the  plate  undersurface  would  then  he  as  iiulioated  In  the  Figure, 


IDEALIZED  RADIAL  JET  SPLITTING  AND  TURNING  UPON 
IMPINGEMENT  ON  A RECTANGULAR  PLATE 


6. 1.1.1  Pressures  Generated  by  Axisymmetrir,  Wall-deta  Impxnglna  on  Kach  Other  - 
It  has  been  assumed  above  that  the  turning  process  when  two  wail-jets  meet  at  a plane, 
of  symmetry  to  form  a fountain  is  the  same  as  that  when  each  half  of  a radiaJ  Jet 
encounters  the  underside  of  an  aircraft  planform.  Consequently,  pressure  data  obtained 
for  the  former  case  can  be  used  to  derive  an  expression  for  calculating  pressures  in 
the  latter  case.  For  a symmstrical,  two  jet  arrangement,  keference  24  contains  measured 
pressure  distributions  along  a line  which  corresponds  to  the  intersection  between  thu 
ground  plane  and  the  plane  of  symmetry  (y  axis  in  Figure  42).  It  has  been  possible 
to  relate  the  measured  pressures  to  the  peak  velocity  which  would  exist  at  each  point 
P in  the  wall-jet,  If  the  jet  continued  as  a free  wall-jet. 
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Thd  basic  concept  Is  that  the  pressure  distribution  caused  by  an  impinging  tur- 
bulent flow  can  be  calculated  from  the  undisturbed  dynamic  pressure  normal  to  the 
Murfacu.  ThuB»  tl.a  peak  preaaurea  generated  when  two  wall -Jets  impinge  on  each  other 
ufo  iL'lnti'il  to  the  peak  dynamic  (ircaaiiros  on  tliu  wall-jet  prof  ilea. 

Referring  once  again  to  figure  42,  the  wall-JutB  develop  along  the  ground  plane 
It  cMHent  tally  couHtant  atutic  presHure,  outnide  of  regions  In  tlie  inum'.diate  neighbor- 
hood of  tile  fmplngement  points  1 and  2.  Near  the  (y,  n)  plane  of  symmetry,  the  wall- 
lets  turn  and  form  tlie  fountatn.  In  this  Lurnl.nc,  process,  the  pressure  on  th(‘  r, round 
plane  increases,  since  the  wall-jet  flows  are  retarded.  The  pressure  reaches  its  max- 
imum vaiue  along  the  y axis,  and  this  maximum  may  be  calculated  by  relating  It  to  the 
peak  velocity  in  the  wall-jet  profile,  • Results  In  References  7 and  8 indicate 

that  for  normal  Impingement  of  a round  jet  on  a plane,  the  peak  velocity  in  the  axis- 
ymraetrlc  wall  Jet  far  from  the  Impingement  point  may  be  calculated  from! 


Je/ 


H 


\-l 


(1.97)  ^ 


(851 


wlierei 


.a 


let  exit  velocity 


il  ■ no!!zle  height  above  ground  plane 

K,  ■ nozzle  exit  r.idlus 

Je 

This  equation  holds  for  an  established  wall-jet  far  from  the  impingement  point. 

At  n point  P along  the  y axis  where  the  radial  wall-jets  interact  (aee  Figure  42), 
it  is  pustulated  tliat  tlie  pressure  is  given  byi 

0 ,,  2 2 rt 

p - p ■ r U cos  fl 
‘ m “ 2 m 


or 


P P 1 “ P V,  ‘ 
oj  J e 


Substituting  Equation  (85) 


-2 


2 ( r ^ , 


- (1*97)  t H / t H 


and  since 


:osfl  “ 


2r 


and  along  che  y axis 


r - [y^  + (S/2)^]^^^ 
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'•S' 


this  may  be  written  in  the  form! 


s'  d' 


C * n.97i  


It  moat  be  noted  tlmt  Equation  (86)  la  independent  of  the  noziile  helglit  H.  'L'hls 
conclusion  la  supported  by  the  data  or  Reference  24.  Figure  6!3  compares  Equation (86 


to  tile  data  of  Refo.rencc  24  for  four  values  of  the  nozzle  spacing  S/D  . Tl>e  agreement 


Lh  very  good,  except  for  tin;  case 


tn  this  instance,  the  nozzles  are  probably  so  close  together  that  the  wall-jot  has 
not  become  fully  established  prior  to  the  fountain  formation. 
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» E^uBtion  (86) 
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FIGURE  63 

PEAK  PRESSURE  DISTRIBUTIONS  BENEATH  FOUNTAIN 
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iha  basic  assumption  that  the  peak  pressure  generated  by  the  impingement  of  two 
radial  wall  jets  is  equal  to  the  normal  component  of  the  maximum  dynamic  pressure  in 
the  wall  Jet  has  therefore  been  aubstantioted  by  the  data. 

6 . 1 . 1 . 2 PreHsurea  Generated  by  Fountain  Implnaemcint  on  an  Airframe  Under- 
surface - rigure  6A  liluatratea  a cwo-Jet  fountain,  represented  by  a radial  Jet 
with  origin  at  a diatanue  below  the  indicated  fuselage  reference  point,  and 
with  initial  radius  Kp  It  is  assumed  that  ambient  pressure,  acts  on  the  planform 
upper  surface,  it  is  further  assumed  that  those  portions  of  the  annular  jet  which 
are  not  blocked  by  the  planform  are  unaltered,  and  do  not  in  any  way  affect  the 
jet  impingement  force  on  the  aircraft.  The  net  force  on  a strip  of  width  dy  will 
then  be: 

X (y) 


dF 


■/' 


( / cp  - p„>  ‘^y 

Kj,(y) 


where  Xj^  and  denote  the  upper  and  lower  limits  of  integration.  Replacing  the 

aircraft  by  its  planform,  and  equating  the  pressure  difference  to  the  .armal  component 
of  the  jet  dynamic  pressure: 

X.  (y) 


dF 


^ [ 
2 ^ 


/ <v 


^ cos^  6)  dx]  dy 


(87) 


x.j(y) 


Since  the  radial  jet  has  a finite  width,  the  radial  jet  momentum  is  distributed  in 
a neighborhood  of  the  plane  x “ x^.  The  crosshatched  area  In  Figure  64  represents 
the "footprint "of  this  region  on  the  planform  located  at  z ■ 0.  It  is  assumed  that 
(for  a configuration  without  stores  or  flaps),  the  pressure  decays  to  ambient  values 
outside  this  footprint  region.  Consequently,  the  integration  limits  x,^  and  Xj^  are 
determined  by  the  shape  and  location  of  the  fountain  footprint,  since  the  only  non- 
zero force  contribution  arises  from  the  area  conmon  to  the  strip  and  footprint, 
IndLc.nted  by  the  double  crosshatching  in  Figure  64, 

The  data  of  Reference  8 are  used  to  analytically  define  the  radial  jet  velocity 
distribution,  which  then  leads  to  an  Integral  for  the  total  fountain  impingement 
force,  Letting  Ax  denote  the  total  jet  width,  the  data  of  Reference  8 have  shown 
that  for  a constrained  radial  Jet: 


Ax  - (R  - Rj)  , 


(88) 


and  from  Figure  64: 


R - , 


97 


I 


I Boundary 
l-mj  j»dy 
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Boundary 
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FIGURE  64 

CALCULATION  OF  FOUNTAIN  IMPINGEMENT  FORCE 


than  substituting  In  Equation  (88) t 


, 21^/2  R 1 

‘'t  r ^ -1^ 


For  P. /R  <<  1 the  footprint  boundary  la  given  by  two  hyperbolas  which  represent 
j ^ 

the  intersection  between  the  conically  shaped  annular  jet  and  the  plane  z • 0,  Tlie 
constant  depends  on  the  location  of  the  annular  jet  boundary,  which  is  presently 
defined  as  that  value  of  x for  which  the  velocity  has  decayed  to  a specified  fraction 
of  the  peak  value  for  the  profile  at  that  point.  The  data  of  Reference  8 indicate 
that  for  the  half-volocity  boundary 

K,  - 0.212 

t 

and  then, 


Axj  - (0.212)  (R  - Rj)  (90) 

This  boundary  la  shown  schematically  ir.  Figure  64.  It  was  felt  that  a greater  fraction 
of  the  velocity  profile  should  be  Included,  liowever,  and  the  boundary  at  which 


- 0.01 


max 


was  chosen.  It  is  shown  In  Figure  8 that  the  self-similar  radial  jet  velocity  profiles 
are  well  represented  by  the  universal  function: 


V 

= sGch^  [(1.763)  -^]  (91) 

m ■ 5 

and  this  equation  can  be  used  to  show  that,  if  Ax'''  denotes  the  99^  velocity  boundary: 
Ax*  ■ (3.394)  AXj 

which  in  Equation  (89)  results  in: 


(0.719) 


(92) 


Tills  equation  permits  calculation  of  the  radial  jet's  footprint  on  the  airframe 
undersurface.  Since  the  annular  jet  velocity  profiles  are  salf-slmllat , it  i.s 
advantageous  to  write  Equation  (87)  in  similarity  variables.  The  equation  is 
first  written  as: 


where 
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5 ' ■ X “ X 


Normalizing  the  force  by  the  total  annular  jet  moinentum  flux 


27TP  Vj 

and  norm.illzlng  all 
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by  R^: 
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and  defining 
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Air 


(94) 


d'7 


(;j)  (;^)(^)  r^(i) 
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Tt  has  been  found  In  Reference  8 that  the  maximum  velocity  in  the  Jet  may  be 
represented  by 


(2.496) 


4-A., 

MR  - Rj) 
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(95) 


Using  Equations  (90)  and  (95)  in  ,”-♦) 


^ ^ R. 


(96) 


where 


'■W  « 


Using  Equation  (91)  for  the  velocity  profile,  the  function  0 may  be  Integrated 
explicitly.  The  result  is 


G - (1.763)"^  tanh  E*  - y tanh^  r* 


where 


5*  ■ (1.76J)  5 


From  the  geometric  relation 
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the  y-integration  may  be  written  out.  The  final  result  for  the  fountain  impingement 
force  is: 
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2tt  J 


y ttanhe*  - -j  tanh^E* 
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and 


f97b) 


5*  - (1.763) 

[n  the  limiting  case  of  an  infinitely  large  plane  (f.^  •+  ®i  '■  ~®.  y '►•  “)  t it  may 
be  shown  that  Equation  (97a)  rnducea  to 


and  consequently  the  total  irapinpament  force  coefficient  on  an  J.nfinita  plane  would 
be 


C 


N 


(98) 


This  is  half  the  value  given  by  Equation  (74)  showing  that  the  normal  force 
coefficient  generated  on  an  infinite  piano  is  one  half  of  the  total  vertical 
momentum  flux  in  the  fountain.  This  implies  that  there  ia  a vertical  momentum 
loss  In  the  impingement  process.  The  difference  occurs  because  the  normal  dynamic 
pressure  I rather  than  the  normal  momentum  flux  has  been  converted  to  normal  force. 

The  dynamic  pressure  leads  to  batter  agreement  with  data  and  will  be  retained, 
although  it  is  realized  that  it  results  in  the  above  anomaly. 

Equations  (97a  and  b)  have  bean  incorporated  into  a computer  program  to  calculate 

the  fountain  impingement  force.  The  aircraft  planform  ia  subdivided  into  strips 

of  width  Ay.  At  a given  value  ol:  , Integrc  cr  ,n  limits  (x,j,  and  at  the  conter- 

llne  of  each  scrip  are  determined  from  a drawing  of  the  aircraft  planform  and 

fountain  footprint,  analogous  to  the  plan  view  in  Figure  64.  Equation  (97a)  and  the 

fountain  width  distribution  given  by  Equation  (89)  scale  with  R , so  as  R incraases , 

0 o 

the  aircraft  planform  shrinks  and  a smaller  section  of  the  annular  jet  ia  captured, 
thus  reducing  the  fountain  impingement  force.  The  integration  limits  are  Input  to 
the  program,  and  the  integral  in  Equation  (97n)  ia  calculated.  For  symmetric  situ- 
ations such  as  indicated  in  Figure  64,  tho  result  may  be  multiplied  by  2 to  obtain 
the  total  Impingement  force,  For  an  offset  fountain  where  the  axes  of  the  opposed 
Jets  do  not  coincide  with  the  \-axls  (y^  ■ 0),  the  computatltins  would  be  carried 
out  separately  for  each  side  y i and  y <.  y . The  hyperbolic  functions  in  the 
Integrand  of  Equation  (97a)  remain  constant  for  > f,*  , so  the  input  can  be 

simplified  by  selecting  tho  wing  trailing  edge  (or  loading  ed('e)  for  the  integration 
limits,  whenever  their  x coordinates  exceed  |Ax 

6.1.7  Calculation  of  Fmmt.sln  Buc.kdwwn  Force  - E.xaminHtion  of  data  olitalned 
from  the  test  program  to  be  described  In  Section  6.3.1  reve.Tls  that  the  cnlculaLcd 
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fountain  impingement  force  dooH  not  account  for  uli  the  forcca  acting  on  the  model. 

Although  the  calculation  predicts  the  shape  of  the  normal  force  curve  obtained  as 

the  fountain  impinges  along  different  fuselage  stations  fairly  well,  the  magnitude  of 

this  force,  and  in  some  cases  its  sign,  were  Incorrect.  Furthermore,  the  calculated  foun 

tain  impingement  force  decreases  as  the  distance  R^j  between  the  virtual  origin  and  the 

airframe  (see  Figure  64)  increases,  whereas  the  data  revealed  an  opposite  trend  over 

a considerable  range  of  R . Tl\e  reason  for  this  behavior  is  that  tha  fountain  Inducos 

o 

a auckdown  force,  in  addition  to  an  Impact  force. 

Figure  65  shows  a symmetric  two-jet  fountain  impinging  on  a rectangular 
plate.  As  bsfore,  the  fountain  has  bean  raplacad  by  an  equivalent 
radial  jet.  In  this  instance,  however,  tha  radial  jet  issues  from 
a slot  cut  in  the  ground  plana,  with  the  plate  located  at  a height  H above  this 
plana.  This  slot  Is  ohown  In  the  section  view  A-A.  The  radial  Jet  and  the  flow 
along  the  plate  bottom  entrain  mnaa  from  the  surrounding  medium.  Consider  a control 
volume  bounded  by  tha  plate  lower  surface,  the  projected  plate  .iroa  on  tha  ground 
plane,  and  the  lateral  area  of  the  volume  bounded  by  tltoso  two  areas.  Continuity 
requires  that  there  be  a net  mass  flux  of  external  air  into  this  control  volume,  to 
balance  the  mass  flux  entrained  by  turbulent  portions  of  the  flowfiold.  This  influx 
from  tha  surrounding  medium  is  responsible  for  reducing  the  pressure  on  the  bottom 
of  the  plate,  and  creating  a suckdown  force. 


Origin  of  Suckdown  Force 
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6.1.2il  Fountain  Suckdovm  Forcg  foi-  a Rectanaular  Plata  - An  approximate  means 
ol!  calculating  the  fountain  suckdown  force  has  been  devised.  The  calculation  ia 
carried  out  in  two  steps.  The  first  step  estimates  the  total  excess  mass  flux  in 
the  turbulent  flow  out  of  the  control  volume  described  in  the  previous  paragraph. 

The  second  step  estimates  the  pressure  in  this  control  volume  by  calculating  the 
Inflow  velocity  for  a one -dimensional  flow  which  supplies  the  excess  mass  flow  across 
the  control  volume's  lateral  area. 

To  calculate  the  mass  flow  entrained  by  the  turbulent  flow,  the  flow  turning 
geometry  discussed  in  Section  6.1  and  depicted  in  Figure  62  is  used.  The  flow 
direction  on  the  plate  bottom  will  than  appear  as  shown  in  the  plan  view  of  Figure 
65,  It  is  assumed  that  after  the  radial  Jat  is  split  by  impingement  on  the  bottom 
of  the  plate,  each  half  continues  to  entrain  aa  if  it  had  not  been  turned.  Con- 
sequently, the  total  maaa  flux  antralnod  on  each  side  of  the  radial  Jet's  plane  of 
symmetry  is  the  same  as  if  the  turning  caused  by  impingemeitt  had  never  occurred.  The 
total  excess  mass  flow  out  of  the  control  volume  is  therefore,  calculated  from  a dia- 
gram such  aa  that  shown  in  Figure  66.  This  diagram  la  obtained  by  "folding"  the 
plate  pianform  along  the  fountain  impingement  line  PQ  in  Figure  65.  For  the  sym- 
metric case  in  Figure  65,  the  diagram  on  both  sides  of  the  fountain  plane  of  symmetry 
is  the  Dame,  If  the  fountain  implngad  at  a distance  x^  from  the  plate  midchord,  then 
for  X > Xg,  the  plate  upper  edge  in  Figure  66  would  be  at: 


z 


x 


0 


I 


and  for  x < x^,  the  upper  edge  would  11a  at: 


To  compute  the  net  mass  flux  entrained  into  the  control  volume,  it  Is  necessary  to 
calculate  the  mass  outflow  normal  to  the  lines 


AB , HC , and  CD 

and  subtract  from  thiSitlie  mass  influx  across  the  line  AD.  Constrained  radial  Jet 
data  obtained  in  Reference  8 are  used  to  calculate  this  mass  flow.  Since  the  flow- 
field  is  symmetric,  it  is  easier  to  work  with  the  right  hand  side  of  Figure  66. 

Tlie  net  mass  flow  normal  to  the  line  OC  is  given  by 


ih 


1 


COS0 


'z 


dxdy 


Vj^  cosO 


dxdy 

" H 

0 


- h 


(99) 
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FOUNTAIN  MASS  FLOW  DIAGRAM  FOR  IMPINGEMENT  AT 
MIDCHORD  OP  A RECTANGULAR  PLATE 


and  normal  to  the  line  CD! 


il\2  ■ 


- 

/ / 

(Rg  - h)  o 


p sine]  dxdK 


(100) 


The  x-integratlona  can  be  carried  out  uelng  the  eelf-oimilar  properties  or 
the  radial  jot  velocity  profile,  and  the  data  of  Reference  8t 

(jO  qO 

/'k 'I- ■ '% 


Using  Equation  (91)  to  reprssent  the  correlated  velocity  profile,  this  bocomes 


104 


Substituti.ng  Equrttiona  (104)  and  (105)  or  (106)  into  (99) i 

^b/2  hj2 

(R  - H)/  , 


ilij  - (0.300)  (pV  ) 


J Ax  R '.(R  + § )/  *r-wi 

^ ^ V ^(R^  + f)'  + y' 


(R  “ H)‘ 
o 


TT-I 


and  carrying  ouC  Che  inCegrnClona: 


*1  ■ “j  { "■  t > *" 

(-sfV)  ‘ 


+ 1 


- (K^  - H)  In 


To  carry  out  the  Integraclon  In  Equation  (100).  note  that 


(107) 


fine 


- 


+ (b/2)' 


Using  this  relation,  in  addition  to  Equation  (104).  Equation  (100)  may  be  Integrated 
explicitly,  and  the  rufulc  iei 


I /b\  ^ ^ ^'’^2)2 

rii,  - (0.300)  pV.  Jix.  R,  f In  • — §■■— [ 

^ J ^ J J '2  l(R  -H)+  y(R  - H)^  + (b/2)^  )■ 


(108) 


Combining  Equations  (107)  and  (108),  the  total  excess  m.aHS  flux  out  of  one  quarter 
of  the  control  volume  is: 


Htj,  - (0.300)  pVj^4xj  R^  K 


(109)  1 


where 


(R  + v)  in 

0 z 


(r^W) 


+ 1 


(110) 
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+ f)  + \ 

/(R^  + c/2)^ 

+ (b/2)^ 

(Rq  ' H)^  + 

(b/2)^ 

ReferrinR  to  Flgura  65,  it  la  aaaumed  that  a one-dlmenalonal  flow  across  one  quarter 
of  th(3  control  volunve's  lateral  area  suppllaa  this  excess  mass  flow,  so 


where 


(111) 


(112) 


AsNumlnK  ambient  pressure  acts  on  the  plate  upper  surface,  the  suckdown  force 
will  be  given  by 


where 


end 


• ^Pb 


P ) A . 
^ W 


be 
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_ V ‘ 
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Using  Equation  (111),  Che  suckdown  force  then  becomes 


•*  (-) 


1/4 

P 


Substituting  from  Equation  (109)! 


(113) 


(U4) 
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F - - (O.OA5)  ^1/4  [p  V,^  (Ax,  R.)  K^] 

” i »»  J .1  J 


0 


i/4 


I'hii  total  auckdiwn  forco.  La  given  by 


F - 4 F 

a s 


1/4 


Defining  a force  coefficient  by  normalieing  by  the  total  radial  jot  momentum 
flux: 


(lib) 


• 2ifp  R,  Ax.V,^ 

J J J 


and  aubstituting  for  F 


C - - (0.0286)  '^'l/A 

^1/4 


Subatituting  for  the  planform  and  lateral  area 


- -(0.0286) 


be 


(b  + c)^ 


(117) 


where  K has  been  given  in  Equation  (110).  The  following  caae  ha:,  been  calculated. 

h - 20  in 

c ■ 8 in 

(R^  - H)  ■ 2 in 

R " 4 - 9 in 
0 

Foi  the  case  where  the  fountain  impinges  at  midchord,  the  reeulta  are  shown  in 
Figure  67.  The  fountain  auckdown  force  calculated  according  to  the  above  equations 
is  Hho:/n.  The  force  is  negative  and  increases  towards  aero  as  R^  increases.  The 
fountain  impingement  force  calculated  by  the  method  detailed  in  Section  6, 1.1,2 
la  also  shown . This  force  is  poaltive  and  decreases  in  magnitude  aa  R^  IncreaHuu . 
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FIGURE  67 

CALCULATED  FOUNTAIN  IMPINGEMENT  AND  SUCKDOWN 
FORCES  ON  A RECTANGULAR  PLATE 


Th«  reiu.lt  of  lunmilnfi  the  two  effect*  in  a positive  force  which  initlrlly  increases 
with  Rg,  and  then  dacrsases  as  the  negative  luckdown  contribution  approaches  zern, 

In  the  fountain  suckdown  caLculatloni  tl\o  radial  Jet  thickness  has  not  been 
taken  Into  account,  and  the  radial  Jet  "footprint"  has  not  bc«n  considered,  as  was 
don*  for  the  ImpingeTnant  force.  Furthermore,  it  has  bean  assumed  that  the  ambient 
air  influx  occurs  acroes  the  entire  lateral  area  of  the  control  volume.  In  reality, 
part  of  this  area  is  occupied  by  the  turbulent  flow,  While  refine- 
ments to  the  method  to  take  account  of  thaaa  factoru  could  bu  dovelopod,  the  main 
objectlva  has  been  to  obtain  approximate  anawcrs  that  yiold  the  correct  trends, 

6. 1,2. 2 Fountain  Suckdown  Forcu  For  a Tr.apezoldal  Wiua  ~ FlRuro  00  ahowa  the 
geometry  relevant  to  tlia  calculation  of  the  fountain  suckdown  force  for  a trapezoidal 
wing.  In  the  sketch  of  tlie  wing  geometry,  the  radlul  .jet  representing  the  fountain 
is  shown  to  impinge  at  n dlstanct*  x^  from  the  origin.  As  in  the  rectangular  plate 
case,  the  mass  flow  dingrama  are  obtained  by  "folding"  the  wing  planform  along  the 
iritarsBction  line  with  the  fountain  plana  of  .symiiiotry,  Thu  excess  mans  flux  out  of 
ti(G  control  volume  bounded  by  the  wing  undersurf  .icei , I ts  projection  on  the  ground 
piano,  and  the  lateral  area  between  those  .surfaces,  is  calculated  from  the  mass  flow 
diagrams  68a  and  6Bb,  Figure  68a  Indicatas  the  geometry  for  x ’ x^.  In  this  ease, 
the  mass  outflow  is  calculated  from  one-half  of  the  radial  Jut  niass  flow  across  the 

in<) 


boundary  A'B'C'D'E'F'G'H' , and  the  Inflow  across  the  boundary  A'H'  is  subtracted  out. 
Note  that  the  shaded  areas  in  Figure  68b  do  not  contribute  to  fountain  entrainment, 
since  according  to  the  assumed  radial  Jet  impingement  process  (Figure  62),  these  areas 
are  not  wetted  by  the  flow. 

The  calculations  proceed  as  in  the  rectangular  plate  case,  except,  for  the  outflow 
at  the  swept  wing  leading  edge.  Referring  to  Figure  68a,  the  mass  outflow  across 
the  Hue  CD  is  given  by* 


rj 

0 o 


pVj^  C0,i  (y“0)  dxds 


Since 


is  cos  Y “ dy 

yf 


/ / 


PV 


cos  (y-6)  dxdy 
R COB  Y 


(118) 


Using  Equation  (104)  for  the  x integrationt 


tfi  - (0.300)  pV, 
e j 


{“’‘i  / 


cos  (v-6) 
cos  Y 


dy  (119) 


The  equation  describing  the  swept  leading  edge  is: 


z - - (tan  y)  y + (R^  + x^) 
Cctnsequently , along  the  leading  edge. 


(120) 


cos  9 


- (tan  y)  y + (R  + x^) 

O t 


|[(tan  y)  y - (R^  + X|.)]‘^  + 


and 


(321a) 


sin  9 


|[(can  Y)  y - (Rj^j  + x^)]^  + y^^ 


(12.1b) 


Expanding  the  cosine  term  in  Equation  (118)  and  using  Equations  (121) 

yf 


Ag  - (0.300)  pVj  ^Xj  Rj  |(R^  + x^)  J 


Jl. 


1 cos 


2^,  - 2tany  (R^  + y + (R^  + x^.)' 
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and  carrying  out  the  integratloni 


■»  (0.300)  pVj 


(.121) 


where 
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^ Yj.  - a li\Y  C03  Y (R  + Xj.)  + 

r~ 


Yj.  - 2 sinY  cos  Y (R^  + x^) 


coa  Y (Rq  + x^)  (I  - sin  y) 


+ cos^Y  (Rjj  + Xj.)^ 


1/2 


(123a) 


and  y , 


(123b) 


tanY 


Equations  (122)  and  (123).  in  addition  to  maEis  outflow  and  Inflow  results 
along  the  ocher  boundaries  indicated  in  Figure  68,  have  been  used  to  calculate, 
the  auckdown  force  due  to  radial  jet  impingement  on  the  undersurface  of  a trape- 
zoidal wing  planform.  The  results  are  presented  in  Section  6,3,  where  they  are 
compared  to  data  from  the  fountain  Impingement  test  program  described  therein. 

The  methodn  for  calculating  fountain  implngemant  and  auckdown  forces  for  two- 
jet  symmetrical  fountains  have  been  extended  to  multiple-jet  fountain  configurations. 
The  general  technique,  and  its  application  to  a four-jet  fountain  configuration, 
are  detailed  in  the  next  Section, 

6 . 2 Multiple-Jet  Fountain  Impingement  and  Suckdown 

As  discussed  in  Section  5,3,1,  the  various  fountain  segments  which  make  up 
three  or  four-jet  fountain  configurations  are  subdivided  into  "inner  region  sectors" 
and  "outer  region  sectors".  The  inner  region  sectors  Interact  to  form  a central 
fountain,  which  is  not  currently  well  defined.  However,  it  has  been  argued  in 
Section  5.3.1  that  in  the  inner  region,  the  total  vertical  momentum  flux  across  planes 
parallel  to  the  ground  plane  must  be  conserved,  provided  the  .■’.entral  fountain  is 
confined  to  the  inner  region.  Wlren  calculating  a multiple-jet  fountain  impingement 
force,  it  is  necessary  to  consider  the  footprint  of  the  inner  and  outer  fountain 
sectors  on  the  aircraft  planform.  Because  of  conservation  of  vertical  momentum  In 
the  Inner  region,  as  long  as  the  inner  region  footprint  lies  wholly  within  the  air- 
craft planform,  the  fountain  impingement  force  due  to  the  inner  region  segments  will 
be  Independent  of  height  above  the  ground  plane.  The  outer  region  sectors  will  behave 


as  In  Che  two-jec  fountain  casei  contributing  a decreasing  amount  to  the  impingement 
force  as  the  aircraft  height  above  the  ground  plane  Increases. 

The  calculation  of  fountain  suckdown  forces  for  multiple-jet  conf Iguraclonn  Is 
greatly  complicated  by  the  presence  of  a central  fountain  and  interaction  between 
various  fountain  segments  after  they  Impinge  on  the  airframe  undersurface.  For 
example,  Figure  69  depicts  an  idealized  flow  pattern  caused  by  Impingement  of  n 
four-jet  fountain  on  the  undersurface  of  a square  plate.  The  planes  of  symmetry 
of  the  four  radial  Jets  which  represent  the  fountain  intersect  the  plate  along  the 
lines  OA,  OB,  OC,  and  OD.  The  segments  00^2i  ®®23’  *^*^41 

the  inner  sectors,  and  the  square  inner  region.  The  idealized 

turnln[.;  process  illustrated  In  Figure  62  yields  the  streamline  pattern  depicted  by 
the  dotted  lines  in  Figure  69.  This  overall  pattern  has  been  highly  idealized, 
for  the  fallowing  reasons: 

1.  In  the  neighborhood  of  the  point  0,  central  fountain  impingement  occurs, 
leading  to  an  unknown  atceamline  pattern, 


FIGURE  69 

IDEALIZED  STREAMLINE  PATTERN  FOR  A SQUARE,  FOUR-JET 
FOUNTAIN  IMPINGING  ON  THE  BOTTOM  OF  A FLAT  PLATE 
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2i  Along  the  diagonal  OPj_,  OP21  OP^,  and  OP^,  the  flows  along  the  plate 
undersurface  interact  with  each  other. 

The  central  fountain  and  the  flow  near  diagonals  OPj , OP2,  OP3,  and  OP^^ 
may  well  develop  as  turbulent  flows  with  ontrainmont  charactoristica  which 
differ  from  those  for  the  radial  jets  which  generate  tliem.  These  entrain- 
ment characteristics  are  completely  unknown,  so  for  computational  purposes, 
the  following  assumption  will  he  made: 

0 The  central  fountain  and  the  flows  generated  by  Interaction  of  flows  along 
the  plate  undersurface  do  not  entrain  any  excess  ambient  air. 

iJith  the  above  assumption,  the  raDial  let  entrainment  rate  is  assumed  to  exist 
along  each  one  of  the  sectors  BOP^,  B0P2»  etc. 

Four-jet  fountain  implngament  and  auckdown  forces  have  been  calculated  for  the 
configurations  in  the  program  detailed  in  Section  6.3.  A description  of  these 
calculations  is  Included  in  Che  sections  following. 

6.2.1  Calculation  of  Four-Jet  Fountain  Impingement  Force  - Figure  70  shows  the 
planform  of  the  model  tested  in  the  fountain  Impingmenc  experiment  described  in 
Section  6.3,  and  a front  view  of  the  modal  orientation  relative  to  two  of  the  four 
radial  jets  which  simulate  the  fountain.  The  planes  of  symmetry  of  the  four  radial 
Jets  intersect  the  plane  of  the  model  planform  along  lines  OA,  OB,  OC,  and  OD . The 
eight-inch  square  defined  by  the  radial  jet  centerlines  bounds  the  inner  region. 

For  the  cases  calculated  the  inlet  region  center  lias  directly  below  the  model 
reference  point.  Calculations  have  been  carried  out  for  varying  values  of 

R^  ■ 5,  6,  7,  and  8 inches  , 

In  Figure  70,  the  inner  and  outer  region  radJ.a],  jet  sectors  are  shown  by  different 
crosshatchiug.  The* 


V„/V 
R max 


O.Oi 


footprint  defined  in  Section  6.1. 1.2  is  shown.  For  the  outer  region  sectors,  the 
impingement  force  la  calculated  by  the  method  outlined  in  Section  6, 1.1. 2.  For  the 
inner  region  sectors,  it  is  assumed  that  the  entire  inner  region  f.iotprint  is  cap- 
tured by  rhe  planform.  Consequently,  the  conservation  of  total  vertical  momentum 
flux, discussed  in  Section  6,2,  permits  calculation  of  the  inner  region  impingement 
force,  without  requiring  det.nilod  knowledge  of  the  central  fountain 'n  flowfit.‘ld, 

The  inner  region  force  la  equated  to  one-iialf  of  the  vertical  momentum  flux 
in  the  inner  jet  sector,  computed  at  z * S/2  (see  the  sketch  in  Figure  70).  From 
Equation  (73)  in  Section  5.3.1,  the  force  transmitted  by  each  of  the  inner  region 
sectors  Isi 
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For  n square  jec  arrangementj  0^^  “ so  lor  all  four  Inner  region  sectors, 

since  the  Jet  radial  momentum  fluxes  are  all  equal! 


(125) 


The  total  fountain  Impingement  force  is  then  calculated  as  follows : 


F -F  +2F  +F  +F 

*total  *^total  *wing  ^forward  *aft 


where 

F 

wing 


- normal  force  induced  on  one  wing  panel  by  the  fountain  outer 
sector . 


F 

forward-  Impingement  force  induced  on  forward  fuselage  by  outer  sector 

F - impingement  force  Induced  on  aft  fuselage  by  outer  sector, 

aft 


The  impingement  normal  force  coefficient  is  then  defined  as  follows! 


^*total 


(126) 


This  normal  force  coefficient  is  combined  with  a suckdown  force  coefficient  to 
obtain  the  total  force  acting  on  the  model,  Thr  calculation  of  fountain  suckdown 
forces  is  outlined  in  the  next  Section. 

6.2.2  Calculation  of  Four-Jet  Fountain  Suckdown  Force  - Figure  71  shows  the 
idealized  flow  pattern  created  when  a four-jet  fountain  strikes  the  undersurface 
of  a trapezoidal  wing  planform.  The  wing  dimensions  correspond  to  the  theoretical 
planform  of  the  model  tested  in  the  fountain  impingement  experiment  described  in 
Section  6.3.  For  simplicity,  the  presence  of  the  fuselage  Is  neglected  in  these 
computations.  Figure  71  shows  one-half  of  the  flow  pattern  derived  for  the  square 
plate  case  in  Figure  69,  superimposed  on  the  trapezoidal  wing  planform.  The  inner 
region  is  defined  by  the  8 in.  square.  As  described  in  Section  62,  it  ic-.  assumed 
that  the  central  fountain  and  the  flow  interactions  near  lines  OB,  OD,  OF  and  OH 
do  not  entrain  any  additional  flow.  Consequently,  these  lines  constitute  boundaries 
at  which  the  total  entrained  mass  flow  is  to  be  computed. 
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IDEALIZED  FLOW  DIRECTIONS  ON  TRAPEZOIDAL  WING  UNDERSURFACE 
AFTER  FOUR-JET  FOUNTAIN  IMPINGEMENT 


Fountain  masa  flow  dlagi'anm  analogoua  to  those  shown  for  two- jet  fountain  suck- 
down  in  Figures  66  and  68  are  set  up.  In  this  casei  however,  the  wing  planform 
is  to  be  "folded"  about  the  two  lines,  AOE  and  GOC . The.  fold  along  AOE  is  performed 
to  calculate  the  entrainment  from  the  fore  and  aft  radial  jets,  whose  planes  of 
symmetry  go  through  this  line<  The  fold  along  GOC  is  performed  to  calculate  the 
entrainment  from  the  left  and  right  radial  jets,  whose  planes  of  symmetry  go  through 
this  line.  Figure  72  indicates  the  mass  flow  diagrams  in  the  x-s  plane,  with  the 
wing  planform  folded  along  AOE,  The  ground  plane  is  located  at  a distance  H below 
the  model.  For  the  rugion  located  at 


I 


i 


X > 0, 

the  entrained  mass  flow  la  computed  as  followii: 

Entrained  mass  flow  • Mass  outflow  at  A'ABOn'  - Masn  Inflow  at  A' O' 
For  the  region  located  at 


X < 0 

Entrained  mass  flow  ■ Mass  outflow  at  O'ODEE'  - Mnas  inflow  at  O'l!' 

These  mass  fliixea  are  calculated  by  formulae  analogous  to  thoae  used  for  single 
radial  jets  striking  rectangular  and  trapezoidal  wings,  hb  described  In  Section  6.1.2. 
The  suckdown  forces  acting  on  each  portion  of  the  wing  are  calculated  by  requiring 


z 


I an  inflow  velocity  acroaa  each  lateral  area  to  supply  the  entrained  mass  flow* 

I ■. 

I and  relating  the  velocity  to  a pressure  as  described  in  SectiDn  6il.2.  The  calculations 

I 

[ ‘ associated  with  Figure  72  permit  an  estimate  of  suckdown  forces  for  wing  area  segments 

I ABO,  ODE  (sea  Figure  71),  and  their  mirror  images  AHO  and  OFE . To  calculate  the 

f rest  of  the  suckdown  force,  the  mass  flow  diagrams  in  Figure  73  are  used.  For  one 

!-  ..  half  of  the  section  at 

i- 

\ X > 0 

; Entrainrd  mass  flow  ■ Mass  outflow  at  O'OBCC'  - Mass  inflow  at  O'C' 

.i  and  for  one  half  of  the  section  at  x < Qi 

1 « Entrainod  roaea  flow  ■ Maas  outflow  at  O'ODD"  CC'  - Mass  inflow  at  O'C 

I;  Note  that  (as  is  the  single  radial  jet  case)  , part  of  the  wing  area  is  not  wetted 

1 by  the  impinged  flow,  and  may  then  be  left  out  of  the.  calculation.  The  calculfttiona 

I associated  with  Figure  73  permit  an  estimate  of  suckdown  forces  for  wing  area  sag- 

r ' 

!'  ments  OBC  and  CD"D0  (.see  Figure  71)  J and  for  their  mirror  images  OHG  and  GF"FO. 

.After  .;he  lengthy  computations  described  above  are  completed,  the  suckdown 
I j forces  for  each  part  of  the  wing  area  are  summed,  and  a total  suckdown  force 

f-  C,  coefficient  is  calculated.  The  results  obtained  are  compared  to  data  in  Section  6.3, 

I. 
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6.3 


U.-)tnPaPlaon  of  Pradicted  and  Experimental  Fountain  Forcaa 

Af,  parK  of  tho  MGAIR  XKAD  program  described  in  Rci'srance  23,  a fountain  gener- 
ator was  designed  and  fabricated.  This  apparatus  is  capablo  of  simulating  the 
fount.iin  gencruted  by  a pair  of  circular,  perpendicular  jets  of  equal  tl\rust,  or 
by  four  circular  perpendicular  Jets  of  equal  thrust  ar r.:,n:.;efJ  in  a square  pattern. 

Tile  simulated  fountains  were  allowed  to  impinge  on  a representative  aircraft  mouel, 
and  the  forces  and  moments  induced  on  the  model  wi.'re  measured  by  a a ix-component 
strain  gauge  balance.  After  a brief  description  of  the  apparatus  and  teat  setup, 
measured  normal  forces  are  compared  to  the  results  predicted  using  the  methods 
described  above. 

6.3.1  Fountain  Impinaement  Teat  Apparatus  - The  fountain  generator  design  was 
based  upon  the  analogy  between  a symroatric  two-Jet  fountain  and  a radial  Jet,  and 
the  postulated  equivalence  between  multlple-Jet  fountains  and  combinations  of  radial 
Jets.  The  fountain  generator  was  used  to  Isolate  the  fountain  forces  from  the  free 
and  wall  Jet  suckdown  forces  which  would  exist  If  the  entire  flowfield  were  simu- 
lated at  once.  Figure  74  is  a schematic  of  the  foumialn  generator  assembly.  Four 
pipes  are  supplied  with  high  pressure  air  from  an  existing  facility.  Each  pipe 
discharges  into  one  of  four  plenums,  which  are  equipped  with  porous  plates  and  screens 
to  smooth  the  flow.  Each  plenum  supplies  two  nosales.  The  four  plenums  and  eight 
nozzles  are  arranged  in  a square  pattern,  creating  four  annular  nozzles  through  which 
the  air  exhausts,  forming  four  at.nular  Jets.  To  simulate  two-Jot  fountains,  three 
of  the  annular  nozzles  were  capped.  To  simulate  four  Jet  fountains  all.  the  nozzles 
were  operated.  As  explained  In  Section  5,  only  a fraction  of  each  annular  Jet's  mass 
flux  should  contribute  to  the  four-jet  fountain.  Consequently,  a detachable  surface 
In  the  shape  of  a pyramid  was  used  to  split  the  annular  Jets  whon  simulating  a four- 
jet  fountain,  as  indicated  in  Figure  74.  The  forces  and  moments  produced  by  the 
simulated  fountain  Impingement  on  a wing-fuselage  model  wero  measured  by  an  internal 
six-component  strain  gauge  balance.  During  the  first  phase  of  the  test  program, 
a tot.ii  pressure  rake  was  used  to  probe  the  Jet  flowfields.  Figure  75  shows  the 
overall  test  arrangement,  including  the  fountain  generator,  model,  rake,  and 
support  structures.  The  model  could  be  translated  in  three  directions  relative 
to  tlia  fixed  fountain  generator,  and  rotated  about  all  three  axes.  lUah  and  low 
wing  configurations  ware  tasted,  with  and  without  trailing  edge  flaps  und  simulated 
underwlng  stores . 

Comparison  of  Two-Jet  Fountain  Results  - An  examination  of  tile  data  ob- 
tained from  this  test  program  soon  revealed  the  necassitv  to  Include  a fountain 
suckdown  force.  For  a alngle  radial  Jet,  Figure  76  compares  some  test' results  to 
a calculation  of  the  fountain  impingement  force  carried  out  by  the  method  described 
in  Section  6. 1,1, 2.  For  the  cases  shown,  a single  radial  Jet  impinged  upon  a clean. 
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FIGURE  76 

CALCUL  ATED  JET  IMPINGEMENT  FORCE  COMPARED  TO  MCAIR  IRAD  TEST  DATA 

SINGLE  ANNULAR  JET,  LOW  WING 
(NEGLECTING  FOUNTAIN  8UCKDOWN) 


low  wing  model  configuration,  with  the  plane  of  the  annular  nozelo  oriented  perpendi- 
cularly Co  the  model  plane  of  ayrametry.  Thi«  arrangement  aimulataa  Che  fountain 
for  a two-lif c-jet  aircraft  with  the  J«tH  mountod  in  the  CuKelago,  Tim  curves 
Indicntu  tlie  variation  in  normal  forcu  coulficLunt  as  the  Jut  impingos  at  different 
Stations  along  the  axis.  As  In'I'ori!,  the  normal  force  coo iri.i  lent  is  deCinod  as: 

,,  Wuiaiiul  force  on  model 

N Annular  lot  monu'ntiim  flux 

For  the  four  Jet  exit  vi'loi:  li  loa  tostod,  tho  ronults  Indicate  that  fairly  good  rlnta 

correlation  is  obtained  using  this  coefficient.  Kach  one  of  the  graphs  corresponds 

to  a different  value  of  the  distance  (r  ) between  the  nozzle  centerline  and  the 

o 

reference  point  on  the  bottom  of  the  fuselaga.  It  is  obvious  that  the  calculated 
Jet  implngamonc  force  does  not  account  for  all  the  furces  acting  on  tho  model.  For 
R^  ■ 8 and  6 in.,  tho  Jet  impingement  calculation  predicts  the  shape  of  the  curve 
fairly  well.  A constant  force  is  created  when  the  Jet  impinges  on  the  fuselage. 

A sharp  increase  in  normal  force  occurs  as  tho  Jet  starts  to  Impinge  on  the  wing 
trailing  edge,  rising  to  a peak  value  whan  lies  In  the  region  of  m.aximum  planfurm 
span.  Thon  a suhsoquont  slower  decay  ns  the  Jet  Impinges  on  a region  of  reduced 
wing  span  containing  the  swept  wing  leading  edge.  Tlic  level  of  the  forces,  and  in 
some  cases  tlielr  sign  1s  not  accurately  predicted  by  the  Impingament  culculation . 

It  appeavs  that  In  addition  to  the  Implngemant  force,  a substantial  fountain  suck- 
down  force  ia  also  .icting  on  the  modol,  This  is  evident  in  IMgura  77,  which  is  a 
croasplot  of  the  data  iu  Figure  76  for  three  Jet  implni’.empnt  stations  One  would 
expect  that  as  R^  Increasoa,  the  Impingement  force  on  the  model  would  decrensu,  ainci? 

a smHlier  sector  of  the  annular  Jut  impinges  on  thn  model.  The  data  Imwaver,  sliow 
the  opposite  crenrt.  In  fact;,  for  x^  ■ -7  tho  force  on  the  racjdel  Is  negative  for 
all  values  of  R^.  Since  the  total  fountain  force  must  eventually  approach  zovo, 
all  th«  curves  in  Figure  77  must  eventually  approach  zoro,  but  tliio  bv.huvlor  was 
not  yet  evidunt  within  the  R^  range  toatocl. 

Wlien  the  fountain  suckdown  force  coefficient  is  innlude.rl,  Che  prediction  is 
substantially  Improved,  As  indicated  in  Figure  67,  for  a rectangular  plate,  the 
fountain  suckdown  force  is  negative,  and  approaches  zero  as  R^  Increasus.  Flguro 
78  compares  tlie  calculated  total  normal  force  coefficlont,  for  the  trapozoidal-w Ing 
model,  CO  tl:t>  data,  Tho  variation  with  R is  IndiciiLod,  fur  two  valuen  of  t.lm.  let 

Impingement  station  x . The  calculated  nonnal  force  coefficient  includeH  a suckdown 

0 

term  obtained  by  the  method  outlined  in  Section  6.  1.2. 2,  Although  l,lu>  aRrcenunit 
between  data  and  theory  is  far  from  perfect,  the  theory  predicts  the  data  treiuls 
reasonably  well.  The  agreement  could  perhaps  he  improved  by  lofin.lug  the  suckdown 


FIGURE  77 

VARIATION  OF  NORMAL  FORCE  WITH  HEIGHT  ABOVE  JET  ORIGIN 


calculation  (e.g.)  including  radial  jet  thickness  effects,  but  there  ia  a more  funda- 
mental factor  which  obscures  comparisons  to  data  from  this  test  program.  As  indicated 
in  Figure  75,  there  was  no  simulation  of  the  ground  plane  during  the  test.  The  cal- 
culation assumes  that  there  In  a ground  plane  which  constrains  the  area  through  v.diirh 
the  entrained  mass  flow  must  be  supplied.  For  the  computed  results  shown  in  Figure 
78,  it  was  assumed  that  the  ground  board  was  loc.ated  flush  with  the  top  of  the 
fountain  generator  plenums  (the  pyramid  was  removed  for  this  single-jet  configuration) 
but  this  It.  only  a crude  approximation  of  the  actual  geometry.  Unfortunately,  this 
problem  must  remain  unsolved  at  present.  A test  program  is  currently  planned,  which 
will  include  the  ground  board  and  allow  the  radial  jets  to  exhaust  through  appropriate 
slots  in  this  ground  board.  This  will  provide  a more  realistic  simulation  of  the 
actual  conditions,  and  will  permit  a more  direct  comparison  of  theoretical  and  exper- 
imental results. 

6.3.3  Comparison  of  Four-Jet  Fountain  Results  - The  calculation  of  four-jet 
fountain  Implngemint  and  suckdown  forces  has  been  detailed  in  Sections  6.2,1  and  6.7.2 
Results  were  obtained  specifically  for  the  configuration  tested.  Figure  79  compares 
the  predicted  force  coefficient  to  the  data.  As  indicated  by  the  sketch  in  the 
figure,  the  four  radial  jet  planes  of  symmetry  go  through  the  model  x and  y axes. 

The  agreement  between  theory  and  data  is  reasonably  good  in  this  case,  which  may  be 
due  to  the  fact  that  the  fountain  generator  pyramid  constrained  the  flow  underneath 
the  model.  When  performing  the  auck<lown  calculations  outlined  in  Section  6.2.2, 
it  was  again  assumed  that  the  ground  plane  was  flush  with  the  tops  of  the  fountain 
generator  plenums. 
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FIGURE  79 

COMPARISON  BETWEEN  THEORY  ANO  EXPERIMENT  • FOUR  RADIAL  JETS 
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7.  EFFECTS  OF  PROTUBERANCES  ON  INDUCED  FORCES 

Frocuberancee  on  the  underside  o£  an  airframe  may  be  used  to  enhance  the 
positive  fountain  force,  provided  they  are  judiciously  placed.  Theae  protuberances 
may  be  underwing  stores  and  pylons,  landing  gear  doors,  speed  brakes,  deflected 
flaps,  or  simply  plates  located  expressly  for  increasing  the  overall  lift  on  the 
airframe.  Methods  to  estimate  the  effects  of  underwing  stores  and  deflected  flaps 
on  the  fountain  force  have  been  derived  as  part  of  this  study.  They  are  exten- 
sions of  te.chniques  to  calculate  the  fountain  impingament  force  which  have  been 
described  in  Section  6. 

It  was  shown  in  Section  6 that  the  total  fountain  force  ie  made  up  of  a 
positive  impingement  force  and  a negative  suckdown  force.  It  is  assumed  that  the 
presence  of  underwing  protuberancee  doaa  not  affect  the  suckdown  fores,  but  doss 
product  an  inersmant  In  tho  impingement  force.  For  underwing  storea  and  for  a 
trai.llng  edge  flap,  the  calculation  of  uase  inevemants  is  outlined  in  the  follow- 
ing subsections.  Comparisons  to  axparimantal  results  from  the  fountain  impingement 
test  described  in  Section  6.3  are  also  shown.  The  methods  apply  only  to  two-jat 
configurations  but  they  could  be  combined  to  include  multiple-jet  configurations, 

7 . 1 Fountain  Fores  Increment  Dus  to  Stores 

In  Section  6.1.1,  a upecific  type  of  flow  turning  procass  has  batm  postulated 
to  occur  when  a radial  jet  etrikae  the  underside  of  the  airframe.  In  this 
turning  procese,  depicted  in  Figure  62,  the  velocity  component  normal  to  the  air- 
craft platform  is  rotated  by  90*,  and  the  tangential  component  is  unchangad.  The 
flow  produced  on  the  underside  of  the  eirfreme  when  the  redial  Jet's  plane  of 
symmetry  is  perpendicular  to  the  aircraft  plana  of  symmatry  than  haa  a velocity 
component  in  the  apanwlee  direction.  This  apanwise  flow  Is  inhibited  by  the  pre- 
sance  of  wing  pylons,  and  an  increase  In  fountain  force  results. 

Figure  80  shows  a drawing  of  the  aircraft  modal  used  for  the  fountain  impinge- 
ment test  program,  with  chordwise  plates  mounted  on  the  wing  undereurfacc  to 
simulate  underwing  pylons  and  stores.  Whan  the  model  undersurface  is  clean,  the 
fountain  Impingement  force  haa  been  calculated  by  equating  the  normal  component 
of  the  jet  dynamic  praeaure  to  the  preaeure  difference  acting  on  the  aircraft 
planform.  Thia  pressure  was  taken  to  act  over  the  radial  jet's  "footprint"  cn 
the  planform.  Since  the  chordwise  plates  restrict  the  jet's  spanwisa  flow,  it  is 
assumed  that  for  the  jet  sector  subtended  by  the  distance  between  the  plates,  the 
total  jet  dynamic  pressure  ie  recovered  as  a positive  pressure  on  the  planform 
undersurface.  In  Figure  BO,  this  occurs  over  the  part  of  the  fountain  footprint 
indicated  by  the  double  crosahatchlng.  Note  that  as  Rg  increases,  the  angle 
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FIGURE  80 

CALCULATION  OF  NORMAL  FORCE  INCREMENT  DUE  TO  UNDERWINQ  STORES 


■ubtandid  by  cha  diatance  batween  tha  atoraa  dacraaaea,  so  that  tha  normal  forca 
incrnmant  will  dacraaaa. 

Utidar  the  above  aesumptiona,  tha  normal  forca  Increment  acting  on  a planform 
atrip  dy  will  bet 


K 1 

) u f £.  V ' 

^ J 2 R 


^ (1  - ooa^6)  dx  dy 


d(AF^)  - ly*  Bin^  dx  dy 


Iha  x-lntegratlon  can  ba  carried  out  aa  detailed  In  Section  6. 1.1. 2.  Ualngt 


aln  6 


TfrP  ■ TTt? 


an  exprasaion  a%ulvalent  to  Equation  (97a)  can  be  derived  lor  the  total  iaplnge- 
meut  force  Increment  due  to  one  of  the  chordwlae  platesi  The  result  la: 


AF 


■ ^ /"  -^T57T 

It  “'o  (.1  +y 


tanh  S*  ” ^ tanh^f.* 


dy  (128a) 


where? , 


X - X 


5*  « (1.783)  [-7^] 


y ■ y/R 


(128b) 

(128c) 


and  the  other  variablea  are  defined  in  Section  6.1. 1.2.  For  the  configuration 
tasted  during  the  fountain  impingement  teat  program,  x^  ■ 3 in.,  and  the  entire 

fountain  footprint  was  captured  in  the  chordwisa  direction.  Then,  taking  the 
limit  in  Equation  (128a)  an 
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Defining  the  total  normal  force  incramant 


AC„  ■ 2AF 
^t  *at 


and  carrying  out  the  integration  in  Equation  (129) 
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For  the  case  te»t«d 

- 10  in. 

and  Equation  (130)  iu  plotted  In  Figure  81  along  with  data  from  the  fountain 
irapingeraent  test.  In  Figure  81,  the  calculated  incremental  normal  force  coeffi- 
cient is  added  to  the  clean  configuration  data,  to  obtain  the  solid  line.  This 
is  to  be  compared  with  the  data  points  obtained  with  the  simulated  stores.  For 
the  case  ahown,  the  chordwiae  plates  extended  from  the  wing  leading  edge  to  72% 
of  the  local  wing  chord.  The  agreement  between  data  and  experiment  is  quite 
satisfactory  except  for  the  lowest  value  of  R . 


Dlitancs  Irom  Nozile  ({.,  Rp  • In, 

FIGURE  81 

THEORETICAL  AND  EXPERIMENTAL  NORMAL  FORCE  COEFFICIENT 
WITH  SIMULATED  STORES 

7 . 2 Fountain  Force  Impinaement  Due  to  Trailing  Edae  Flap  Deflection 

Figure  82  Illustrates  the  reasoning  used  to  predict  the  increment  In  fountain 
Impingement  force  produced  by  deflecting  a trailing  edge  flap.  Once  again,  the 
geometry  corresponds  to  that  tested  In  the  fountain  impingement  axperimenuai 
program.  The  effect  of  the  flap  la  to  redirect  the  flow  along  the  wing  bottom 
surface  downward,  and  in  ao  doing,  to  creata  an  upward  reaction  force  on  the  air- 
frame. For  the  case  considered,  the  flap  extends  along  the  entire  exposed  wing 
span,  the  flap  hinge  line  is  at  72%  of  local  chord,  and  the  flap  deflection  angle 
ie  40”.  As  before,  the  radial  jet  1h  assumed  to  split  and  "fold"  upon  Impingement 
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FIGURE  82 

CALCULATION  OF  NORMAL  FORCE  INCREMENT  DUE  TO  FLAP 


on  the  wing  undeeaurfece . The  normel  force  la  then  related  to  the  momentum 
flux  normal  to  the  wing  trailing  edge. 

A "momentum  flux  diagram"^  analogoua  to  the  "maae  flow  diagrama"  uaed  to 
calculate  fountain  auckdown  forcea  in  Section  6 la  ahown  in  Figure  83.  Referring 
to  Equation  (72),  the  momentum  flux  normal  to  the  wing  trailing  edge  ie  given 
by 

M ■ [sin  - ain  6^]  • (135) 

A factor  of  (1/2)  haa  been  introduced  to  account  for  the  fact  that  only  one  half 
of  the  radial  Jet  momentum  la  deflected  towarda  the  trailing  edge,  In  addition, 
another  factor  of  1/2  ahould  be  included  to  account  for  the  fact  that,  aa 
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MOMENTUM  FLUX  DIAGRAM  FOR  COMPUTINQ  FOUNTAIN  FORCE 
INCREMENT  DUE  TO  TRAILING  EDGE  FLAP  DEFLECTION 


diicuiiid  In  Sactlon  th«  dynamic  praaiura  Impact  relation  Impllae  that 

one  half  of  the  radial  Jet's  momentum  la  lost  whan  It  Implngeo  on  the  wing  under- 
aurfaca.  On  the  other  handt  Equation  (13S)  would  have  to  ba  multiplied  by  2 
to  account  for  tha  other  half  of  the  aymmetrlc  momentum  flux  dlagrami  Conae- 
quently,  Equation  (135)  will  be  taken  to  repreaant  the  momentum  flux  neaded  to 
calculate  tha  total  fountain  force  Increment.  Introducing  tha  momentum  component 
normal  to  the  planform,  the  normal  force  increment  ia  given  by 

AF^  - ain  6^  [ain  0^  - ain  ^ 

So,  in  tarma  of  a normal  force  coefficient 

Bin  i. 

£•  [ain  - Bin  8^]  . (136) 

Note  chat  in  Figure  83  the  flap  only  extenda  acioan  the  exposed  wing  span.  The 
angles  9^  and  0^  arei 
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° J(b/2)^  + (R^  + 2.7)' 


(137b) 


There  le  an  extra  force  increment  due  to  the  momentum  flux  acrosa  the 
flap  tip  I but  this  term  hae  been  neglected. 

Equations  (133)  and  (137)  ware  used  to  calculate  an  Incrementaj.  normal 
force  coefficient  due  to  flap  deflection.  This  Increment  wa.)  added  to  the 
clean  configuration  normal  force  coofflclenti  and  the  result  Is  Indicated  by  the 
solid  line  In  Figure  84.  The  figure  also  Includes  the  experimental  results 
obtained  with  the  defected  flap.  Agreement  between  theory  and  data  is  satlof ac- 
tor y,  although  the  theory  appears  to  overpredict  the  force  at  tha  higher  values 
of 
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FIGURE  84 

THEORETICAL  AND  EXPERIMENTAL  NORMAL  FORCE  COEFFICIENTS 
WITH  PULL-SPAN  FLAPS 
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8.  SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS  FOR  FURTliER  STUDY 
A methodology  for  the  prediction  of  propulsive  Lift  system  jet  Induced 
Aerodynamic  effects  for  multi-jet  VTOL  aircraft  hovering  In  and  out  of  ground 
effect  has  been  dovnlopcd,  To  implement  the  anaLysls,  the  viscous  (turbulent) 
fluwfleld  below  the  aircraft  In  ground  effect  has  been  divided  Into  five  major 
elements i 

1.  Free- jet  development, 

^ 2,  Jet  Impingement, 

3,  Wall-Jet  development, 

>'4.  Stagnation  and  fountain  base  lino  formation. 

S.  Fountain  development  and  impact  on  aircraft  under-aurfoce, 

The  forcea  (auokdown)  and  moments  on  the  fuaelage/wing  Induced  by  the 
entrainment  of  the  turbulent  free-Jots  and  wall-JetH  both  in  and  out  of  ground 
effect  (where  appropriate)  are  computed  by  a Douglas  Neumann  potential  flow 
procedure  which  has  been  modified  apeclflcally  for  VTOL  hover  computations. 

The  entrainment  of  the  Jets  is  modeled  by  praacribingi  from  empirical  data, 
a sink  diatrlbutlon  on  the  Jet  surfacos  which  simulates  thu  inflow  due  to 
turbulent  Jet  entrainment.  TVila  requires  the  panelling  of  pertinent  free-jet 
and  wall-jet  aurfacea.  Forces  and  moments  from  the  induced  potential  flows  are 
determined  on  airframe  aurface  panels  from  local  surface  pressures  computed  by 
application  of  the  Bernoulli  equation. 

The  contribution  to  the  overall  forces  and  moments  on  the  airframe  due 
to  fountain  impingement  la  computed  eeparately  and  added  to  the  induced  potential 
flow  computation  results,  Since  the  fountain  Impinges  on  the  airframe  directly, 
the  prediction  methods  developed  are  mure  empirical  in  nature.  To  calculate 
the  fountain  forces,  an  analogy  between  fountains  and  radial  jets  with  virtual 
origins  located  below  Che  ground  plane  surface  was  (’si.'.ibl  ishod . The  dynamic 
preesure  distribution  in  these  Jets  was  used  to  calculate  the  fountain  impingement 
force  on  the  air  frame,  and  the  ontralnment  distribution  on  these  jets  was  used 
to  calculate  a fountain  suckdown  force, 

Liniitod  comparluons  of  the  results  of  tlio  complete  pvtullction  methodology 
with  expurlmunt/i I data  in  ground  effect  for  a 1.13%  scale  model  of  a four  .lift 
jet  VTOL  atlicls  aircraft  wore  made.  These  comparisons  showod  acceptable  agree- 
ment between  menuured  and  computed  suckdown  for<'es  and  overprcilLctton  of  tlie 
magnitude  of  thn  fountain  impingement  force,  The  overpredict I on  of  tlu>  foun- 
tain ImpinKument  force  is  thouglit  to  be  diie  to  the  Inability  of  the  mcthoilology 
to  predict  fountain  suppresHinn  due  r.o  purtlal  merging  of  tlie  lift  jets  |)flor 
to  ground  tmpln.gement  or  overlap  of  implngc'meril  roglouH  on  the  ground  piano. 
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8 . I Fcom  Jot  Uavelopmunt 

The  data  obtained  on  non-cirenlar  .jeta  show  a c;onrtl.derable  amount  of 
Hoattur  and  are  rather  limited  Jn  seope,  Hpeci  f i e.a  1 1 y , more  inrormatJon 
on  tlie  entrainment  rate  diatrlhutlon  Tor  roctaiiKU lar  .)etH  v«;lth  non-uni lovm 
exit  veleelty  prorilen  Is  n.‘t|uired,  and  information  on  the  entrainment 
eh.iraeter  I at  Les  of  eompreaalhlo  reetanRular  Jets  la  needed, 

8 . i Jet  Imp  1 ne.eiiient  and  W.i  I I -.lets 

The  methodology  deveioped  herein  assumes  that  wail.-Jet  properties  are 
established  by  the  impingement  of  individual  jets  on  the  ground  plane. 

Plow  propurtles  at  tlie  exit  of  the  impingement  region  (ns  investigated  in 
P.eferenoe  2 for  circular  Jets)  arc  needed  for  Jets  from  non-olrcular  noz- 
ziuB,  especially  high  aspect  ratio  (AR  > A)  rectangular  nozzles.  Addi- 
tionally, in  cases  whore  two  or  more  Jets  merge  or  partially  merge  before 
Impingement,  or  where  two  Jets  Impinge  with  overlapping  implngenent  region 
rndll,  R^,  the  formation  of  a fountain  between  the  two  Jets  may  be  sup- 
pressed . The  magnitudu  of  this  suppression  and  Its  depondnneo  on  Jet 
spacing  and  nozzle  exit  taiape  must  be  investigated. 

^ Kraintain  Hovelopiiient 

Thu  fountain  upwauh  ilowfield  models  devehiped  in  this  study  have  been 
restricted  to  fountains  generated  by  a combination  of  circular  or  nearly 
circular  Jets,  all  oi  wltich  Impinge  at  90“  to  the  ground  pliine,  are  located 
at  the  same  lielght  above  the  ground  plane,  mid  have  equal  thrust.  Thero 
nre  no  data  on  the  fountains  generated  when  the  Jet  interaction  flowflcld 
Is  asymmetric  even  for  tlic  simplest  case  of  two  perpendicular  Jets  of  dif- 
ferent thrust.  The  fountain  flowfiuld  development  for  non-circuluc  nozzles 
is  another  area  tliat  requires  further  invustigat lun,  In  the  calculation  of 
fountain  induced  forces  and  moments,  there  are  areas  which  require  further 
development,  even  within  the  restriction  of  the  preseiu.  fountain  flowfiuld 
models.  The  fountain  suckdown  force  Is  currently  calculated  hv  a rather 
crude  method,  and  this  slinrtcomlng  could  be  improved  by  Including  the  founlaln 
entrainment  dlstrihiitinn  In  the  overall  Hiickdown  calculation  performed  with 
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APPENDIX  T 

AIRFHAME  AND  JET  COORDINATES 

Probjtmi:  To  clofino  airframu  f^eomotry,  Jot  poKltlonH,  jol  voi-torJtig 

anfiU'H,  iincl  Jol  Implnpemc'nl  polnfH  for  an  arbitrary  VTOi.  con  f I Rural  Ion 
hovcrini’  In  ground  of  fort. 

Approach!  Derive  transformation  from  airplane  fixed  to  ground  fixed 
coordinates.  Locate  Jet  exit  positions  in  ground  fixed  coordinates  and  Jet 
inclination  angles  relative  to  ground  plane.  For  each  nozzle,  locate  effec- 
tive Jet  impingement  point  on  ground.  Thia  point  is  defined  as  the  inter- 
section with  the  ground  plane  of  a straight  line  incllnnd  at  the  Jot  inclina- 
tion angle  and  going  t.hrough  the  nozzle  center  of  thrust. 

Purpose ! To  provide  basic  geometric  input  to  (static)  aero-propulsive 
Interaction  fjroblem.  Input  geometric  parametuts  are  assumed  to  be! 

0 Alrpliuu!  c.g.  locution  (fusclugu  station,  waUirllno,  Initt  line). 

0 Location  of  center  of  thrust  for  each  nozzle. 

o Thrust  vectoring  and  splay  angles  for  each  nozzle, 

0 Airframe  angle  of  pitch,  roll  and  yaw.  (The  latter  is  only  important 
for  cross  wind  cases.) 

0 Height  of  airplane  c.g,  ahovo  ground  plane,. 

Transformation  Equations 

Transform  unit  base  vectors 


a.  j.  k) 

along  airplane-fixed  axes  (x,  y,  z)  in  Figure  1-1,  to  vectors 

(t".  .f",  k") 

along  iine.s  (A,  B,  C)  respectively. 

The  reference  .angles  arc  the  .standard  Euler  angles! 

H : pitch 
ijj  : yaw 
<li  : roll 


!‘j 


i 
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i 

i 
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FIGURE  1-1 
EULER  ANGLES 
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Defining  flrMt  the  Intermediate  vectors  and  along  lines  OR  and  OS 
reupectlvuly I 

0^  • 1"  cos  ij'  + I"  flln  iji 
0^  “ “I"  sin  i|j  + ,1"  cos  il' 

Using  thcac  vectors,  the  relation  between  the  uniJrlmed  and  doubJ e-pr lined 
base  vectors  may  be  shown  to  be: 

f - {"  [cos  i|i  cos  + 1"  I Gin  I)'  cos  0]  - it"  [sin  0] 

J O'  i"  [cob  i|’  sin  0 sin  1'  - sin  i|)  cos  li] 

+ [I"  [sin  Bln  0 ain  <ti  + cos  i|i  cos  i|i]  + it"  [sin  i|>  cos  n] 

it  P'  [cob  sin  coe  1'  + sin  i|-  sin  'lO  + 

+ [sin  i(i  ain  h cos  '|i  - cos  i|i  aln  i|>]  + It"  [cob  o cos  ll 
Now  Introduce  a third,  ground  fixed,  coordinate  system  such  that 


1'  - i" 

r - 1" 

k'  - -k" 


x'  - x" 

Y ! y n 

- (z"  - h ) 
eg 


m V&xy  • n 


where  h^^  Is  the  height  of  the  center  of  gravity  above  the  ground  plane | posi- 
tive upward. 

In  the  airplane  cncirdinatoa,  the  nth  jet  noKftle  location  will  be  defined 
by  a vector 


■^^jn  ‘ "in  ^jn  "ja 

Kquations  (l-i),  (I.-2)  and  (1-3)  may  be  used  to  obtain  the  nostale  coordinates 
in  the  ground  fixed  ayateni 


Jn 


•r- 


jn 


' y ' . + k ' 

Jn 


,1n 


where 


“ Xj^[cos  iji  COB  0]  + yj^[coa  lii  ain  0 aln  i|)  - sin  ij;  cos  ji]  (I-Aa) 

+ Zj^^[cti«  s'l  sin  0 cos  i|)  + sin  i)j  sin  i|)] 

y'jj^  ■-x^jjsin  iji  coa  0)  - y^^tsin  ij’  aln  0 nin  i|)  4 cos  cos  i))]  (T-Ab) 
- Zj^(Bin  ^ aln  0 cos  (|i  - cos  i|/  sin  ili] 

I 

* yj  isib  i}i  cos  0]  - z.  [cos  0 cos  i|)]  (J.-4c) 

jn  eg  jn  jn  Jn 

The  thrust  vector  direction  for  the  nth  nozzle  is  defined  by  two  angles, 

Thrust  dsflotLlon  angle  A 

n 

Nozzle  splay  angle 

Referring  to  Figure  t-2,  a unit  thrust  vutor  in  tlu>  alrplmu!  coordinate  sysLuiii 
is  then  given  by 

t„  " r (cos  i"  cos  A_)  + J (sin  t ) It  (cos  i:  sin  A ) (1-5) 

n n n n n u 

In  the  pvound  coordinate  uystom,  the  unit  thrust  vector  for  the  nth 
nozzle  is  denoted  by 

t’  • ■?'  t'.^,  4^  1'  t'  , + it''  t'  , 

n ^ n n 

and  using  F.quatlons  (I-l),  (1-2)  and  (1-3),  Equation  (1-5)  yields 


t'  , «~(coB  i,  QOS  A )fuoH  i|)  coa  0]  '■  sin  r;  [coa  i sin  0 sin  i - sin  iji  cos  (j>] 
« n n n 


f (cos  aln  A^)[coh  i|^  sin  0 cos  ji  + sin  I sin  iji] 


(T-6a) 
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FIGURE  1-2 

NOZZLE  THRUST  DEFLECTION  ANGLE  AND  SPLAY  ANGLE 


t'yi  «-(coB  e;^  con  fi^^)[aln  <|)  cos  0]  ~ sin  e^[sln  iJj  sin  6 sin  (|>  + cos  i|i  cos  (j)] 


+ (cos  sin  6^)  [sin  \|i  sin  (■)  cos  -ti  - cos  i|i  sin  i}il 

t'  , - (con  r cos  6^)  [sin  0]  - sin  r fsin  i cos  0'| 

« ^ n n w ^ 

+ (cos  c_  sin  ^ )[cos  o cos  xl 
n n 

The  effective  Jet  lmpir,Roment  point  is  defined  by  the  intersection  of  u 
straight  line  emanating  from  the  nozzle  center  and  parallel  to  the  vector, 
and  the  ground  plana  z'  "■  0.  From  the  equation  for  the  line 

x'-x'  y'-y’.  z'-b'. 

„J.n  . ' _ljn  . Jn 


(l-6b) 


(l“6c) 


it  can  be  shown  that  the  effective  Implngemenc  point  coordinates  are  given  by! 


(I  •■7a) 


(l-7b) 
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Al’PENDlX  II 

wIET  PATH  nEFt.ECTlON  DUE  TO  CROSS-FLOW 


The  piiHltion  and  direction  of  tlic  thrii.sl'  vector  lor  the  nth  nozKle  of  an 
aircraft  ut  un  arbitrary  attitude  may  be  Mpecifiud  by  the  trunafurmutlon 
eciuations  In  Appendix  1.  The  reference  coordinate  syHteni  is  fixed  to  the 
ground  and  defined  by  the  (x'l  y'*  z*)  axes  Hhown  in  Figure  11-1.  The  position 
of  the  thrust  vector  for  the  ntli  jet  Is  denoted  by  * 


r 


Jn 


x' , + J'  y ' , + !<' 

jn  ^ Jn 


(11-1) 


and  its  direction  by  the  unit  vector  • 

t'  - t'  t'  , + J'  t'  , + it'  t’  , (n-2) 

« x'n  y'n  z'n 

In  the  above  equations  (!',  J',  k')  are  unit  vectors  In  the$t',  y'l  z')  direc- 
tions, respectively.  In  u crossflow,  the  Jet  path  for  the  nth  Jet  is  defined 
by  the  empirical  equation  presented  in  Reference  12  t 


f.  ■ 


6 

+ (:  tan 

°n 


(II-3) 


where  (,  and  are  distances  which  have  been  normalized  by  the  jet  diameter, 
and  €Iq^  is  the  complement  of  the  angle  between  the  Jet  velocity  vector  and  the 
free  stream  velocity  vector.  This  equation  has  been  used  with  good  results  in 
studies  of  power-induced  effects  in  transition  flight,  (Reference  13). 

Equation  (II-3)  is  considered  to  apply  in  the  plane  defined  by  the  free  stream 
and  Jet  velocity  vectors.  By  definition  of  the  ground-fixed  coordinate  system. 
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and  since  the  Jet  velocity  vector  direction  is  opposite  to  that  of  tliu  thrusi 
vector 


(Tl-5) 


U9 


The  croea  product  of  these  two  vectors  defines  the  normal  to  the  plane  which 
contains  them,  iind  the  dot  product  dofinas  tlie  unnlt!  hatween  the  free  stream 
and  Jet  velocities.  Thus 
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where  the  unit  vector  e,  and  the  angle  Oq  are  shown  in  Figure  ll-l. 


FIGURE  n-1 

JET  PATH  PLANE  DEFINITION 


Using  Equations  and  (11-5)  in  Equations  (M-6),  it  l'«  posslhle  to 

show  that 

Oq  - tf  . ] (II-7A) 
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and 
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From  Equation  (ll-7b)  It  Is  then  possible  to  calculate  the  angle 
<f>on  “ ^ 

Thla  ang.lfi  ia  niao  ahown  In  FlRuro  11-1. 

In  the  (x',  y',  ?,')  coordinate  Hyatcm,  pointa  along  the  jot  path  are 
given  by 

>‘'Pn  - -‘‘Jn^  >‘'jn 

y'pn  " y'jn  “ '’Jn^  *°n 

*'Pn  " *';ln  “ '^°n 

where  C,  and  C are  related  by  Equation  (tl-l)  and  d^^  is  the  exit  diameter  for 
the  nth  Jet.  Intersection  with  the  ground  plane  is  defined  by  the  condition 

«'pn  - 0 

Using  this  in  (II-9c)  and  substituting  In  (tl-9b) 
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and  using  (11-3)  and  (II-9a) 
■''"n*  ''Jn 


Equations  (11-10)  and  (11-11 
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(11-11) 


define  the  point  at  which  the  Jet  path  intercepts 
the  ground.  To  define  the  angle  of  Intsrsectlon,  Equation  (II-3)  is  differen- 
tiated to  calculate 
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(11-12) 


where  Equation  (ll-9c)  for  - 0 haa  been  us.)d,  and  where  0„  is  the  angle 

between  the  Jot  path  and  a perpendicular  in  the  (/..  O plane . as  shown  in 
Figure  11-1.  With  0„  calculated  from  Equation  (1.1-12),  It  la  then  posal.ble  to 
define  a unit  vector  tangent  to  the  Jet  path,  and  the  components  of  this  vector 
in  the  (x* , y'l  s')  system.  The  result  is 

t'o^  " '^n^  "n  j'On^  *^n  (11-13) 
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As  in  Appendix  I the  effective  impingement  angle  for  the  Jet  la 
On  - cos"^  tk'  • t'o^J 
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In  additloni  the  angle  of  the  olane  in  which  the  above  angle  lleN  in  given  by 
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